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NOTE THE USE NON-ORTHOGONAL WAVE FUNCTIONS 
PERTURBATION CALCULATIONS! 


ABSTRACT 


method first approximation proposed for handling collision problems 
involving chemical ionic change, that the interaction energy function 
different the beginning and end the collision. 


There are many transition problems mechanics for which the 
initial and final states belong two different sets orthogonal states, that 
the calculation one has deal with mixture wave functions belonging 
two different orthogonal sets. then becomes necessary take into account 
the lack orthogonality these wave functions order get adequate 
basis for the calculation. 

example may consider collision between two atoms resulting 
the transfer electron from one atom the other. For the initial state 
have neutral atom and positively ionized atom Bt, each moving with 
definite momentum, and for the final state have the atoms and each 
again having definite momentum, such that the total momentum conserved. 

now have consider the total energy from two points view, 


where H(A), are the energies the atoms Bt,... alone, 
the interaction energy between and and the interaction energy 
between and Initially the state the system eigenstate H(A) 
and for the beginning the interaction the perturbing energy. 
But the final state eigenstate and towards the end 
the interaction the perturbing energy V2. 

The wave equation 


where number, the total energy the state are interested in. The 
initial wave function satisfies 


received August 25, 1956. 
from St. College, University Cambridge, Cambridge, England. 
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Thus, putting 


have 


Let denote the eigenfunctions H(A)+H(B*), representing the atoms 
and each stationary state and each having definite momentum, 
such that the total momentum has the required value. Thus 


where number. Different eigenfunctions are distinguished different 
values for Let denote similarly the eigenfunctions 
that 


with also number. 

The set eigenfunctions far from being complete set with the given 
total momentum, since does not include any wave functions representing 
states with free electron. Even the two sets joined together will prob- 
ably not complete, that arbitrary wave function cannot expanded 
terms them. can, however, expand any wave function, such 


where the residual term, made small possible, which requires that 
shall orthogonal every and every x,. 
Substituting [6] get 
The work far exact. now try get approximate solution [7]. 
discard altogether, that [7] may written 
neglecting the perturbation applied get 
from [4]. Similarly, neglecting the perturbation applied x,, get 


from [5]. Thus approximate solution [7] 


The solution [10] represents ingoing well outgoing waves. can ob- 
tain from solution representing only outgoing waves the standard way, 
replacing the reciprocal and 
similarly for The various terms the solution will then 
the various scattered waves that are produced the collision process, and 


NON-ORTHOGONAL WAVE FUNCTIONS 


their coefficients and will give us, when take the squares their moduli 
and them the appropriate factors, the transition probabilities for 


the various final states. 
calculate the coefficients and multiply [6] and integrate 
over all values the variables that occur the wave functions, thus 


Assuming that the are normalized, that 


this reduces 


Similarly, multiplying [6] and integrating, get 


The equations [11], [12] may solved for and convenient 
look upon these equations matrix equation. Introduce the square matrix 


0 Dry 
and consider the numbers forming single-column matrix and the 


Then [11] and [12] may written 


The solution thus 


essential for the success the method that the matrix shall have 
reciprocal. The condition for this that there shall not exist set numbers 
not all zero, such that 


such that 


This means that the wave functions must all independent. this 
condition does not hold, the expansion [6] not unique and the coefficients 
are not determinate. 

There does not seem any simple way proving the independence 
the examples practical importance. question for further study. 
one should find particular example that the are not independent, 
one might still able save the method working with smaller set 
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functions which are independent. One could get smaller set, for 
example, discarding those that refer one the atoms being 
highly excited state. 

reduction the number that one works with has the advantage 
not only increasing the prospect for them independent, but also 
simplifying the calculation, since reduces the number rows and columns 
the matrix whose reciprocal has worked out. has the dis- 
advantage that the accuracy the method probably diminished, since 


makes larger. 


The foregoing work arose from discussion with Dr. and others 
the National Research Council, Ottawa. 


RADIO OBSERVATIONS THE IONOSPHERE OBLIQUE 
INCIDENCE! 


ABSTRACT 


description given automatic method using the radio pulse tech- 
nique oblique incidence investigate the structure the ionosphere. Sweep 
frequency photographic records were taken for the path 

2355 km.). The records obtained confirm the validity the theory relating 
vertical incidence oblique incidence propagation. Sweeps showing such 
phenomena layer transmission the Pedersen ray, multiple reflec- 
tions, are shown. found that the layer the important propagating layer 
rather than the layer under certain conditions. Hence conventional methods 
predicting frequencies over long distances require modification these 
circumstances. 


INTRODUCTION 


During the last thirty years most high frequency radio observations 
the ionosphere have been carried out normal incidence. The data ob- 
tained from such observations have been applied the calculation maxi- 
mum usable frequencies for radio communication purposes and some measure- 
ments have been made check the relationship between vertical incidence 
and oblique incidence data. appears that, general, the theories are correct 
first approximation. Since most oblique pulse measurements have been 
made short distances (less than 1000 km.) was decided make investiga- 
tions greater distances and especially high latitudes, where the iono- 
spheric characteristics differ from those low latitudes. 

The present soundings, made over distance 2355 km., are part 
comprehensive investigation oblique incidence ionospheric radio propaga- 
tion Canada (2). Owing the nature much the country difficult 
obtain vertical incidence ionospheric observations many points high 
latitudes and oblique incidence soundings offer means investigation 
the ionosphere over these relatively inaccessible regions. 


EXPERIMENTAL PROCEDURE 


The equipment used each station was slightly modified LG17 Automatic 
Ionospheric Recorder (manufactured the Marconi Company Canada). 
block diagram, illustrating the principle operation this recorder, 
reproduced Fig. variable frequency oscillator serves tune both the 
transmitter and receiver. variable frequency Mc.) the range 30.6 Mc. 
about Mc. mixed with fixed frequency 30.0 Mc. and the difference 
frequency (f—30 Mc.) amplified and used control and excite the trans- 
mitter the desired frequency. The received signal again mixed with the 
frequency Mc. and the difference frequency f—(f—30) Mc. Mc. 
amplified and mixed with frequency 28.6 Mc. yield second intermediate 


received July 1956. 
Contribution from the Radio Physics Laboratory, Defence Research Board, Ottawa. Work 


carried out under Project D48-95-11-02. 
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ANTENNA 
FIXED FREQUENCY 
OSCILLATOR 


TRANSMITTER 


FROM ANTENNA 


FREQ. 


EQUEN 

SWEEP 


Fic. diagram LG17 Automatic Ionospheric Recorder. 


frequency 1.4 Mc. which amplified, detected, and displayed cathode 
ray tubes the conventional manner. 

this experiment pulses duration were transmitted the rate 
per second, which was also the repetition frequency the time base. 
All the timing features (repetition frequency, time base, height markers, etc.) 
were obtained counting down from 100 kc. crystal standard having 
stability few parts 10°. enable the pulses from the distant trans- 
mitter set convenient positions both time bases, variable frequency 
oscillator Ottawa could switched permit adjustment the equip- 
ment triggering time. The frequency this oscillator could varied continu- 
ously from about kc. about 101 kc. 

keep the frequencies both equipments exactly the same the variable 
frequency oscillator the LG17 was replaced oscillator which the 
frequency was increased steps about 100 kc. ensure very accurate 
frequency control, auxiliary circuits were included which ‘pulled’ the frequency 
that the nearest harmonic the 100 kc. frequency standard. descrip- 
tion these auxiliary circuits will given elsewhere. The oscillator remained 
any single frequency for period two seconds and then ‘stepped’ quickly 
the next frequency. The period for one complete sweep the LG17 was 
320 sec. that frequency sweep 16.0 Mc. could obtained. The frequency 
ranges generally used were: 


2.0 Mc. 18.0 Mc. night, and 
5.0 Mc. 21.0 Mc. day. 


this way the high daytime and the low nighttime maximum usable fre- 


= | 
| 
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quencies could covered. Sweeps were taken once every min. The trans- 
mitting and receiving antennas were horizontal rhombics. The sweeps were 
recorded photographically mm. film; each photographic sweep record 
was either in. length. 
The computation the equivalent path signal between sender and 
receiver can made follows. Referring Fig. 
let the interval (expressed kilometers, where km. 3.3 usec.) 
between the triggering the equipment and the reception 


the interval between the triggering the equipment and the 
the equivalent path between and 
the equivalent path between and 
Then 


STATION 


STATION 


SIGNAL 


Fic. Timing cycle. 


Assume that the equivalent paths are equal the two directions. Then 
[2] Pas 


assist obtaining comparable records both ends the path, and 
were kept nearly equal. practice, since approximately 2500 
the Ottawa equipment was adjusted that the signals received from Sas- 
katoon were within the range 2200 km. about 2800 km. 


TYPICAL SWEEPS 


Experiments were conducted the summer and December 1954. Some 
typical sweeps are shown Figs. Fig. was taken near noon and shows 


| Pa | 

! 
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BA 


CANADIAN JOURNAL PHYSICS. VOL. 


716 


SS 5 


Fic. typical winter day equivalent path versus frequency sweep. Saskatoon, Decem- 


ber, 1954, 13:21 E.S.T. 


Fic. Winter day sweep showing layer cutoff phenomenon. Ottawa, December, 1954, 


14:01 E.S.T. 


Fic. Sweep taken during slightly disturbed period. Saskatoon, December, 1954, 17:21 


E.S.T. 


PATH 
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ARBITRARY 
ORIGIN 


800 
600 
MEGACYCLES 


Fic. (a) Sweep showing the influence sporadic and relatively strong Pedersen ray. 
Saskatoon, December, 1954, 17:01 E.S.T. 
(b) Vertical incidence sweep near mid-point path. Note the 
presence echoes. Kakabeka Falls, December, 1954, 17:05 E.S.T. 


reflections from the layer with both the once-retlected signal and the twice- 
reflected signal present. The once-reflected signal extends from Mc. 
Mc. and, its high frequency end, shows the more strongly retarded 
Pedersen high angle ray. The twice-reflected signal not transmitted 
frequencies above 13.5 Mc. owing the acute angle incidence the 
ionosphere. The changes equivalent paths near the lower frequency limits 
both traces are probably due the laver winter, not sepa- 
rated from the layer. may noted that the double-hop trace more 
spread near the maximum frequency than the corresponding portion the 
one-hop trace. This feature appears true general. 

somewhat similar sweep shown Fig. Here there interesting 
feature, namely, the retardation the layer the once-reflected echo 
near its low frequency limit. The lower cutoff frequency about Mc. agrees 
quite well with the values calculated from the vertical incidence critical 
frequency near the mid-point the circuit. 

more complicated sweep that shown Fig. difficult determine 
the maximum usable frequency (MUF) with accuracy because weak signal 
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Sweep with reflection from the region. Saskatoon, December, 
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MEGACYCLES 


Fic. Typical September day sweep. Ottawa, September, 1954, 13:2 


Fic. Sweep taken September evening. Note the close resemblance winter sweeps. 


tember, 1954, 21:01 E.S.T 


ARBITRARY ORIGIN 


t 
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present frequencies above the so-called ‘classical MUF’. This phenomenon 
was noticed frequently and was found that the maximum observable frequency 
depended, some extent, the characteristics the equipment. much 
spread two-hop reflection also present. intermediate echo can seen 
about 2290 km. between 8.0 and 9.0 Mc. 

Another sweep considerable interest shown Fig. 6(a). will noticed 
that the single-hop Pedersen ray* has appreciable strength frequencies 
well below the maximum frequency. Between 5.0 Mc. and 8.0 Mc. there 
little sign the once-reflected signal whereas the multiple reflected signals are 
quite strong. addition the one-hop, two-hop, and three-hop reflections, 
there are several traces which are produced reflections from the layer 
followed reflections from the layer. Evidence that this the case 
provided vertical incidence sweeps obtained near the mid-point the 
circuit (see Fig. which show type reflections, i.e., echo received 
after reflection from the layer followed reflection the ground and 
finally from the layer. 

Fig. shows one the relatively rare winter sweeps when reflections were 
obtained below the main echo. was found that such echoes were only 
obtained when sporadic ionization was present heights greater than about 
130 km. 

typical September day sweep shown Fig. This sweep typical 
summertime conditions and more complicated than most winter sweeps. 
The maximum frequency observed, about Mc., reflected from the layer. 
The signal received after two reflections from the layer can also seen, 
having maximum frequency near Mc. will also noticed that the 
maximum frequency reflected from the layer greater than the maximum 
frequency reflected from the layer. The importance the Pedersen rays 
relative the low angle rays clearly evident Fig. Note also the splitting 
the layer Pedersen ray into ordinary and extraordinary components 
near its lower frequency limit. 

record taken September evening shown Fig. The and 
layers had disappeared, thus giving the sweep the appearance winter’s 
day sweep. 

RECIPROCITY 

has not yet been adequately demonstrated that propagation radio 
signals the ionosphere completely reciprocal. Indeed, the influence the 
earth’s magnetic field may such that the principle reciprocity may not 
always hold. Hence this point worthy investigation. 

section was assumed that reciprocity was valid, i.e., 

(i) that the transmission times are the same both directions, and 

(ii) that the sweeps obtained both ends the circuit are identical. 


Unfortunately, not possible make direct test (i) since there 


*This ray, also known the high angle ray, has smaller angles incidence the reflecting 
layer than the low angle ray and reflected from regions higher electron density. the frequency 
increased these two paths come closer together and are identical the maximum frequency 
transmission. 
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nie i 
2600 
2200 
MEGACYCLES 


Fic. 10. and (6) Simultaneous sweeps taken Ottawa and Saskatoon illustrate 
the principle reciprocity. December, 1954, 10:41 E.S.T. (a) Saskatoon, Ottawa. 


common time reference. However, taking the low angle one-hop trace 
reference possible compare the time delays the one-hop Pedersen ray 
and the multi-hop modes transmission. Figs. and show that, 
this respect, the patterns Ottawa and Saskatoon are very similar. 


OBSERVED AND CALCULATED MAXIMUM USABLE FREQUENCIES 


From the point view practical radio communication the maximum 
frequencies for any circuit are great importance. Various methods 
developed for relating vertical incidence h’f data oblique trans- 
mission problems. For present purposes the theory developed Millington 
(3) has been used convert vertical sweeps oblique sweeps. Beynon (1) 
England has shown that over distance 715 km. there fairly close corre- 
spondence between vertical and oblique sweeps. Sultzer and Ferguson (4) 
have verified this for distance 1150 km. 

Owing the inaccessibility the country the mid-point the circuit 
was decided locate ionospheric recorder Kakabeka Falls, Ontario, 
situated about 200 km. south the mid-point. 

Examination the records obtained Kakabeka Falls showed that, near 
midday, the minimum height reflection the layer was between about 


PATH 
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2600 2600 
2400 
2200 220 
4 
2000 2000 
MEGACYCLES 
(a) DECEMBER 1954 12.20 EST. (b) DECEMBER E.S.T. 
——CALCULATED 


Fic. 11. (a) and (6) Comparison between observed sweeps and sweeps calculated from 
vertical incidence data. 
190 km. and about 210 km. Since the true height reflection must less 
than this, value for 180 km. was used calculate the transmission 
curve, where the height the bottom the layer. Two examples 
oblique sweeps calculated from vertical sweep data are shown Figs. 
and The observed curves were obtained averaging the values for 
Saskatoon and for Ottawa (see Section 2). Considering that Kakabeka Falls 
about 200 km. from the mid-point together with the simplifying assumptions 
that are used the theory, the agreement between deduced and observed 
curves reasonable. The maximum frequencies are usually fair agreement, 
but the observed path delays appear consistently greater than those 
calculated. 

24 


CALCULATED MUFs 

© OBSERVED MUFs 
----MINIMUM OBSERVED FREQ, 
VISUAL OBSERVATIONS 


——CALCULATED MUFs 

© OBSERVED MUFs 
----MINIMUM OBSERVED FREQ. 

VISUAL OBSERVATIONS 


(a) (b) 


MEGACYCLES 
B 
MEGACYCLES 
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Fic. 12. (a) and Comparison between observed and calculated maximum frequencies. 
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Observed and calculated MUFs are shown Figs. 12(a) and will 
seen that, the whole, the agreement reasonably good, again bearing 
mind the limitations involved. should noted that the maximum fre- 
quencies generally were scaled from the photographic records. However, the 
frequencies which the signals disappeared were noted the operator 
during the sweep and the latter values are included Figs. and 
for cases. where photographic records could not scaled. plotting the 
minimum frequency which echoes were observed, the diurnal variation 
the useful frequency spectrum emphasized. During the day, communications 
could established frequencies between about 6.0 Mc. and 20.0 Mc. 
whereas night the usable frequency spectrum compressed between about 

has been noted that the station Kakabeka Falls about 200 km. south 
the mid-point the path. comparing the values 
obtained Ottawa, Kakabeka Falls, and Winnipeg was found that, 
the average, the difference between the values obtained stations sepa- 
rated 1180 km. was 0.38 Mc. and distance 570 km. the difference was 
0.2 Mc. Hence reasonable expect that for separation about 200 km. 
the average difference would about 0.07 Mc. That is, the data from Kaka- 
beka Falls are not significantly different from those which would have been 
obtained the precise mid-point. The differences between the Ottawa data 
and the Winnipeg data are sufficiently great render these data unsuitable 
for determining ionospheric conditions the mid-point. 

Owing the curvature the earth there limit the distance over which 
wave may propagated the single-hop mode transmission. This 
limiting distance depends the height reflection. interest note 
that, over much Canada, summer days the heights reflection the 
layer are about 200 km. and winter days the heights reflection the 
layer are about 210 km. Now for height reflection 210 km. the 
limiting distance for single-hop transmission about 3300 km., whereas 
problems concerning the prediction maximum usable frequencies, 
customary assume that one-hop transmission effective 4000 km. 
This implies height reflection about 320 km., value appropriate 
equatorial, but not polar, regions. 
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MATHEMATICAL ANALYSIS SOLUTE REDISTRIBUTION 
DURING SOLIDIFICATION! 


ABSTRACT 


mathematical analysis made the redistribution solute which occurs 
during the solidification, without convection, single-phase alloy. Calculations 
solute distributions are carried out for steady state conditions, and for the 
transient changes solute concentration, both the solid and the liquid, which 
occur the beginning and end solidification and result change the 
speed solidification. The solute distributions the transient regions are con- 
sidered for the case two solutes present the melt. shown that, under 
certain conditions, junctions may produced the transient regions soli- 
dification germanium silicon melts containing two suitable solutes. 


INTRODUCTION 


well known that, melt impure metal uniform composition 
solidified under non-equilibrium conditions, redistribution the impurity 
occurs, that the resulting solid not uniform composition; that is, segre- 
gation occurs during solidification. The problem calculating the impurity 
distribution which will obtained solidifying single-phase alloy under 
specified conditions has been treated Pfann (6) and others under the 
assumption that mixing complete all points the melt all times during 
solidification. This problem has also been treated Tiller a/. (10) under the 
alternative assumption that there convective mixing; rather, assumed 
that the only mixing due diffusion the impurity the melt. Wagner (11) 
has given analysis the problem, taking into consideration both convection 
and diffusion the liquid. has calculated the effect convection 
segregation, and has shown that, for wide range solidification rates 
commonly used the freezing castings and the growth single crystals, 
natural convection has little effect the segregation. These calculations have 
recently received some support the work Walton al. (12) who have 
studied experimentally the conditions solidification corresponding the 
interface the type studied Rutter and Chalmers (8). The conditions 
solidification corresponding this transition were predicted from the theoreti- 
cal considerations Tiller (10), which are based the assumption 
convection the liquid. The fact that the experiments were agreement with 
this theory indicates that convection was absent during solidification, thus 
substantiating the work Wagner. 

Since natural convection appears play small part the segregation over 
wide range practical conditions solidification, is, therefore, useful 
treat the problem considerable detail under the assumption that convection 
negligible. Accordingly, the present paper constitutes extension the 
work Tiller giving more precise and more complete calculation the 

received August 11, 
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segregation which will obtained under various conditions solidification. 
The results will plotted for various conditions solidification practical 


interest. 
MATHEMATICAL ANALYSIS 


Assumptions 

order treat the problem quantitatively, necessary employ more 
specific assumptions than those given the introduction. The problem which 
will considered that the end-to-end segregation which accompanies the 
solidification impure metal melt which uniform cross-section (Fig. 1). 


Interface 


Direction Interface Movement 


Fic. solidification melt having uniform cross-section. 


The case which single-phase alloy formed will considered. The follow- 
ing assumptions will used: 

Diffusion the solid negligible. 

Convective mixing the liquid negligible. 

The distribution coefficient, defined the ratio the concentration 
impurity the solid being formed the concentration the liquid 
from which forms, constant. Attention will directed the case 
less than unity. 


The interface separating solid and liquid plane and perpendicular the 
axis the specimen. 


The rate diffusion impurity the liquid specifically taken into account. 


The Equation 


Under the assumptions made, growth proceeds the movement 
interface separating solid and liquid. This interface considered plane 
and perpendicular the axis the specimen. The problem here clearly 
linear one, since growth takes place one direction only. 

Since the concentration impurity appearing the solid related that 
the liquid adjacent the interface the factor the distribution the 
solid resulting from solidification may obtained calculating the concen- 
tration the liquid the interface for all positions occupied the interface 
during freezing. has been assumed that the only factor causing movement 
solute the liquid diffusion and, therefore, the equation which must 
solved the one-dimensional diffusion equation. 

With regard fixed co-ordinate system, this equation may written as: 
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where the coefficient diffusion the solute the melt, 
the concentration solute the melt, 
distance, 
time. 

However, the solute distribution relative the moving inter- 
face required and is, therefore, convenient use co-ordinate system 
which moves with the interface. Transformed into this co-ordinate system, 
equation [1] becomes: 


[2] 


where rate movement the interface, assumed constant, 
distance measured from the interface into the liquid. 


solution for equation [2], giving the distribution solute both the 
liquid and the solid, will obtained for semi-infinite sample the follow- 
ing cases: 

The initial transient rise concentration solute which occurs the 

first part the sample freeze. 

The steady state which reached after the solidification has proceeded 

for some distance from the beginning. 

The terminal transient change concentration which occurs the last 

portion the sample freeze. 

The transient change concentration which occurs result change 

the speed solidification. 


The Steady State Solution 

For the case distribution coefficient which less than unity, the solid 
which forms will purer than the liquid from which freezes and, result, 
rejection solute into the liquid will occur the interface. 
Owing the finite rate diffusion impurity the liquid, layer high 
impurity concentration will build adjacent the advancing interface. 
the concentration rises the liquid the interface, must also rise the 
solid being produced, since the two concentrations are related the distribu- 
tion coefficient the original melt uniform concentration then, 
clearly, the layer high impurity concentration must build until the con- 
centration appearing the solid also the concentration the liquid 
the interface will then When this condition has been attained, the 
layer high concentration will have reached steady state, and further 
variation with time will take place long the conditions growth are not 
altered. Since time has been eliminated variable, the steady state distribu- 
tion impurity the liquid may found from solution the equation: 


subject the boundary conditions: 
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atx’ =0. 
The solution has been obtained Tiller (10) and may written: 


[4] 


where (1—). 

natural variable this equation the dimensionless quantity 
Equation [4] plotted the form vs. for various values the 
parameter Fig. will noted that the concentration falls off from the 
interface into the liquid more rapidly for higher rates growth and that the 
amount impurity carried ahead the interface, obtained integration 
equation [4], greater for lower rates growth. 


— 


D 


Steady state solute distribution the liquid ahead the advancing 
interface. 


The Initial Transient 

order calculate the way which the region high solute concentra- 
tion builds adjacent the interface during the first part the solidification 
and the resulting transient rise concentration the steady state value, Co, 
the solid produced, necessary solve the time-dependent equation [2]. 
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The solution will carried out for the case semi-infinite sample which 
solidification initiated one end time and proceeds constant 


speed, 
Equation [2] must solved subject the following boundary conditions: 


The boundary condition given equation [7] expresses the fact that the 
concentration must remain finite the interface all times. The 
solution equation [2] which satisfies the above constraints will describe the 
distribution solute the liquid all times. 

the use Laplace transform, equation [2] may transformed from 
partial differential equation and into ordinary differential equation 
x’. This accomplished replacing by: 


0 
Substitution [8] equation [2] gives: 


dx” 


Rearranging terms, get 


Equation [9] may now solved for and operating with the inverse 
Laplace transformation, the time dependent solution, restored. 
The solution carried out Appendix and gives the concentration, 
throughout the liquid any time, 

The impurity concentration appearing the solid equal times the 
concentration the liquid any time. Therefore, the impurity 
distribution resulting from solidification the initial transient region may 
the co-ordinate Rt, which represents distance measured from the begin- 
ning the sample. derived Appendix the concentration, the 
solid, given by: 


Again, natural variable the solution the dimensionless quantity 
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(R/D)x. plotted against (R/D)x for various values the pa- 


Fic. Initial transient solute distribution the solid. 


Sudden Change the Rate Solidification, 

Suppose that the steady state distribution, solute the liquid has 
been established growth rate The rate growth then abruptly 
changed Since the amount solute carried ahead the interface 
the steady state case smaller for than for there will transient 
rise the concentration the liquid near the interface, and hence the solid 
formed, the adjustment the new steady state made (10). time 
measured from the instant which the speed change occurs, this transient 
may found solving equation [2] subject the following boundary 
conditions: 


This problem may solved the use Laplace transformation, 
manner similar that employed for the initial transient case. The solution 
derived Appendix II. The concentration distribution produced the solid 
result the speed change given by: 


where distance from the point which the speed change occurred. 
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the case the initial transient distribution, this equation may put 
dimensionless form plotting against the indepen- 
dent variable. However, since equation [11] contains two parameters, and 
obtained for each value Three such families are shown Figs. 
and 4(c). 

The Terminal Transient 

Let the advancing interface, carrying ahead the steady 
state impurity distribution, approach the end the specimen constant 
speed, The solute represented the area under the curve must appear 
this last part solidify, producing transient rise concentration. The 
impurity distribution this terminal transient portion the sample may 
calculated the following manner. 

Take the origin the end the specimen, and let represent distance 
measured from this origin. Let the interface distance from the origin. 
This co-ordinate system fixed with respect the specimen, rather than 
moving with the interface and, therefore, the differential equation satis- 

order that the concentration shall remain finite the interface and the end 
the specimen during solidification the terminal transient region, the 
following boundary conditions must satisfied: 


Concentration 


The first these boundary conditions may met extending the sample 
the left, shown Fig. and assuming that two interfaces exist, one 
right with the same speed. clear that for this the same 
terminal transient function will obtained for the speci- 
men’’. The solute distribution ahead each interface known the steady 
state distribution until the interfaces get close enough interfere with the 
second boundary condition. The steady state distributions are given by: 


=> — — 
I3R j2R Ry RI 
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and 


order simplify the considerations, these distributions may regarded 
waves strength moving speed the constant concentration 
the two waves approach each other, time will come when they seriously 
disturb the second boundary condition. This may corrected placing 
the left and right respectively with velocity 2R, indicated Fig. The 
double speed keeps the spacings equal and also balances the double distance 
the original sources. later stage, the sources interfere, and neces- 
sary add sources +3y, strength moving speed This 
argument may repeated indefinitely and ultimately yields infinite 
sequence sources for which 


outlined Appendix III, the coefficients can determined such that 
the second boundary condition, equation [14], satisfied; the first boundary 
condition, equation [13], satisfied for any values C,. The analysis given 
Appendix III results expression for the concentration solute the 
liquid any point for any interface position. considering the concentration 
the interface, was done the cases the initial transient and the change 
speed, the concentration produced the solid result the freezing 
the terminal transient region may obtained. given by: 


less than unity, this series diverges infinity for For all other 
values the series converges, and does quite rapidly for larger values 
Xo. 

the series converges even for though perhaps quite slowly. 

the case the initial transient, equation [18] contains the dimension- 
less quantity natural variable. plotted against 
for various values the parameter Fig. 


DISCUSSION 


The Extent the Transient Regions Solidification 

frequently desirable produce, solidification from the melt, 
specimen which uniform composition. this accomplished, all 
transient regions solidification must avoided subsequently removed 
from the solid sample. Consequently, interest know how large the 
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Fic. Terminal transient solute distribution the solid. 


transient regions solidification will be. basis for comparison the 
lengths the transients, the distance, measured from the beginning (or end, 
the case may be) the transient region, which half the transient change 
has occurred may used. For example, the case the initial transient, the 
area between the curve and the line represents the 
solute deficiency the transient region. The value (R/D)x corresponding 
one-half this area will used measure the extent the distribu- 
tion. number such are given Table 


TABLE 
COMPARISONS EXTENTS SOLUTE DISTRIBUTIONS 


Distribution Steady state, Initial transient, Terminal transient, 
coefficient (R/D)x’ (R/D)x 
0.5 0.693 1.4 0.8 
0.1 0.693 7.6 <0.01 
0.01 0.693 <0.01 
0.001 0.693 670 


will noted that the half-area value for the steady state distribution 
equation [4]. The half-area values for the initial transient were obtained 
graphically. the case the terminal transient, precise computation half- 
area values not feasible for values smaller than about 0.5 because the 
slow rate convergence equation [18]. 

interest compare the precise equation for the distribution the 
initial transient region, equation [10], with the approximate equation given 
Tiller (10). Their equation for the initial transient distribution may 
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—k(R/D)z 
written: Cs(x)/Co 


This equation plotted Fig. along with equation [10] for comparison. 
will seen that the equation is, fact, quite accurate; the maxi- 
mum error appears about per cent. Hence, the approximate equation 
may used with good accuracy. 


Equation Tiller etal. 


Equation 


Fic. equation [10] for initial transient distribution with that given 
Tiller al. 


Two methods may used for the production specimen uniform com- 
position solidification. mixing the liquid complete all times, the 
process zone levelling (6) may used. there convection the 
liquid (or constant amount convection) growth under steady state condi- 
tions outlined this paper may employed. However, this latter method 
chosen, care must taken avoid the formation the prismatic substruc- 
ture described Rutter and Chalmers (8), since the formation this sub- 
structure accompanied segregation impurities into the boundaries 
the elements the structure. practice, this means that the steady state 
method restricted cases distribution coefficients close unity low 
solute concentrations (10). 


The Nature the Transient Regions Solidification When Two More Solutes 
Are Present the Melt 

two more solutes are present low concentrations melt, may 
assumed that they will behave independently. this the case, the theory 
presented this paper may applied each solute individually order 
calculate the distribution resulting from solidification under given set 
conditions. For simplicity, the discussion will limited the case two 
solutes. 

Let two solutes, having different distribution coefficients, distributed 
uniformly throughout the melt, the concentration one being greater than 
that the other. will shown below, possible, during transient 
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regions solidification, produce solid which the major solute the one 
which was present the smaller extent the original melt. the parent 
material the melt germanium and one the solutes will produce n-type 
germanium while the other produces p-type germanium, then such reversal 
the predominant impurity the solid will result change conductivity 
type with the consequent production junction. Because the general 
interest the production junctions, the effect will discussed this 
context. 


(1) The Initial Transient 


Suppose that two solutes, which shall numbered and having different 
distribution coefficients, are added melt pure germanium. Let: 


will assumed that greater than 

proper choice concentrations, possible have less than 
and the same time Co(1) greater than this the case, 
impurity will appear the solid greater concentration than impurity 
during the first part the initial transient region solidification; impurity 
will excess throughout the rest the initial transient well the steady 
state which follows. the impurities are, for example, gallium 0.12) (7), 
which produces p-type germanium, and antimony 0.007) (7), which 
produces germanium, then the sample will p-type the first portion 
the initial transient and thereafter. junction will formed 
the sample the point where the atomic concentrations antimony and 
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Fic. Production junctions crossover initial transients for two solutes (gallium 
and antimony). 
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gallium are equal. The distributions antimony and gallium, calculated from 
equation [10], the portion the initial transient region interest here are 
shown Fig. The antimony and gallium concentrations are plotted 
fraction Co(2), the initial concentration antimony the melt, vs. 
for various values the concentration ratio Fig. shows, both 
the thickness the p-layer and the difference between the concentration 
gradients the two solutes the junction depend upon the rate solidifica- 
tion and the starting concentrations. result, the rate change the 
difference concentration, which may called the differential concentration 
gradient, across the junction may controlled suitable choice the rate 
solidification and the starting concentrations the melt. The difference 
between the antimony and gallium concentrations, obtained subtracting 
corresponding values Fig. plotted Fig. the slope this curve 
the junction gives the differential concentration gradient across the 
junction. 


So 
> 


‘ 


Fic. Concentration difference across junctions Fig. 


number junctions may produced one sample the following 
method. Suppose that solidification has passed the initial p-region, described 
above, that the solid being formed germanium. Distributions 
gallium and antimony will have built the liquid adjacent the interface 
the extent corresponding the stage solidification which has been 
reached. these solute distributions the liquid are now destroyed, for ex- 
ample stopping solidification and stirring the liquid, the first material 
solidify when growth continued will p-type germanium. That is, the 
initial transient region solidification will repeated and p-region will 
formed between two indicated Fig. 10. The junction formed 
the point where growth resumes will step junction, while the other 
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Point which solidification 
begins after mixing 


Distance along specimen 


Fic. 10. Production distribution repetition initial transient. 


junction, formed the manner described above, will graded junction. 
This process could repeated several times sample sufficient length, 
thereby producing series junctions. The thickness the p-layers 
produced this method could controlled varying the speed solidifica- 
tion and the starting concentrations employed. Theoretically, there reason 
why the p-layer could not made thin enough produce sample from which 
transistor could made. producing p-layer requisite thickness for 
this purpose, the possibility choosing solutes having more suitable distri- 
bution coefficients must considered. The requirement here have the 
point equal atomic concentrations the two solutes occur close 
possible the point which solidification resumed. clear that this 
will favored rise the concentration impurity from Co(2) 
which rapid possible and rise impurity from Co(1) which 
slow possible. the basis the approximate equation Tiller al., 
may shown that the maximum rate rise will occur for distribution 
coefficient 1/2. Therefore, production thin p-layer will favored 
the use two impurities, one which has distribution coefficient 1/2 
and the other which has distribution coefficient greater than unity. the 
case germanium, the closest available approach this situation appears 
(1) p-type impurity. 

example p-layer produced the method described above sample 
germanium containing antimony and gallium impurities shown Fig. 
11. This sample was solidified horizontal graphite boat approximate 
the conditions the theory. The step junction appears line running across 
the etched sample. The location the graded junction not revealed 
etching, but can found using hot copper probe. The operation the 
probe based the fact that the thermal e.m.f. generated between copper 
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graded step 
junction junction 


Fic. 11. germanium crystal containing antimony and gallium. 


and n-type germanium opposite polarity that generated between copper 
and p-type germanium. The width the p-layer contained the sample 
Fig. approximately 0.005 in.; was measured using the hot probe 
mounted micromanipulator. Similar samples have also been produced 
this technique using boron and arsenic impurities germanium. 

Owing, first, lack accurate information concerning such quantities 
the diffusion coefficients, the distribution coefficients, and the concentrations 
the solutes and, second, the impossibility producing the ideal experi- 
mental conditions required the theory, has not been possible obtain 
accurate comparison theory and experiment. The theoretical values for 
p-layer width are larger than but, the average, within one order magni- 
tude those obtained experimentally. 

Layers of-p-type material between regions material have also been 
produced germanium Hall (4, and silicon Tanenbaum al. (9). 
The method used these workers has been called the growing” tech- 
nique. consists essentially varying the crystal pulling rate the melt 
temperature both while growing crystal from properly doped melt. 
These workers have attributed the production regions different conduc- 
tivity types their samples variation the distribution coefficients the 
impurities with speed solidification, the amount variation depending upon 
the particular solute. This explanation based the experiments Hall 


ball 
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(4, 5), who found such variation germanium crystals pulled from melt 
which was stirred rotation the crystal 200 revolutions per minute. 
was believed that this would produce complete mixing the melt all times. 
However, doubt has been cast upon the validity this explanation the 
work Burton (1, 2), which shows that Hall probably did not have 
complete stirring the melt during his experiments. 

addition the method crystal growth employed, the technique used 
these workers produce p-layers differs other respects from that 
described the present paper. Hall (4, and Tanenbaum al. (9) used 
compensated (9) which were stirred during growth, whereas 
the present technique employs melt which not compensated and which 
precautions are taken avoid convective mixing during growth. 


(2) Change the Rate Solidification 

shown Fig. the extent the transient change concentration pro- 
duced given change speed depends upon the distribution coefficient 
the solute concerned; smaller distribution coefficient results larger 
increase concentration for given change speed. Therefore, possible, 
when solidifying melt containing two solutes which have different distribu- 
tion coefficients such antimony and gallium germanium, change the 
major impurity appearing the solid from one the other over limited 
region the transient produced change speed. This illustrated 
Fig. 12, where impurity has 0.5, and impurity has 0.1; the curves 
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Impurity 
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12. Production distribution change rate transients two suitable 
solutes. 


are plotted for 0.2. the case the change rate transient, both the 
junctions produced germanium sample would graded junctions. 
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APPENDIX 
The differential equation solved the transformed equation 


derived from the original equation, [2], the substitution 


0 


The boundary conditions, equations [5] and [7], transform 


and 

Try solution the form 

where and are constants. Substitution [22] the equation [19] gives 
_1R 

using the second condition, [21], found 


Substitution the values [23], and [24], equation [22] gives 


The time-dependent solution may now obtained transforming back 
the required through the inverse Laplace transform. The required 
transform has been given Campbell and Foster (3). Incorporating their 
solution, the solute concentration any point the liquid for any time 


Since the convergence domain the transform violated this problem, 
necessary check the results the computations. That equation [26] 
the correct solution the problem has been verified showing that 
satisfies the following conditions: 


= 
i 
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(4) may shown direct differentiation and substitution 
that 
(5) direct differentiation and substitution, may shown that the 
solution satisfies the original differential equation. 


The concentration the solid may found putting equation 
[26] and multiplying putting and rearranging terms, the con- 
centration the solid may expressed function distance, x,.from the 
beginning the specimen, the following form: 


The. following additional checks the correctness this solution have 
been verified: 
(3) The quantity solute carried ahead the interface must equal the 
solute deficiency below the initial transient region the solid. 
Therefore, the following equation must hold: 


0 0 
direct integration, has been shown that this relation satisfied. 


APPENDIX 


this case, the differential equation satisfied the concentration, 


The boundary conditions are: 


Appendix the equation may transformed using 


0 


Substituting [28] yields the transformed equation 


2 
72 
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The integral this equation 
[34] 


Therefore, try solution 


Applying the third boundary condition, equation [31], form the quantity 


x r’=0 


This gives 

Solving equation [36] for yields 

[38] 


Substituting [37] and [38] [35] gives the form 


will seen that the first and third terms this expression are, except for 
place the initial transient given Appendix The second and fourth 
terms give rise new transients. applying the appropriate inverse Laplace 
transforms, given Campbell and Foster (3), the terms equation 
the time-dependent solution may obtained. The fourth term gives 


q 


744 CANADIAN JOURNAL PHYSICS. VOL. 


and, general, 

n—1 n—1 
Therefore, substituting [48] [44], 


n=co n—1 n—1 

0 n=1 0 0 

may easily shown that each term the above series satisfies the differen- 

tial equation, noting that dy/dt Hence, the above expression satisfies all 


the requirements the problem. 

Putting gives the concentration the interface for any interface 
position Multiplying gives the concentration the solid the position 
from the end the specimen. Replacing agree with previous 
notation: 


Combining terms [50] which have the same exponentials: 
n=cp n—1 n—1 


n=co 


Written out full, [51] becomes: 
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C2 k(1+k) 


and, general, 

n—1 n—1 
Therefore, substituting [48] [44], 


may easily shown that each term the above series satisfies the differen- 
tial equation, noting that Hence, the above expression satisfies all 


the requirements the problem. 

Putting gives the concentration the interface for any interface 
position Multiplying gives the concentration the solid the position 
from the end the specimen. Replacing agree with previous 
notation: 


Combining terms [50] which have the same exponentials: 


Written out full, [51] becomes: 
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COSMIC-RAY EXPERIMENTS WITH PROTON VELOCITY 
SELECTOR! 


OGILVIE? 


ABSTRACT 


Experiments with selector slow heavy cosmic-ray particles are described. 
The pressure coefficient protonsat sea level the momentum region 700 1100 
Mev./c found 0.64+0.08% per mb. Points are obtained the differential 
momentum spectrum protons sea level between 650 and 1250 Mev./c, 
agreement with the results Mylroi and Wilson. This spectrum discussed 
connection with the experiments Filthiith between 200 and 900 Mev./c. 


INTRODUCTION 


Duerden and Hyams (8) 1952 described apparatus designed study 
the heavy particles the cosmic-ray beam. This utilizes the fact that, under 
the appropriate conditions, liquid Cerenkov counter will detect with very 
high efficiency particles passing through with velocity greater than the 
Cerenkov critical velocity. This critical velocity c/n, where the velocity 
light free space, and the refractive index the medium. Fig. shows 


Fic. the principle the velocity selector. 


the arrangement principle. The tank containing the liquid used the 
Cerenkov radiator placed between the two trays Geiger counters and 
Light, generated the passage fast particles through the tank, 
detected the photomultiplier and gives rise output pulse from 
coincidence with the pulses from the Geiger counters. anticoincidence 
the type denotes the passage particle moving with less than 
the critical velocity. can arranged that mesons, electrons, etc., having 
great enough range penetrate the apparatus, will have velocities greater 
than the critical, that the recording anticoincidence denotes the 
passage heavy particle. This apparatus thus selects particles able 
penetrate it, but unable emit Cerenkov light. The majority such particles 
the cosmic-ray beam sea level are protons. 

their paper, Duerden and Hyams show that the velocity meson 
greater than the Cerenkov velocity water its residual range water 
greater than about gm./cm.? They also show that the effect knock-on 
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electrons produced water the passage protons moving with velocity 
below the critical velocity negligible. anticoincidence 
with the arrangement Fig. thus indicates the passage heavy particle. 
shall call the range equivalent the apparatus and the residual 
range particle velocity c/n. is, thus, particles with ranges between 


DESCRIPTION APPARATUS 


practice, necessary collect sufficient fraction the light, 
due the passage fast particle, the photocathode the phototube. 
also necessary sure that light emitted fluorescence, excited 
particles moving slowly through the tank. The refractive index the 
liquid must low enough include between the limits and sufficiently 
wide band energies. These three requirements place stringent conditions 
the substances used Cerenkov radiators. One such suitable liquid 
distilled water, the energy protons corresponding the range being then 
480 Mev. 

The absorption coefficient water for light the spectral region the 
maximum response the photocathode such that 10% absorbed 
path length meters (10). Studies the spectral distribution the Ceren- 
kov light excited water the gamma rays from indicate that fluo- 
rescence negligible (Rich and Slovaceck, 16). this apparatus, the tank 
containing the water made brass, and two E.M.I. type 6262 phototubes 
are mounted watertight glands, similar those described Harrison (9), 
the center opposite sides. The photocathodes these tubes are the end- 
window type, and also cover strip in. wide around the periphery adjacent 
the end. The tubes are therefore mounted project in. into the tank. 

The output pulses from the phototubes are applied conventional 
amplifier, gain approximately 80, which they are connected parallel. 
the output this amplifier the pulses have rise time and 
duration approximately The inside surfaces the tank are coated 
with white waterproof paint, developed the paint research laboratories 
the National Research Council (Harris and Ogilvie (8)). Between the 
wavelengths 4250A and 7000A, the reflection coefficient this paint 94%, 
the pigment used being zinc oxide. 

Examination those pulses from the phototubes that are coincidence 
with the Geiger counters shows that they have considerable spread ampli- 
tude, the average pulse corresponding the emission about electrons 
from the cathode, Fig. Since meson, moving with velocity emits 
250 photons per cm., this indicates collection efficiency per tube, 
compared Duerden and Hyams value 15% for tank having volume 
45% this one and containing one multiplier, and agreeing closely with 
calculations based the theory Mando (12). the present experiment two 
tubes were found necessary, and the combined gain each, together 
with the following amplifier, sufficient deliver the coincidence circuit 
pulse from one photoelectron that will operate it. 
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COUNTING RATE- ARBITRARY UNITS 


PULSE SIZE ARBITRARY UNITS 


Fic. height spectrum pulses from the phototubes coincidence with triple 
Geiger counter coincidences. 


block diagram the apparatus given Fig. the form which 
was used study the heavy particle beam. The circuits used were for the most 
part conventional, the coincidence mixer between the photo and Geiger 
tube circuits, where large spread input pulse sizes encountered, being 
the type described Garwin (5). The Geiger counter trays were in. 
square, and the extreme ones were separated distance in. 


COINCIDENCE 


IGER T 
GEIGER COUNTER DOUBLE 
COINCIDENCE 


GEIGER COUNTER 
CERENKOV 
TANK 
ABSORBER ———— sc. ALER— PROTONS 


SCALER—QUADRUPLE 


Fic. diagram the experimental arrangement. 
ACCIDENTAL ANTICOINCIDENCES 


anticoincidence could occur because accidental triple coincidence, 
leading too high proton rate. Accidental quadruple coincidences can 
neglected since they cannot recorded protons. The expression for the 
number accidental coincidences between three counters 


where the resolving time the coincidence circuit and and 
represent the real coincidence rates between pairs counters. Owing 
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the high efficiency the Cerenkov counter for the detection mesons, 
the contributions the third and fourth terms are negligible, since each 
case particle must cross the tank, leading quadruple coincidence. The 
term likewise very small, with triple coincidence circuit resolving 
time Experiments without the lead absorber below the tank (see 
below) show that relativistic particle must cross only 20% the depth 
the tank recorded. The second term thus the only one interest 
us, and has been shown (15) small. further check this the 
proton anticoincidence counting rate has been measured while the resolving 
time was varied, and the results indicate that accidental triple coincidences 
account for less than the proton counting rate where the resolving time 
used the measurements described here. They will, therefore, 
neglected below. 
RESULTS MEASUREMENTS 

With the in. thick layer lead below the Cerenkov tank, which neces- 
sary exclude slow mesons, the apparatus records particles with mass 
greater than 550 electron masses and, particular, protons the energy 
range 230-480 Mev. was decided study the nature the anticoincidence 
counts both measuring their barometric pressure coefficient and com- 
paring the flux with the proton flux this energy range, measured 
Mylroi and Wilson (14). 


THE BAROMETER COEFFICIENT 


The apparatus was run for period several months Ottawa, (altitude 
320 ft.), and the results linear correlation counting rate with pressure 
are given Table below. 


TABLE 
BAROMETER COEFFICIENT DETERMINATION 
No. readings (hourly) 940 348 476 
Average pressure (mb.) 1003 .54 
Counting rate per hour 47.1 47.15 
Barometer coefficient, %/mb. 


The weighted mean barometer coefficient per mb. 


identical piece equipment was also used aircraft altitudes 
25,000 ft., Ogilvie (15), and the absorption length determined. The result 
was 128+12 agreement with the absorption length found 
Tinlot (19) for the variation penetrating shower rate with atmospheric 
depth. The results are set out Table II, together with the pressure coefficient 
and corresponding absorption length obtained Simpson (18) for the slow 
neutron component. 

will noted that the pressure coefficients for protons obtained the 
two different experiments are not equal. Such small discrepancy could prob- 
ably explained change the angular distribution the protons with 
altitude. 
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TABLE 
ABSORPTION LENGTHS 


Method Pressure coefficient Absorption length, gm./cm.? 
Slow neutrons (<700 gm./cm.?) (18) 140+3 
Protons, sea level 0.64+0.08 158+20 
Protons, aircraft (15) 0.78+0.06 128+10 


(calculated) 


suppose the whole the discrepancy between the sea level slow 
neutron and proton pressure coefficients due the admixture some 
slow mesons with the latter, taking the meson pressure coefficient 
(Rose, 17), the minimum proportion protons 85%. 
plot some results obtained during one the barometer coefficient experi- 
ments given Fig. showing daily mean values and the slope the best-fit 


line. 


PRESSURE mb. 


COUNTS PER HOUR 


Fic. Some results pressure coefficient experiments. Daily mean results are plotted. 


THE PROTON FLUX 


The proton spectrum, obtained Mylroi and Wilson (14), has been re- 
drawn differential energy spectrum Ballam and Lichtenstein (1). 
the region interest, 230 480 Mev. (700-1100 Mev./c), the proton 

the present experiments the observed anticoincidence rate was 46.9 per 
hour, when corrected barometric pressure 1000 mb. This value the 
mean the results obtained during the pressure coefficient experiments, and 
during the same period the mean meson, triple coincidence, rate was 16,100 
per hour. Thus, taking the absolute intensity the vertical hard component, 
with cm. water and 3.8 cm. lead the telescope, par- 
obtain uncorrected value for the proton 


4 
1020 
1000 
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for the scattering loss between 230 and 350 Mev., detailed Ogilvie (15), 
amounting (19+5)%, arrive the value particles/ 
cm.?/sec./steradian, agreement with Mylroi and Wilson (14). view 
the foregoing, may conclude that the apparatus counts heavy particles, 
that excludes slow mesons with high efficiency, and that large proportion 
the particles counted are protons. 

estimate for the proton flux the energy region 150 230 Mev. can 
derived, and the same time estimate the minimum distance that 
relativistic particle must move the tank order detected can 
made. 

The minimum energy required proton pass through the telescope 
with lead absorber being 150 Mev., order find the number 
protons recorded the range 150-480 Mev. therefore require the counting 
rate with lead absorber together with estimate the number mesons 
which produce anticoincidences when lead present below the tank. 
the anticoincidence counting rate the apparatus with the 3.8 cm. 
lead absorber below the tank, and with above the tank, then 


number protons the energy range 480 540 Mev. 
number such protons lost the absorber owing interactions 


Here the number mesons counted protons where the absorber 
the upper position, namely, those which are slowed the absorber 
travelling with 0.75 for large fraction their path through the tank. 
The number protons the range 480 540 Mev., Mylroi and Wilson (14), 
accounts for approximately seven counts per hour. The observed counting 
rates with the lead different positions are shown Table 


TABLE III 
COUNTING RATES FOR DIFFERENT ABSORBER POSITIONS 
Position absorber Counting rate, 
per hour 
Above tank 54.5 
Below tank 41.6 
None 98.3 


Thus counts per hour, and assuming that the number 
protons lost the 3.8 cm. small, per hour approximately. 
Noting that the differential range spectrum low energy mesons sea level 
obtained Lichtenstein (11) closely constant this region, obtain the 
number slow mesons recorded with absorber, need only multiply 
the ratio the range intervals the two cases. this way arrive 
counting rate protons per hour the energy range 150 540 Mev., 
and thus per hour the region 150 230 Mev., particles/ 

From the value and the meson spectrum referred above, can 
deduce that the highest energy with which undetected meson enters the 
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water from the absorber about Mev. 0.8), and that traverses one 
quarter the tank before going below the Cerenkov threshold. emits 
the average about 120 photons per cm. this path, total 600 photons. 
Thus relativistic particle, emitting 250 photons per cm., must only traverse 
15-20% the depth the tank, cm., order detected. 


ABSORPTION CURVE 


After these experiments had been performed, increasing quantities lead 
were placed above the bottom Geiger-counter tray, which had been moved 
about in. below the bottom the tank. The uncorrected column Table 
shows the counting rates then obtained. When there large gap between 
the tank and the lower tray, the error due accidental triple coincidences 
and showers increases, and the correction for them was found the following 
way. The position the bottom tray, and thus the solid angle, remained the 
same during the whole the measurements the absorption curve. Thus 
the anticoincidence rate this position, should equal 
where the proton rate with the same quantity, 3.8 cm., lead and the 
bottom tray immediately under the tank, and and are the triple coinci- 
dence counting rates. Actually where respresents the 
accidental and shower correction. this way was found that 12.5 
counts/hour. This will remain substantially the same lead added above 
the lower tray. Table the counting rates for various thicknesses lead 


TABLE 
ABSORPTION MEASUREMENTS 

Inches Uncorrected Corrected 
lead absorber rate, rate, 

6.5 
6.6+1.5 
5.5 
20.4+1 
4.5 20.2+1 
22.3+1 

3.5 23.8+1 

25.3+1 12.8+1 

2.5 30.0+1 17.5+1 

24.6+1 

1.5 31.1+1 


are given. The corrected rates, obtained subtracting from the observed 
rates, are shown the table and plotted Fig. The rate 31.1 counts per 
hour, conjunction with meson rate 11,147 per hour, corresponds 
proton flux 3.0+0.1 

The range protons 480 Mev. lead 225 gm./cm.? Thus, with total 
absorber 7.8 in. lead equivalent comprising the telescope, tank, and lead 
absorber (6.5 in.), protons should counted. The residual rate, 5.6+2 
counts/hour, equal 15% the proton flux, provisionally assumed 
represent upper limit for the flux deuterons the appropriate energy 
range, 460 Mev. 960 Mev. 
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COUNTS PER HOUR 


Fic. Absorption curve for protons obtained with the apparatus. 


The above procedure for obtaining neglects scattering, which the 
same the two positions the bottom tray, since the 3.8 cm. lead absorber 
moved with it. also assumes that the difference between the angular depen- 
dence protons and mesons does not introduce appreciable error. set 
measurements the angular dependence protons, made with apparatus 
the present type, gives cos"w law with approximately equal 


THE PROTON FLUX HIGHER ENERGIES 


order obtain information about the proton spectrum energies 
above 480 Mev., increasing quantities lead were placed above the apparatus. 
The usual 1.5 in. lead was place below the tank. The results these 
experiments are summarized Table 


TABLE 
COUNTING RATES WITH LEAD ABOVE THE APPARATUS 


Inches Energy band, Mev. Counting rate, Intensity, 


Lower Mean Upper sec. 
280 415 550 34.0+1 0.92 
300 437 575 32.0+1 0.915 
320 459 597 29.2+0.9 0.85 
345 477 610 27.7+0.9 0.88 
370 500 630 0.81 
395 522 650 0.82 
415 542 670 21.9+0.8 0.83 
440 565 690 19.5+0.8 0.80 
465 582 700 0.81 


N.B.—Intensity corrected the basis the absorption length being 300 
(Bernardini al. (2) and George and Jason (6)). 


480 MEV. 
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THE DIFFERENTIAL MOMENTUM SPECTRUM PROTONS SEA LEVEL 


Several authors have examined the proton spectrum sea level and the 
results obtained are shown Table and Fig. The results 1952 


TABLE 
THE PROTON MOMENTUM SPECTRUM 


Momentum 
interval, 
No. Observer ev./c X108 Date 
Rochester and Bound* 0.34 -0.40 1940 
Ballam 0.42 -0.91 
0.51 -0.96 +14 
Merkle, Goldwasser, 0.48 -0.88 1950 
and Brode 
Mylroi and Wilson* 0.53 -0.71 +5.3 1951* 
0.71 10.7+1.9 
0.93 -1.35 6.6+1.1 
York 0.55 -0.75 1952 
0.75 -0.88 10.7+0.8 
0.18 -0.21 2.2+0.3 1954 
0.21 3.7+0.5 
0.275-0.312 
8.0+1 
0.432-0.446 
469 
0.615 10.0+1 
0.730 8.4+1 
0.887 4.8+1 
0.70 -1.10 
0.85 -1.240 6.5+0.5 
0.95 6.0+0.7 


*Corrected York (20). 

for absorption and production secondary protons. This would increase three 
high energy points about 20%. 

all points Table are plotted, interest clarity. 


have been discussed York. Those Ballam; Merkle, Goldwasser, and 
Brode (13); and Goldwasser and Merkle (7) are early measurements. The other 
points divide themselves into two classes, namely, those due Filthiith (4), 
and the others. Both suggest maximum the region 500 600 Mev./c, 
but since the results Filthiith constitute almost the only information 
available below about 500 Mev./c, worth while enquiring the cause 
the discrepancy between the two classes higher energies, amounting 
factor about two. view the good agreement between the results 
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° 100 200 300 =— 400 500 600 700 800 900 1000 00 1200 1300 
MOMENTUM 
The differential momentum spectrum protons. Numbered points refer Table 


Mylroi and Wilson, York, and the present author, shall consider the 

experiment Filthiith see possible systematic error might found. 
The value the proton intensity, this type experiment given 


where the meson intensity, and C,, and are the meson and proton 
counting rates under the same conditions. the energy (momentum) 
range covered the experiment. Thus accurate knowledge 
decisive for the determination the proton intensity, but careful considera- 
tion leads Filthiith the conclusion that his points would not error 
owing this cause more than 10%. Since the whole curve drawn through 
the points marked Fig. appears displaced the direction lower 
intensities, the cause must least approximately independent incoming 
proton momentum. estimate the scattering protons the absorbers 
the experimental arrangement used Filthiith indicates that scattering 
might sufficiently large account for the observed discrepancy. 


CONCLUSIONS 
The momentum spectrum protons sea level considered agree with 
that given Mylroi and Wilson the momentum range between 600 and 


3.0 
2.0 
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1250 Mev./c, when measured means Cerenkov counter velocity selector. 
The spectrum given Filthiith between 200 and 700 Mev./c considered 
displaced the direction low intensity. 
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FORBIDDEN TRANSITIONS DIATOMIC MOLECULES 


CARBON MONOXIDE! 


ABSTRACT 


Two examples forbidden electronic transitions heteronuclear molecules, 
and are studied the spectrum the molecule. The 
bands have been obtained absorption the region 1750 1230 with 
absorbing paths 400 cm-atm., using the fourth and fifth orders meter 
vacuum spectrograph. most the bands all four predicted branches 
are observed. For the bands, somewhat unusual structure predicted: 
two branches and form and three branches form. the best-resolved 
bands, four branches are observed which closely fit the predicted branches, two 
the branches being very nearly coincident and unresolved the present 
spectra. The rotational and vibrational constants the upper states, and 
the two band systems have been determined. Some these data had 
previously been obtained from perturbations other band systems CO. 
the whole, the perturbation data agree satisfactorily with the more precise data 
obtained directly the present study. 


INTRODUCTION 


earlier papers this series (13, 15, 17), examples four types 
forbidden transitions homonuclear molecules were discussed and analyzed. 
Several other types had been discussed previously various authors (for 
summary see (16)). Thus far, only two types forbidden transitions have 
been analyzed for heteronuclear diatomic molecules, viz., and 
The best-known example transition the Cameron bands 
first observed absorption and recognized transition Hopfield 
and Birge (19) and later analyzed detail Geré, Herzberg, and Schmid 
(10), (8), and Rao (22). Several other examples this type are known, 
particularly for heavier molecules for which the selection rule less 
strong. The only example transition the red bands GeH 
analyzed recently Kleman and Werhagen (20). 

investigation the ultraviolet absorption spectrum carbon mon- 
oxide under high resolution, have analyzed examples two further types 
forbidden transitions heteronuclear diatomic molecules, namely 
and 

The absorption spectrum the vacuum ultraviolet was previously 
investigated Leifson (21) and Hopfield and Birge (19). Leifson (21) observed 
only the simple progression the fourth positive group starting 1546 
while Hopfield and Birge, with higher resolution and longer absorbing paths, 
observed number additional band systems including the very much 
weaker Cameron which start 2063 Between the start the 
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Cameron bands and that the fourth positive group, they observed another 
very weak progression which they designated distinguish from the 
Cameron bands a—X. Later, Asundi (1) found the red and infrared the 
transition (Asundi bands). From work these bands and perturbations 
other emission bands the visible and ultraviolet regions, was shown 
various authors (4, 25, 11) that thea’ state state. 

The dispersion used Hopfield and Birge was not sufficient resolve the 
bands nor follow them through the much stronger A—X (fourth 
positive) bands. This was done the present investigation and this way 
absorption band system was found among the much stronger bands the 
fourth positive group. The bands this new system have very peculiar 
structure and have been assigned transition. The upper state 
identical with the state whose existence had already been derived from 
perturbations the state Coster and Brons (4) and Schmid and 
(26).4 Much improved rotational and vibrational constants both the 
and the states have been derived. 


EXPERIMENTAL 


The absorption spectrum was photographed the fourth and fifth 
orders the meter vacuum spectrograph previously described (3). The 
reciprocal dispersion was 0.63 and 0.50 A/mm., respectively. The light source, 
previously, was Lyman discharge through argon mm. pressure. 
The lower orders were eliminated the use LiF prism (5, 3). 

For most spectra absorption tube cm. length was used which was 
filled with pressures ranging from 760 mm. order obtain 
longer path lengths for the weaker bands the forbidden transitions, the 
whole spectrograph was filled with pressure 400 mm., thus reaching 
absorbing path 400 cm-atm. 

few experiments were also carried out with absorption tube cm. 
length entirely cooled liquid nitrogen. diagram this tube shown 
Fig. Considerable difficulty was experienced obtaining tight seal 


Ground gloss 
joint cemented 


Liquié nitrogen with Zeluphen 


J 


“—Lead 
spectrograph j>-Fluorite window Copper To lightsource 
Copper glass seal 


Fic. Low temperature absorption tube. 


between the LiF windows and the absorption tube low temperature because 
the difference the coefficients expansion LiF and the material the 
tube. This difficulty was finally overcome using Cu-glass joints, the 


state should not confused with the state called Fig. 197 Ref. (14). The latter the 
upper state the Kaplan bands and should perhaps renamed 
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part which was fused lead cylinder which turn was sealed Zeluphen 
cement (supplied Zeh and Company, Chemische Fabrik, Ludwigsburg, 
Wiirtt., Germany) the windows. 

the previous work oxygen (3), lines second order served 
wavelength standards. The wavelengths were converted vacuum means 
Edlén’s tables (6) except the case the exposures with the spectro- 
graph. For these exposures, since the index refraction the vacuum 
region not known, correction curve was established from spectral features 
that could measured accurately both with proper absorption tube and with 
the absolute values some wave numbers, but the relative accuracy 
estimated well within 0.1 

The most conspicuous feature the absorption spectrum the 
region 1800 1200 the progression the fourth positive system. 
This progression masks considerable part the region the absorbing 
paths necessary observe the forbidden transitions. addition, the 
progression the fourth positive system which appears with considerable 
intensity overlaps number the forbidden bands. was the object the 
low temperature absorption spectra reduce this overlapping eliminating 
completely the progression the fourth positive system and 
reducing the extent the individual bands the progression. this 
way, some the new bands could recognized much more clearly. However, 
most the measurements were made from the plates taken normal tempera- 
ture since these plates were somewhat better quality. 


The transition forms one the band systems known 
the Hopfield—Birge bands. These bands, mentioned earlier, were discovered 
absorption Hopfield and Birge (19). They assigned the first band 
observed them, 58927 the band the system. the basis 
additional band found Birge 57763 Estey (7) raised the 
vibrational quantum numbers the state one unit. The final assign- 
ment was deduced Beer (2) from rotational perturbations the 
system caused the state. According the perturbation data, the 
vibrational numbers proposed Hopfield and Birge have 
raised three units. The expected position the 0-0 band then 55380 
Some fairly accurate data about the rotational structure the 
state have already been obtained Geré (9) and Geré and (11) 
from study the fine structure few the Asundi bands 
addition, from the numerous perturbations produced the state 
the and states, Schmid and (25) and Beer (2) have 
obtained further information about the constants the state. Beer 

the present investigation, the state reached absorption from 
the ground state single long v’-progression (with single- 
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headed, red-degraded bands observed. starts with very weak band 
55363 which, according the perturbation data, the band. 
The intensity increases fairly rapidly shorter wavelengths, reaches flat 
maximum about the 15-0 band, and then drops again. The last band that 
could observed under the available experimental conditions was the 23-0 
band 78190 The band heads, with their estimated relative intensities, 
are listed Table The absorption intensity the overlapping fourth 
positive system, which allowed electronic transition, approximately 
thousand times higher. 

Fig. greatly enlarged spectrogram the 5-0 band 1632.4 
reproduced. shows clearly the presence four branches, two branches 
and two branches giving the appearance doublet structure. shorter 
wavelengths, because the smaller wave number dispersion, the doublet 
branches (which have the smaller splitting) are longer resolved. 

The observed the bands immediately understood the basis 
the diagram Fig. which indicates the possible transitions the basis 
electric dipole selection rules. According the selection rules AJ=0, 
obtain two branches form (AN +1) and two branches 
form (AN —1), but branches form (AN 0). The designations 
the branches follow from the values and AN. They are (see Fig. 


Fic. band. The quantum number was formerly called 


RR, the two branches having and form respectively. The 
structure the bands entirely similar that the atmospheric oxygen 
bands except that the symmetries are reversed since the oxygen 
bands represent magnetic dipole transition with 

shown Fig. the separation the branches and corresponds 
the separation the sublevels and the state while the separation 
the bands. The order the levels F2, corresponding N+1, 
respectively, normally not that shown for simplicity Fig. but 
lies either below above the two others which, turn, are close together. 
Which case applies here can ascertained from the combination differences: 
only and not give values 
that agree with those known for the ground state. this way, the numbering 
given Fig. was established. 


RQ RQ’ RQ Rg Rg 
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other states, the triplet splitting very nearly constant for 
values above Near the band origin, for the bandsare not well resolved 
and the deviations from constant value, expected from the Schlapp (24) 
formulae, cannot detected. However, Geré and (11) from the analy- 
sis the Asundi bands have established the variation low The sum 
the splittings corresponding members and branches gives the 
constant (see (15)). the present case, this constant found 
for and for 19, that is, unlike the 
case the and states (15, 3), there only very slight, 
any, change with vibrational quantum number. The negative sign 
indicates that, unlike the ground state the and levels are 
above with the highest. 

Table are given the wave numbers the lines all bands that 
could measured with reasonable accuracy. Lines marked asterisk 
are blends are very faint and are therefore not accurately determined 
the others. 

Because the unfavorable wave number dispersion the vacuum region, 
the accuracy with which the rotational constants B’’ and the lower 
state could determined much inferior that obtained from infrared and 
microwave spectra. Therefore, attempt was made evaluate them. 
Rather, the and values were determined from the combination sums 

spectrum (18). For all bands, plot the left-hand side 
the above equation against gave very nearly straight line indicating 
that within the present accuracy (probably the same 
23. The slope the straight line gives and therefore B’; the 
intercept gives The band origins and the rotational constants determined 
this way are included Table The and values can represented 
the following two equations: 

Here the quantity 1081.74 the zero point energy the ground state 
CO. The and values calculated from these equations well the 
values obtained Schmid and (25) and Beer (2) from the perturbation 
data are included Table for comparison with the observed values. 

Even though the state causes many perturbations the 
and states, only one clear perturbation was found it, namely the 

state which has not been observed directly. The reason for the apparent 
absence the other perturbations the bands the fact that these 
perturbations occur mostly high values which were not observed 
absorption. 
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P(J) 


3-0 band, 1697.2 


58912.81* 
907 .62* 


58922.11* 58913.39 
922.11* 908.48 


TABLE 

750.73 752.06 916.63 
718.53 748.43 719.48 912.81* 
728.57* 684.58 730.00 685.53 
720.04 670.74 721.40 671.75 
699.05 700.37 640.37 
673.06 674.40 604.07 835.83* 
658.18 584.04 820.49* 
642.05 561.66 643.40 562.75 
606 .06 515.43 516.46 766 


61191.29* 61183.55* 61192.50* 61184.28 
039.49* 040.50 164.19 187.85 164.89 
030.25* 031.44 182.37 154.89 155.66 
052.81 019.75* 021.06 176.76 144.29 
028.85 980.78 982.02 151.86 104.47 105.26 
018.24 965.12 019.75* 140.82 142.04 089.30 
948.10 007 .68 071.20 129.69 071.94 
929.85 931.12 114.80 052.45 115.96 
978.49 910.22 979.79 911.58 099.72 032.40 100.93 
962.64 889.33 890.64 011.07 011.89 
927.00 28.12 844.91 046.42 964.24 965.12 
736.78 840.93 956.41 854.43 855.37 
813.44 707.69* 707 930.54 931.88 

62297.08* 62289.85 62290.85 63384.15* 
292.06 293.49 270.96 355.63 378.73 356.29 
281.80 250.83 335.16 367.91 335.87 
274.60 275.87 359.18 322.77 360.36 323.50 
265.97 267 225.11 350.27 351.44 309.77 
192.90 245.89 193.90 328.20 329.41 278.01 
231.85 175.22 176.16 314.99 316.27 260.11 
217.64 156.19 218.97 157.11 300.38 301 .66 240.75 
202.08 135.76 203 .08 284.41 218.94 285.61 219.85 
185.20 113.89 186.52 114.90 196.76 268.15 
147.19 148.57 148.04 229.19 148.98 
126.10 040.16 127.48 041.13 121.39 
012.83 105. 013.81 183.10 184.35 094.58 
054.52 955.04 132.78 134.14 034.21 

939.78 819.78 941.06 820.79 

907 908.91 

873.93 875.37 


902. 23* 
885.79* 
837 .04* 
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TABLE (concluded) 


The second system here under discussion has not previously been observed 
even under low resolution. The reason that this system wholly submerged 
the far stronger fourth positive group and can separated from only 
spectra high resolving power. Some the bands can only seen the 
low temperature plates since they are overlapped bands the 
fourth positive group and others are unobservable because overlapping 
the bands this system. 


447.01 448.02 474.41 475.46* 
450.55 405.25 475.46* 431.06 476.65 
368.39 369.48 448.18* 449.51 395.92* 
326.68 391.80 327.73 414.84* 352.11 416.28 
375.32* 376.99 304.64* 
331.49 
305.56 
| 13-0 band, A, = 1439.0 A 16-0 band, A, = 1383.2 A 
0 | 69493 .65* 
494.09* 
490.44 491.96 475.81 278.67 294.69 
486.33 466.19 487 .63* 269.73 
480.63 455.89 481.91* 456.80 
464.03 430.70 465.32 431.48 231.88 
415.42 454.62 416.30 216.53 
441.09 398.81 442.36 239. 198.94 
411.73 360.50 412.99 361.48 158.84 
394.63 338.98 395.90 339.69 136.59 
375.92 377.24 169.81 112.19 
333.60 334.94 124.33 058.42 
043.74 
013.47 014.63 
19-0 band, 1334.5 20-0 band, 1319.7 
748.19 
721.25 722.46 


the and 20-0 bands, the and branches are longer properly resolved. However, 
each the unresolved broad blends has been measured two lines. 


§ 
4 
f 


HERZBERG AND HUGO: FORBIDDEN TRANSITIONS. 765 


Like the system, the system consists single 
progression weak bands which are degraded the red and show single 
heads only. The progression starts with relatively strong band 64816 
which only the head could observed because the rest the band 
masked the band the fourth positive system. Because the relative 
strength this band and the absence further bands this progression 
toward longer wavelengths, was assumed that the band the 
system. sure, longer paths and higher pressures, extremely weak 
feature observed 63686 which might band preceding 64816 
the progression. However, off the expected position and 
was therefore not considered belong it. The progression the 
system could identified with certainty far the 10-0 
band 74709 shorter wavelengths, the bands are either completely 
masked badly overlapped the fourth positive and systems. 
Band heads could, however, recognized the expected positions the 
12-0, 13-0, and 17-0 bands the progression. The band heads, with their 
vibrational assignments and estimated intensities, are listed Table 

casual inspection the fine structure the bands the e—X system, 
exemplified the 3-0 band Fig. gives the appearance 
transition with and branch. turns out, however, that this 
assumption agreement with the well-known combination differences the 
ground state can obtained. Rather the differences 
are about twice what they should be. This, together with the observation that 
the best-resolved bands the branch double, led search for other 
interpretations the fine structure. was found that transition 
accounts all details for the observed band structure. Fig. diagram 
the expected rotational transitions shown. Considering the electric dipole 


Fic. Branches band. The quantum number was formerly called 


selection rules +1), five branches are expected, viz., 
three branches form (AN and and one branch each 
S(AN +2) and O(AN —2) form: and Since state the 
and levels N+1) lie close together and somewhat above below 
the levels N), the and branches will nearly coincide, but will 
separated from the branch. This conclusion agrees with the observation 
the best-resolved band shown Fig. most other bands, the three 
branches are entirely unresolved. 
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readily seen from Fig. that the following combination relations hold 
for bands: 


Table shows these differences obtained from the 3-0 band well the 
values from the infrared bands. Considering the facts that the 


TABLE 


COMPARISON VALUES OBTAINED FROM THE 3-0 BAND AND THE 4-0 
BAND 


From From rot. vib. spectrum 

11.39 11.53 
19.22 19.06 19.22 
27.20 26.91 
34.56 34.50 34.60 
50.02 
57.60 57.69 
65.34 65.14 65.35 
72.99 73.06 
80.61 80.66 
88.25 
96.19 
111.28 111.36 111.36 
119.02 119.11 118.99 
126.77 126.66 
134.35 134.31 
142.05 141.95 


branches and have not been resolved and that the wave number 
accuracy the vacuum region not far good the infrared, the 
agreement shown Table must considered excellent. supplies 
definite proof that the transition really 

Schlapp (23) has given the intensity factors for the branches both 
and bands. For the latter these intensity factors contain, 
addition two constants and which depend the detailed nature 
the electronic states. For the purposes this paper, since quantitative 
intensity measurements have been made, only the intensity factors the three 
form branches are interest. Schlapp gives 

for 
for 
for 


According these formulae, were small compared the branch 
would about four times intense (at large the branch while, 
that the sum and about strong (see Fig. thus not 
disagreement with theory. 


| 


768 CANADIAN JOURNAL PHYSICS. VOL. 


Table the wave numbers the observed lines the bands are listed. 
Blended and very weak lines are again indicated asterisks. 
According Schlapp’s formula (24), only the levels state follow 
the standard formula 


However, this level, according Fig. not involved the and form 
branches, but only the branch. Therefore, the best way determine 
and from the branch which represented 


Plotting against J(J+1), one obtains very nearly straight line 
indicating that D’—D” small, i.e., D’’. The slope the line gives 
and, since known, B,’ obtained. The intercept gives vo. The 
and values obtained this way are given Table 

The magnitude the doubling the form branch, that is, the separation 
between and approximately constant and nearly equal 
The average this quantity from the bands which the branch 
resolved 0.84 order get the correct for the bands which the 
branch not resolved, the quantity 0.42 has added the value 
obtained from the lines. This assumes that the measured position the 
unresolved lines half-way between the true positions and 
that is, that has the same intensity (see above). The correc- 
tion 0.42 has been taken into account for the 1-0, and 7-0 bands 
Table III. 

The B,’ values and the band origins can represented the formulae 


vo = 65340.o5+ 1093.993(v’ + 3) —9.578(v’ +3)? 


where again the quantity 1081.74 the zero point energy the ground 
state CO. The values given the above formulae are included Table III. 

mentioned before, the existence the state has previously 
been derived Coster and Brons (4) and Schmid and (26) from per- 
turbations the Angstrém bands and the fourth positive group. From these 
data, approximate positions the vibrational levels and approximate 
values were derived the authors mentioned. These values are included 
Table III. Considering that these data are often based very few perturbed 
lines, the agreement with the more precise values here obtained quite 
striking, particularly for the This agreement gives support other mole- 
cular constants which have been derived from perturbation data only. 

Apart from the a’—X, e—X, and the strong A—X bands (fourth positive 
group), two other long progressions appear our plates, the a—X bands 
(Cameron system) and the d—X bands, both 
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TABLE 


HERZBERG AND HUGO: FORBIDDEN TRANSITIONS. 


1-0 band, 1517.7 


2-0 band, 1493.7 


2 872.70 890.66 929.13 928 .31* 945. 28* 

3 869.18 890.66  65856.89* 924.92 924.07 946.13 66911.73* 
4 863.62 890.66 846.18 919.16 918.42 945. 28* 901. 19* 
5 856 . 89* 888.53 834.37 911.73* 911.73* 943. 13* 889 .39 
6 848.71* 885.54 821.22 903 . 86 902.93 939 .75* 875.94 
7 839.11 880.91 806.76 894.01 893 . 27 934.92 861.30 
8 828.14 874.99 790 .83 882.81 882.05 928.31* 845.38 
9 815.98 867 .54 773 .68 869 .78* 869 .78* 920.80 — 

10 802.23 858 . 89 755.07 856.14 855.48 911.73* 809 .05 

13 753 .03 24.46 691.26 805 .78* 805 .78* 743.90 
14 734.31 810.26 666 .93 786 . 23 785.30 

15 713.71* 794.65 641.70 764.71* 764.71* 

16 692.01 777 .88 615.67 742.64 741.63 

17 668.99 759.48 586.54 

18 643.58 739.84 557.14 

718.76 

561.70 672.78 

22 531.77 

3-0 band, A, = 1471.0 A 5-0 band, A, = 1428.9 A 

0 67976 .53 

1 67968.84  67968.10 979 .91* 

3 961.80 960.85 982.48 948.88 

4 955.88 955. 28* 981.81* 938.31 69971 .13 

6 940.15 939 .37 975.40 912.89 954 . 26 

7 930.10 929. 26 969.98 897 .96 943.89 

918.54 917.80 963.51 881.71 932.07 

9 905.66 904 . 84 955. 28* 863 .92 918.29 

10 891.38 890.45 945.85 844.81 903 .58 

11 875.55 874.67 934.93 824.21 886.74 

13 839.66 838.74 908 . 55 878.51 849.42 

750.67 838.74 

724.84 723.99 

697.51 
770.94 
745.70 
22 718.76 

7-0 band, 1390.1 9-0 band, 1355.2 

910.35 

767.40* 

8 869 .06 732.99* 732.99* _ 699.75 
855.33 718.51 717.69 
822.22 683.92* 636 
664.11 
13 643. 16* 643. 16* 
14 620.40 619.42 
15 596.01 595.19 
16 569.98 569 .08 
17 442.70 41.27 


systems. The fine structure the former has been well studied (10, 22) while 
that the latter has not yet been studied detail. 
addition these band systems, one further band progression was observed 
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extending from 1348 1234 This progression was found represent the 
progression the fourth positive group has been observed 
from tov’ 13. Most the isotopic bands could only measured 
the low temperature plates which the strong fourth positive bands 
are reduced width. 


DISCUSSION AND CONCLUSION 


The rotational and vibrational constants the and states 
obtained the foregoing analysis are summarized Table VI. both cases, 


TABLE 
VIBRATIONAL AND ROTATIONAL CONSTANTS THE AND ELECTRONIC STATES 
Constant 
1.266;+0.0010 
0.01872+0.00010 
0.000205 0.0001 
A+4y —1.1; +0.8, cm7! 
Were —0.0011; 


the and values are considerably smaller than the ground state. This 
accounts for the length the observed progressions. Further bands the 
and e—X systems would most probably found shorter wavelengths 
this region were studied. This was not possible the present work because 
the lack transmission the LiF windows wavelengths shorter than 
1200 

The transitions and heteronuclear molecules are 
many ways complementary each other: The branches possible one transi- 
tion according the electric dipole selection rules (AJ +1, 
cannot occur for the other. Both transitions are made possible 
interaction. The selection rule that states not combine with states 
does not apply intercombinations was shown Schlapp (23) and does, 
therefore, not introduce additional degree forbiddenness for 
Indeed the intensities the two band systems are the same order. 
transitions between states the same multiplicity are much weaker, like 
the bands (see (15)). bands have yet been 
found. Examples homonuclear molecules (or 
transitions electric dipole type have been observed for (see the paper 
III this series (17)) and (Vegard—Kaplan bands). They are very much 
weaker than the bands and, the case have not yet been 
observed absorption spite many attempts. 
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The ground state electron configuration can written 


where stands for the two shells. The two lowest unoccupied orbitals are 
and o*2p, both antibonding type. Since, shown the spectrum 
the orbital has energy not much below the lowest excited 
states would expected belong the following configurations 
(omitting which the same for all) 


[5] 


These configurations give rise the electronic states given the right. 
There question that the and states [2] are identified with 
the observed and states (upper states the fourth positive and 
Cameron systems). The upper states the two band systems here under 
discussion and must, all probability, correlated with configura- 
tion [3]. The state the same configuration presumably the 
upper state the Angstrém bands, (compare the energy level diagram Fig. 
197 (14)). The and states [4] may readily identified with the 
observed and states and the ‘II and states [5] with the 
observed and states. 

The electronic configurations are entirely analogous the low-lying 
configurations the molecule (see (12)). particular, the upper 
state the bands the analogue the state 
here considered. remarkable that spite the great similarity electron 
configuration, there such large difference intensity. Further theoretical 
work the intensities forbidden transitions indicated. 
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AURORAL HEIGHTS OVER WEST-CENTRAL CANADA! 


ABSTRACT 


complete analysis the parallactic photographs aurora taken 
Chesterfield, N.W.T., (63.3°N., 90.7°W.) 1933 show marked maximum 
the heights the lower limits both homogeneous arcs and bands and 
auroral structures with rays between 100 and 110 km., (b) mean height the 
lower limit arcs and bands 108.2 km., (c) mean height the lower limit 
auroral structures with rays 113.4 km., decrease height with 
increasing intensity, and (e) definite evidence change height with 
time after sunset. The results suggest that heights aurora observed from 
Northern Canada are slightly higher than aurora observed from Norway. 
This may due the locations the two regions relative the zone maxi- 
mum auroral occurrence. exceptionally low heights were observed 
Chesterfield. Remeasurements the 1933 photographs for Saskatoon (52.1°N., 
106.6°W.) indicate that one arc was definitely below the 70-km. level, and 
possibly below the 60-km. level. 


INTRODUCTION 


The numerous measurements height made Vegard and 
and Harang (8, 10, 11) aurora observed from Norway have 
shown characteristics that are assumed generally apply all regions. 
Data support this assumption are exceedingly scanty, only small numbers 
measurements having been made over other regions. During the Second 
International Polar-Year, 1932-33, about thousand pairs parallactic 
photographs aurora were taken Chesterfield, N.W.T., (63.3°N., 90.7°W.). 
Height measurements were made shortly after 200 pairs that were selected 
roughly random from the group. The results showed heights for the lower 
limit aurora that were substantially the same for Norway (2). the 
same time, Alty and Wilson (1) made small number height measurements 
Saskatoon (52.1°N., 106.6°W.). number these were much lower 
(50-70 km.) than any reported previous since that time. view the 
approaching International Geophysical Year, seems worth while 
report the complete analysis the Chesterfield photographs, and 
remeasurements the Saskatoon photographs. The Chesterfield measure- 
ments are also the only ones any number from place well inside the zone 
maximum auroral occurrence. 

The length the Chesterfield base-line was 32.14 km. The orientation 
(6.3°E. N.) was particularly good for height measurements arcs and 
bands because their horizontal extent from east west. They were taken 
during three periods major auroral activity—January 22-27, February 
26-28, and March mostly the evening hours. definite effort 
was made get photographs auroral structures satisfactory for height 
measurements. This meant that rapidly changing forms and ones with poor 
orientations for measurement relative the base-line were not photographed. 
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result, only about one out every two pairs had rejected because 
unsatisfactory orientation, too low angular altitude, motion during 
photographing, unsatisfactory exposure times, and on. 

The length the Saskatoon base-line was 112.5 km., orientated slightly 
north east. The long base-line resulted large parallactic displacements 
even for forms low angular altitudes. The orientation was unsatisfactory 
for structures with extent, the parallactic displacements running 
nearly parallel the edges arcs and bands. Photographs were taken 
January 30, February 25-26, and April and 20. Only were considered 
suitable for height measurements. 

The methods used measuring heights from parallactic photographs have 
been described detail auroral publications (6). They consist essentially 
projecting the auroral photographs tracing paper cm. the tracing 
corresponding one angular degree) and using the star backgrounds fit 
the tracings which are projections plane surfaces orthogonal 
systems small and great circles for the celestial sphere. means these, 
the altitude, the azimuth, and the parallactic displacement for selected auroral 
points are determined. Graphical methods are used also determine the 
height and the distance along the earth’s surface from the primary photo- 
graphic station the vertical projection the earth’s surface each auroral 
point. 

The heights several points (three four five depending the angular 
extent) along the lower edge each arc band, and the apparent lower 
and upper limits rays, are usually measured. The upper limits arcs 
and bands cannot determined with any accuracy because the lack 
brightness contrast between the upper edges and the sky background. Inter- 
mediate heights auroral structures are sometimes useful for determining 
the orientation space auroral structures, but serve other purpose 
far height characteristics are concerned. 

fairly critical attitude must taken toward the measurements each 
form. the values along the lower edge arc band vary irregularly 
and large amounts, ray structure not approximately parallel 
the lines the earth’s magnetic field, probable that the set photographs 
was not suitable for height measurements. Sometimes the parallactic angles 
are fault because low altitudes poor orientations relative the base- 
line the auroral structures. More frequently volume effect responsible, 
edges that apparently correspond the two photographs being actually 
different ones because folds and depth the auroral structures. set 
photographs yielding consistent values that are much smaller larger than 
normal more difficult eliminate. This situation most likely occur 
with rapidly changing auroral forms when failure take the two photographs 
simultaneously can result large errors the parallactic angles. 

Both the Chesterfield and Saskatoon measurements were examined with 
these sources error mind. All values for particular structures were dis- 
carded whenever there was definite evidence that one more these errors 
had appreciably affected the results. Some measurements were retained 
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that were suspect, since discarding them could mean the loss significant 
information. Cases with heights departing appreciably from the expected ones 
will discussed later. 


RESULTS 


The pertinent information 1378 heights the lower limits homo- 
geneous arcs and bands observed from Chesterfield are given Table 


TABLE 


FREQUENCIES OF MEASURED HEIGHTS OF THE LOWER LIMITS OF HOMOGENEOUS ARCS AND 
BANDS, AND OF AURORAL STRUCTURES WITH RAYS 


Homogeneous arcs and bands Structures 
km. Jan. Feb. Mar. Total total 
95-99 106 168 
100-104 126 208 
105-109 143 228 
110-114 108 171 


and are shown graphically Fig. preparing the table, heights were 
grouped ranges km., starting upward from km. The number occurring 
each height interval plotted against the median height for the interval. 
Included Table are the height frequencies the lower limits arcs and 
bands for each the three periods during which photographs were taken 
(January, February, and March), and the lower limits auroral structures 
with rays. Because the smaller number measurements ray structures 
they are not subdivided into monthly groupings. 

Heights the ranges and km. are definitely more numerous 
than for other ranges. The average height the lower limits the arcs and 
bands 108.2 the structures with rays, 113.4 km. Both the frequencies 
occurrence and the average values suggest heights for Canada the north 
the auroral zone that are slightly greater than for Norway (6, 10). 
Since less energetic particles can penetrate greater depths the atmosphere 
inside than outside the auroral zone this difference may actually exist. Other 
factors, such tidal effects the upper atmosphere, intensity the aurora, 
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Fic. The frequency distribution the heights the lower limits homogeneous arcs 
and bands. The number measurements 5-km. ranges, starting from km., are plotted 
against the median height for the range. 


and time from sunset, affect auroral heights such extent that not much 
weight can attached the slightly larger values from the Chesterfield 
measurements. 

The frequency distributions for some the Norwegian measurements 
show two maxima (6, 11). Egedal and Harang (4) have shown that the 
maximum for the larger height associated with flood-time, and the maximum 
for the smaller height with ebb-time the lunar tide. The Chesterfield photo- 
graphs for January were taken almost entirely during times when the tides 
were low; for February, when high; and for March, partly when high and 
partly when low. The average height the arcs and bands for January 
106.1 km.; for February, 112.0 km.; and for March, 110.0 km. These values 
support the belief that lunar tide present the atmosphere auroral 
levels. Since most the Chesterfield heights are for times low tide, the most 
frequent and the average values that have been deduced may smaller than 
normal for the region. should also noted that this supports the possibility 
larger auroral heights inside than outside the auroral zone. 

The Norwegian measurements indicate lower heights for the brighter 
more intense auroral forms (5). This also evident from the Chesterfield 
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measurements. visual intensity figure for aurora, ranging from for the 
weakest for the strongest, had been entered the photographic logbook. 
All forms with three more height determinations along their lower edges 
were divided into two groups according intensities—weak and moderate. 
Measured forms strong intensity were too few worth separating from 
the ones moderate intensity. Each these groups was again divided into 
height ranges km. form was considered belonging particular 
height range the majority the heights for fell within the range. Forms 
with heights scattered over more than two ranges were not included. This 
procedure eliminated from consideration both individual heights and forms 
where the measurements were probably error. The results these groupings 
are shown Table II. apparent from the table that forms moderate 
intensity are relatively more frequent the lower heights than ones low 
intensity. 
TABLE 


NUMBERS OF HOMOGENEOUS ARCS AND BANDS, AND OF AURORAL STRUCTURES WITH RAYS OF 

WEAK AND OF MODERATE INTENSITY HAVING THEIR LOWER LIMITS LARGELY WITHIN THE 

INDICATED HEIGHT RANGES. FORMS OF MODERATE INTENSITY ARE RELATIVELY MORE FREQUENT 
AT LOWER HEIGHTS THAN ONES OF WEAK INTENSITY 


Homogeneous arcs and bands Structures with rays 
Heights, 
km. Weak Moderate Weak Moderate 


progressive decrease height with time after sunset, shown Egedal 
and (page Ref. 6), could not deduced from the Chesterfield 
measurements. Figure shows plot the average heights arcs and bands 


HEIGHT 


HOURS AFTER SUNSET 

Fic. mean heights the ‘lower limits homogeneous arcs and bands for hourly 
periods centered the indicated number hours after sunset. The quantities brackets are 
the number heights used evaluating each mean. 
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for each hourly interval centered the indicated hours after sunset. The 
variations are possibly due intensity changes. mentioned previously, 
the photographs were taken during periods above normal auroral activity. 
Each night the aurora went through several periods maximum and minimum 
activity, each lasting for hour more but occurring irregularly from night 
night. Subdivisions the data for Fig. according intensity were equally 
inconclusive, and are not included. 

attempt was also made detect change height with geomagnetic 
latitude, has been demonstrated (9, 10) for latitudes the south 
the auroral zone. Heights arcs and bands the south Chesterfield 
were compared with those the north Chesterfield. significant differences 
could detected. Again, intensity variations may obscuring factor. 
Displays Chesterfield developed much the same way they the 
south the auroral zone—building the direction toward the auroral 
zone and spreading overhead and into the part the sky away from the zone 
the peak their activity. Obviously, any increase height due decreas- 
ing geomagnetic latitude could balanced the decrease height due 
the higher intensities such times. 

There was one case unusually low aurora. This occurred January 
for arc-like structure intensity and with pink and red shades. 
The heights for five points separated about 5-degree intervals azimuth 
were 77, 76, 71, 72, and km. The parallactic angles varied from 4.5 5.2 
sufficiently large give reasonably reliable height measure- 
ment. The only criterion that might used for discarding these measurements 
the possible lack simultaneity taking the two photographs, since such 
forms develop and change rapidly. Similar forms have been measured 
heights below km. Harang and Bauer (7). The values for this arc are not 
included Table and Fig. 

Unexpectedly large heights were not numerous. Most those associated 
with arcs and bands were probably due volume effects which could not 
detected with certainty inspection the photographs. addition there 
were cases forms, described variously diffuse bands, corona-like sur- 
faces, and surfaces, with heights between about 150 and 300 km. All these 
had large parallactic angles and all occurred after periods much auroral 
activity. Lack clearly defined edges, movement and volume effect undoubt- 
edly reduce the reliability the measurements, but the fact that similar types 
structures were not observed lower heights indicates that the measure- 
ments are not sufficiently error discarded. Presumably some the 
diffuse structures that form active display dies away are much greater 
heights than the sharply defined forms preceding and during the display. 

The upper limits rays associated with arcs, bands, and draperies varied 
from 130 over 300 km. There were cases rays exceeding 400 km., but 
these values are not considered particularly reliable. None were sunlit. Photo- 
graphs rays completely isolated from other structures were not obtained. 

The Saskatoon values were checked having assistant remeasure all 
the plates. Where his values the heights departed appreciably from the 
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original values, the writer measured the heights once again. The agreement 
with the original values was reasonably good for all cases where twists and rays 
the auroral structures had been used determine the parallactic angles. 
For structures with regular lower edges, nearly parallel the plane the 
parallactic displacements, there was notable lack agreement. the same 
criteria for the Chesterfield photographs had been used, these structures 
would have been discarded. Three sets successive photographs the same 
arc with heights reported Alty and Wilson (1) between and km. 
are not exclusively this group. The exposures were long for this type work 
(between and sec.), but the lower edges the photographs are well 
defined and show evidence motion during the exposures. The directions 
the edges are not particularly favorable the determination reliable 
parallactic angles, but their values are internally consistent for the three 
sets photographs. Repeated measurements indicate that this arc was 
definitely below km. and possibly below km. Unless some unsuspected 
reason can found for discarding these measurements, must accepted 
that auroras occasionally occur the south the auroral zone much lower 
levels than usually assumed. 
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THE THE BOTTOM THE FRAUNHOFER 
LINES! 


ABSTRACT 


The line profiles were measured scanning photomultiplier and exit slit 
the second order 21-ft. grating. The large distortions resulting from 
finite resolution were evaluated means absorption cell containing 
sodium vapor with helium pressure 170 mm. This produced artificial 
line whose central intensity was known zero. matching its profile 
carefully that the solar line the central intensities and were found 
(5.0+0.6) and (4.8+0.6) terms the nearby continuum. 


INTRODUCTION 


One the most valuable results coming from studies the twilight sodium 
emission estimate the number free sodium atoms per square centi- 
meter column the upper atmosphere (1, 5). The data required are the 
scattered intensity, the scattering power sodium atom, and the intensity 
the incident sunlight. The first two are known with reasonable accuracy, but 
the last very uncertain. The intensity the continuous spectrum well 
known (6), but the radiation active exciting the twilight scattering the 
bottom the Fraunhofer lines. The ratio the central intensity the con- 
tinuous intensity (which will call not known with any accuracy. 
The best and most recent measurement that Priester (8) who used 
Lummer plate crossed with high dispersion spectrograph observe the line 
profiles. The recording was done photographically, and separate measurements 
were made for various parts the sun’s disk. Substantial corrections had 
made for finite resolution and scattered light. 

was felt that independent measurement different method, using 
light from the whole sun, would desirable. This consisted measurement 
the line profiles using 21-ft. concave grating, and second measurement 
under similar conditions but with the central intensity reduced zero 
sodium vapor absorption tube. The difference then gave the central intensity 
directly. Direct photoelectric photometry was used order eliminate 
uncertainties the standardization photographic plates. The agreement 
our results with Priester’s highly satisfactory. 


PRELIMINARY MEASUREMENTS 


was thought first that direct photometry with the grating would yield 
the required accuracy. The spectrograph Eagle type using 15,000 
lines per inch grating in. wide. photomultiplier and exit slit were mounted 
that they could scanned for short distance tangent the Rowland 
circle the second order. coelostat was arranged bring the sunlight 
without forming image, that the results were for the whole sun. 
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simple double monochromator was set ahead the grating instrument 
and adjusted pass about band. Errors caused scattered light were 
thus greatly reduced. The detector was 1P21 photomultiplier followed 
simple but stable and linear direct-current amplifier. similar system 
fixed position the first order served monitor for changes the incident 
intensity. The circuit the amplifier was similar one already published 
(4), but with the electrometer stage omitted since the photomultiplier load 
resistance was only megohms. two-pen recording potentiometer allowed 
both signals recorded one chart. 

The theoretical resolving power the grating the second order 180,000; 
practice were able achieve 100,000. The resulting distortion the 
line profiles should however the values were about 20% rather 
than about 6%. Because the double monochromator this could not 
attributed scattered light. The instrumental line profile was therefore 
measured using the 5790 line mercury 198. (The loan the source was 
kindly arranged Dr. Douglas the National Research Council and 
Dr. Vallance Jones this Department.) The width halfway down was about 
0.06 giving the resolving power already quoted, but there were faint wings 
extending nearly 0.4 out from the center. These would have very little 
effect the sharpness lines emission spectroscopy, but would result 
serious filling-in the absorption lines, was observed. 

attempt was made correct the observed profiles method resembling 
that Hardy and Young (3). For and the resulting values were 
8.5% and 7.7%. These cannot considered all accurate, since the correction 
about twice the final result. 


THE SODIUM VAPOR TUBE 


The correction the center line was evaluated directly forming 
absorption line whose shape matched that the solar line but whose 
central intensity was known zero. Any observed intensity was spurious 
and was course equal the required correction. The artificial line was 
produced passing sunlight through tube containing sodium vapor and 
helium the correct pressures. The helium was required make the line 
broad enough and smooth enough, since pure sodium vapor gives very 
narrow line the pressure range that was used. 

Since glass attacked hot sodium vapor, the tube was lined with 
glaze used for vapor lamps the General Electric Company (generously 
given Mr. Schmies). This normally used with special glass 
matching thermal expansion, but was found satisfactory with Pyrex. The 
glaze must fired temperature the softening range the glass, and 
difficult keep the flat glass for the windows from distorting. This problem 
was solved making layer about mm. thick powdered graphite 
iron plate and packing pressure with the window blank. Heat could then 
applied from above and below. The uncoated surface remained adequately 
smooth and kept its polish. The coated surface was finely pebbled and had 
polished with rouge felt pad. The two windows were found remove 
about half the light from collimated beam. 
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The diameter the tube was mm. and the length mm. Many the 
first attempts failed when the tubes cracked after the windows were sealed 
on, because strains introduced when the glaze was fired. This difficulty was 
removed annealing the parts oven before assembly. The completed 
tube was also annealed. Because considerable pressure helium was 
used, the tube was not carefully outgassed. was evacuated single- 
stage oil diffusion pump, flushed with helium, evacuated again, and filled 
the final pressure 170 mm. Sodium purified vacuum melting and distilla- 
tion was introduced from thin-walled ampoule. 

The profile the absorption line different temperatures was calculated 
from equation (97) (book Mitchell and Zemansky, Reference 7). 
This equation assumes Lorentz broadening, but should give good approxi- 
mation the pressure used. For the pressure helium used and for tempera- 
ture near 180°C., equations (104) and (96) give Lorentz damping 
ratio 1.5. Some lines were also calculated for 0.5 and 1.0 find the 
effect the line shape. very good approximation temperature 
giving the same shape could found for each value. Since there optimum 
value the largest was used, that the temperature could kept rela- 
tively low. Calculations were made with the aid Appendix MZ. 
The vapor pressures used were those Edmondson and Egerton (2). 


MEASUREMENTS AND INTERPRETATION 


The observing procedure was very simple once the spectrograph had been 
adjusted give the sharpest possible lines. The two lines were observed 
different runs, since too much time would otherwise have been lost 
scanning between them. The line was scanned several times the sodium 
tube slowly cooled. Another measurement was made with the tube removed 
from the optical path. The position the continuum was estimated 
drawing line through the trace about from the center. The observed 
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Fic. The observed central intensities for and function the temperature 
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residual intensity was the ratio the minimum deflection the deflection 
for the continuum. The results are plotted Fig. function tempera- 
ture the tube. Since there break the curves indicating the desired 
point, necessary estimate this point with the aid the calculated 
absorption line profiles. 

The line desired has the same shape the solar line but displaced down- 
ward that its central intensity essentially zero. The shape used was that 
found Priester (8) for sin 0.60; the angle between the line sight 
and the normal the sun’s surface. average shape for the whole disk 
would preferable, but did not seem worth the considerable labor calculat- 
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Fic. the line 5896 The upper solid curve with points Priester’s 
measurement. The lower the same displaced down. The three dotted curves are calculated 
for sunlight passing through the sodium vapor tube 160, 170, 180°C. (left right). 

Fic. Profiles the line 5890 The upper solid curve with points Priester’s 
measurement. The lower the same displaced down. The three dotted curves are calculated 
for sunlight passing through the sodium vapor tube 150, 160, 180°C. (left right). 
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ing it. However, the average central intensity was calculated, allowing for 
the contribution each region the total area and for limb darkening. The 
latter was represented the approximation suggested Minnaert (6), witha 
coefficient 1.6. For and the results were 6.6% and 5.8%. 
would expected, these are close the values observed sin 0.60. 

The line shapes actually achievable are shown the dotted curves Figs. 
and Also shown are the shape the solar line found Priester and 
the same curve displaced down until touches the horizontal axis. The dotted 
curve most nearly matching the desired curve the region the instrumental 
profile the one desired. The best temperatures are 167°C. for and 158° for 
with uncertainty about +2° both cases. Reference Fig. then 
double the error figures allow for uncertainties temperature measurement, 
the vapor pressure relation for sodium, and the curves Fig. Account 
must also taken the fact* that the intensity from the line center 
only about 92% the continuum (9). The final result therefore 


DISCUSSION 


There little doubt that should less than satisfactory that our 
result agreement with this. The results Priester already quoted 
6.6% and 5.8%) differ somewhat more and are little higher. 
Our values may taken giving general confirmation Priester’s, but 
indicating that his may slightly too high. This not impossible, since the 
corrections for instrumental effects, though far smaller than ours, were still 
appreciable. 

possible that slightly better accuracy might achieved averaging 
several measurements this type, but far greater improvement could 
made starting with spectrograph higher resolution. Another possibility 
use double-beam system with one beam passing through the sodium 
cell and with rapid alternation the two beams the phototube. measure- 
ment this method now progress. 
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ENERGY AND ANGULAR DISTRIBUTIONS PHOTONEUTRONS 
PRODUCED 70-MEV. 


ABSTRACT 


The angular distributions the photoneutrons produced 70-Mev. X-rays 
seven elements have been measured with zinc sulphide lucite scintillation 
detector. For the heavy elements the distributions are essentially isotropic, 
indicating the predominance the evaporation process, while for the light 
elements there also anisotropic component peaked 90°. The energy 
distributions the photoneutrons produced copper and lead targets have been 
determined with nuclear emulsions. shown that most the neutrons can 
attributed evaporation process which governed constant nuclear 
temperature, the temperature found for copper being 1.2 Mev., and for lead being 
1.0 Mev. high-energy tail the energy distributions attributed direct 
interactions. 


INTRODUCTION 


Studies photonuclear reactions can conveniently divided into studies 
the reaction cross section function photon energy, and studies the 
energy and angular distributions the emitted particles. The latter type 
experiment first all expected yield information the relative impor- 
tance compound-nucleus formation and direct ejection. The study the 
compound-nucleus processes particular value that the excitation the 
compound nucleus may more with photons than with incident 
particles. The direct ejections are course intimately connected with the 
photon-absorption process, and their study expected give more detailed 
information the absorption process than can obtained from total cross- 
section measurements. 

The energy and angular distributions particles emitted compound 
nucleus are characteristically different from those ejected direct interaction. 
The energy distribution particles evaporated from compound nucleus 
expected somewhat maxwellian shape, while directly-ejected particles 
are expected follow approximately inverse power law (5, 13). Direct 
particles thus form high-energy tails the energy distributions. The angular 
distribution for evaporated particles expected isotropic, unless there 
are only few states the residual nucleus available, which case however 
the angular distribution still symmetrical about 90° (10, 25). For direct 
ejection (electric dipole absorption) the angular distribution expected 
the form where the angle which the neutron emitted 
with respect the X-ray beam and and are constants. For protons, 
interference can give rise peak forward 90°. 

Although the energy and angular distributions are characteristic the 
interaction process, should noted that unambiguous separation the 
evaporation and direct processes often experimentally impossible. Indeed 

1Manuscript received August 
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may futile seek sharp delineation, since processes intermediate 
between conventional compound-nucleus formation and direct interaction 
undoubtedly occur. 

There have been several studies the energy and angular distributions 
photoprotons, but relatively few photoneutrons, probably owing the 
difficulty measuring neutrons the unfavorable conditions which exist 
experimentally. Section angular-distribution measurements made here 
with zinc sulphide lucite detector are described and Section III energy 
distributions for copper and lead are presented. 


II. ANGULAR DISTRIBUTIONS 


Experimental Arrangement 

measuring photoneutrons one must contend with very large background 
intensities secondary electrons and gamma-rays which originate the target 
placed the X-ray beam produce the photoneutrons. These intensities 
exceed the neutron intensity factor roughly 1000. However possible 
obtain the required discrimination against gamma-rays from the fact that 
the specific ionization recoil protons much higher than that secondary 
electrons, provided that counting rates are kept low enough avoid 
pile-up pulses during the high-intensity bursts X-rays from the synchro- 
tron. suitable detector has been described Hornyak (11). Grains zinc 
sulphide diameters the range microns are molded lucite, which 
serves produce recoil protons and transmit the emitted light photo- 
multiplier. The sensitivity this detector found increase almost linearly 
with neutron energy (24). This increase seems due the fact that the 
average range the recoil protons increases with neutron energy greater 
rate than the cross section for neutron—proton scattering decreases. 

The experimental apparatus used here consisted target placed the 
X-ray beam, zinc sulphide lucite detector distance cm. from the 
target, EMI 6262 photomultiplier, turntable rotate the detector about 
the target, and associated electronic apparatus. large shield lead and 
paraffin was used shield against the field stray X-rays and neutrons from 
the donut the synchrotron. lead collimator with half-angle about 
was used define the X-ray beam. Lead shielding thickness approximately 
in. was used about the photomultiplier. 

order reduce the pile-up small pulses, lead absorbers were also used 
between the target and detector. Aluminum absorbers alone, intended 
reduce the electron intensity incident the detector, were not sufficient. 
Some integral bias curves illustrating the effect absorbers are shown 
Fig. Since the electron and gamma-ray intensities are strongly forward 
their angular distributions, the thickness absorber required must deter- 
mined experimentally for the smallest angle for which proposed take 
measurements. For the present work this angle was 45°; from Fig. was 
decided that in. lead absorber was ample this angle. 

further measure against pile-up was taken maintaining the length 
the X-ray pulse large value (about 200 and thereby reducing the 
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Counts 


Discriminator bias 


Fic. Integral bias curves with the zinc sulphide lucite detector angle 45° the 
X-ray beam. The source photoneutrons was lead target. The curve for absorber 
1/4in. aluminum shows the effect pile-up small pulses. 


instantaneous intensities. This was accomplished cutting off the synchrotron 
r.f. when the magnetic field was near its peak value. circuit 
was also tried, but did not prove necessary. 

Large targets were required produce reasonable neutron counting rates. 
The targets used are listed Table Two different beryllium targets were 


TABLE 
TARGETS USED ANGULAR DISTRIBUTION MEASUREMENTS 


Radius Length Weight 
Element (cm.) (cm.) (gm.) 
Lead 1.21 7.9 328 
Tin 1.30 9.1 306 
Copper 1.44 408 
Iron 1.26 306 
Aluminum 1.25 8.0 108 
Carbon 2.37 8.4 246 
Beryllium (1) 10.2 490 

Beryllium (2) 5.08 91.3 


used. the first, which had dimensions in. in. in., there was con- 
siderable smearing the emergent neutrons due scattering the neutrons 
within the target. second target, the form hollow cylinder length 
in., inside diameter 2.4 cm., and outside diameter 4.3 cm., was then obtained 
order reduce this effect. 


i 10% ° 
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Results 

The angular distributions the photoneutrons from the lead, tin, copper, 
iron, aluminum, and carbon targets are shown Fig. and the distributions 
for the two beryllium targets are shown Fig. The counting rates shown 
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Fic. The angular distributions measured with the zinc sulphide 
lucite detector. 


Beryllium 


min. 


Angle 
Fic. The angular distribution photoneutrons from beryllium targets measured 


with the zinc sulphide lucite detector. The upper curve for the rectangular target, and the 
lower curve for the hollow cylinder. 
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for the various elements have not all been normalized one another. Angular 
distributions for most the targets were obtained several different days, 
and these results were then averaged. single day’s run consisted two 
three complete surveys the angular positions, this being preferred one 
survey with longer counting periods order that any systematic drifts the 
counting rate would apparent. The X-ray beam intensity was monitored 
with ionization chamber, and the output recorded moving chart. 
The counting rates with the X-ray beam on, but with target, have been 
subtracted from the total counting rates give the net rates shown Figs. 
and The no-target background rate amounted about half the total 
counting rate for the carbon target, and was relatively less important for the 
other targets. The errors assigned the experimental points Figs. and 
are root mean square deviations which include the effects statistical un- 
certainties both the total and background counting rates. 

The solid curves shown Figs. and have been obtained from least- 
squares fits the experimental data curves the form The 
accuracy the data does not warrant attempts fit them more compli- 
cated curves, although forward shift the peak seems present most 
the angular distributions. (This shift may due part pile-up gamma- 
ray pulses.) There theoretical well practical reason for fitting the 
data simple curve. This based the theorem that the absorbed 
photon has angular momentum the outgoing intensity distribution may 
contain spherical harmonics only, provided that both the incident 
beam and the target are unpolarized. (See for example Sachs (22).) Since only 
photons small angular momentum are absorbed appreciably, simple 
angular distribution the outgoing particles expected. 

quantitative measure the anisotropy provided the ratio 
although the accuracy this quantity low, particularly when small. 
The values obtained from the present results are summarized Table II. 


TABLE 
EXPERIMENTAL VALUES FOR b/a 
Correction factor Corrected 

Target for self-scattering b/a 
Lead 1.10 —0.08+0.08 
Tin 1.08 
Copper 1.48 0.23+0.15 
Iron 1.35 0.09+0.25 
Aluminum 0.36+0.29 
1.8 1.6 +0.8 

eryllium 2.6 
Beryllium (2) 1.35 1.2 +0.4 


Corrections for the self-scattering the photoneutrons within the targets 
have been made assuming that one can replace the actual problem with 
point source neutrons located the center spherical target and emitting 
according The corrections obtained this manner are also 
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given Table II. can shown that the particular geometry employed 
here the loss angular resolution due the finite size the target, the 
finite size the detector, and scattering the lead absorbers small compared 
the effect self-scattering. 

the heavy elements the angular distributions are seen nearly 
tropic, indicating the predominance the evaporation process. These results 
are general agreement with the results Price and Kerst (20) and 
Palmer (17), both obtained with detectors involving the moderation the 
fast neutrons and the subsequent capture thermal neutrons. Geller, Halpern, 
and Yergin (7) have obtained similar results with zinc sulphide lucite 
detector. the other hand, Price (19) and Poss (18) with Al(m, 
and Johannson (12) with zinc sulphide lucite detector have 
observed b/a some heavy elements. order explain these various 
results, only necessary suppose that the anisotropic neutrons the 
average have much higher energy and are much fewer numbers than the 
isotropic neutrons. Then the detector sensitive all low-energy 
neutrons, the observed ratio essentially zero, while the detector 
sensitive only high-energy neutrons, the observed ratio characteristic 
the anisotropic process. 

The variation the observed b/a with mass number can also under- 
stood this basis. With detector Price (19) found increase 
with increasing Since the evaporated neutrons tend have lower energies 
heavy nuclei, owing the lower nuclear temperatures, this trend probably 
caused the rapidly decreasing sensitivity his detector evaporated 
neutrons, rather than change the value for the direct interactions 
alone. the present results can said show trend with all, 
that b/a decreases with increasing Since the response the detector used 
here not critically dependent the energy the evaporated neutrons, the 
present results can interpreted giving qualitative indication the 
relative importance the evaporation and direct processes over wide range 
nuclei. 

large anisotropic component carbon has been found the present 
investigation, contrast with previous results obtained with the moderator 
type detector (17, 23). Both this large anisotropic component and the long 
tail observed the high-energy side the cross section curve (1) 
probably can attributed rather large direct process. The remarks 
Mann, Stephens, and Wilkinson (15) the angular distribution the photo- 
protons produced 24-Mev. X-rays should also noted this connection. 

large peak has been found for beryllium the present work and 
Johannson (12). Since the value found with sulphide lucite 
greater than that found with the moderator type detector 
appears that the anisotropic neutrons are the average more energetic 
than the isotropic neutrons. This result suggests that the effect the low- 
energy anisotropic component, known present upon the absorption 
2.76-Mev. photons (9), has been masked the usual behavior high photon 
energies. The angular distributions neutrons resulting from the absorption 
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photons intermediate energies are not known, although the disintegration 
processes found Carver, Kondaiah, and McDaniel (4) make seem unlikely 
that there would marked anisotropy, least for photons. 
Photoneutron Yields 

The experimental data the angular distributions were also used give 
information photoneutron yields function atomic number. making 
this analysis was necessary normalize the various counting rates. Correc- 
tions were applied for the absorption the X-ray beam the targets and for 
the lateral variation the beam across the target. attempt was made 
correct for the effect inelastic scattering the photoneutrons within the 
target; inelastically scattered neutrons were detected with reduced efficiency. 
was found that the yield per mole varied approximately with 
This result accord with other measurements: with the moderator 
type detector, the yield per mole found vary approximately 
(20, 17, 23), while only high-energy neutrons are counted, the yield found 
vary —Z) (12). 

ENERGY DISTRIBUTIONS 


Experimental Procedure and Results 

many methods for the determination the energy distribution 
neutrons necessary measure the total ionization loss recoil protons. 
For photoneutrons also necessary provide sharp discrimination against 
gamma-rays. These requirements often place lower limit the neutron 
energy which can measured. For example proton telescope which both 
and dE/dx are measured cannot successfully used below about Mev. 
Roughly the same energy limitation applies two-crystal delayed-coinci- 
dence spectrometer which the discrimination against gamma-rays achieved 
arranging that the second pulse due scattered neutron occurs time 
appreciably later than the second pulse due scattered gamma-ray. Nuclear 
emulsions seem the author the best means presently available for 
measuring the complete spectrum the photoneutrons produced high- 
energy X-ray beams. 

Ilford plates dimensions in. in. and emulsion thickness 
were exposed neutrons produced copper and lead targets irradiated 
the X-ray beam. The plates were exposed angles 60°, 90°, and 120° 
and distance cm. from the target the center the plate. Back- 
ground exposures were made with the target removed. lead shield thick- 
ness 1/8 in. was interposed between the target and plates order reduce 
the electron intensity incident the plates, correction the observed 
energy distributions being required for this thickness. The large shield lead 
and paraffin and the collimator mentioned above were used. 

The targets were chosen small enough that correction the observed 
energy distribution would necessary because inelastic scattering them. 
The copper target was cylinder length in. and radius 3/8 in. The lead 
target was roughly spherical shape, with the radius corresponding about 
0.05 mean free path for inelastic scattering. 
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The following criteria were used for the selection recoil-proton tracks. 
Each track had lie square pyramid half-angle 20° the unprocessed 
emulsion, the axis the pyramid being parallel the edge the plate, and 
the apex being located the beginning the track. The track had start 
and stop the emulsion. Individual calculations were made each track, 
both because the rather large acceptance angles, and order make sure 
that tracks were lost. The relation Rotblat (21) 
was used ascertain the proton energies, and the neutron energies were then 
calculated using the known angle the track. All tracks due neutrons 
energies less than 0.5 Mev. were ultimately discarded, because the non- 
agreement different observers this energy range. The experimental 
distributions were corrected for the escape proton tracks through the sides 
the emulsion, and for the scattering cross section. The 
areas scanned and the total numbers acceptable tracks are given Table 


TABLE III 
TRACKS USED ENERGY-DISTRIBUTION MEASUREMENTS 
Area scanned Number 
Target 
Copper 217 991 
Lead 140 259 
None 166 


Copper 


Neutron energy (MEV) 


Fic. energy distribution photoneutrons from copper. The histogram represents 
the present results, and the circles the results Byerly and Stephens. The solid and dashed 
curves are calculated for the evaporation first and second neutrons, respectively, their sum 
being normalized the histogram. 
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Number neutrons 


@ 


Neutron energy (MEV) 


Fic. The energy distribution photoneutrons from lead. The solid and dashed curves 
are calculated for the evaporation first and second neutrons, respectively. 


The experimental neutron energy distribution for copper* shown Fig. 
and the distribution for lead shown Fig. These distributions have been 
corrected for the no-target background, which was found appreciable 
only the energy interval Mev. The errors shown are root mean square 
deviations based the number tracks observed. The data shown for copper 
are averages results obtained angles 90° and 60°, there being 
significant differences between the two distributions the energy range 
0.5-5 Mev. For higher energies there seemed more tracks the 90° 
plate. The data shown for lead were obtained 120°. 


The Evaporation Process 
From formulae given Blatt and Weisskopf (2), the energy distribution 
for evaporated photoneutrons expected the form 


where the neutron energy, the neutron mass, the cross section 
for formation the compound nucleus absorption neutron energy 
the photoneutron threshold, the photon energy, Emax the maximum 
photon energy, the spectral distribution X-rays, the sum the 
partial cross sections for which the first particle evaporated neutron, 
w(U) the level density the nucleus excitation and finally 


*The author wishes acknowledge the generous assistance Dr. Livesey measuring 
the tracks used for this distribution. 
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assumed that the nuclear temperature independent the excitation 
energy variety experimental results seem suggest (14), then the 
level density given 


[3] w(U) w(0) exp(U/T). 
Using equation [3] [1] and [2], one obtains 
where 


Here the cross section for the inverse reaction has been assumed indepen- 
dent neutron energy. This independence cannot expected hold strictly, 
especially view the decreasing mean free path nuclear matter for 
nucleons increasing energy which envisaged Feshbach, Porter, and 
Weisskopf (6) and calculated for example Morrison, Muirhead, and Rosser 
(16). must also remembered that this cross section refers the capture 
neutrons the residual nucleus when excited state. 

Equations [4] and [5] can used determine the nuclear temperature 
from the experimental data process successive approximations. The 
first approximation consider the incident photons monoenergetic, 
that varies linearly with From plot vs. one obtains 
first approximation the temperature, which can then used evaluate 
and hence obtain the next approximation The integral not 
sensitive that the process rapidly convergent. 

The procedure can modified include the second neutrons emitted 
(y, 2n) reactions, assumed that the residual nucleus after the emission 
two neutrons has the same temperature the nucleus left after the emission 
the first neutron. One obtains double integral similar [5], which denote 
fe; for the determination the nuclear temperature one now replaces 
higher order reactions probably can neglected. 

The nuclear temperature found this way for copper 1.2+0.1 Mev., 
and for lead Mev. These values the temperatures seem 
significantly higher than those found Graves and Rasen (8) for the inelastic 
scattering 14-Mev. neutrons, namely Mev. for copper and 
0.76+0.08 for lead. 

The neutron energy distributions, calculated from equations [4] and [5] 
with the above values the nuclear temperatures, are shown the solid 
curves Figs. and The dashed curves represent the results the analo- 
gous calculations for the second neutrons. The sums the two curves have 
been normalized the experimental distributions. The results Byerly and 
Stephens (3) for copper are shown circles Fig. The agreement with the 
present results very good, except possibly high energies; this can 
attributed their relatively low maximum X-ray energy (24 Mev.). The only 
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other photoneutron energy distributions which have been reported are those 
Price (19) for silver, based 129 tracks, and for bismuth, based 196 
tracks. 


The Direct Effect 

evident from Figs. and that there are high-energy tails the 
energy distributions which cannot accounted for conventional evapora- 
tion process. These tails can attributed direct interaction process, 
possibly the Courant type. The distribution for copper also appears have 
broad peak about 6-7 Mev. similar peak can seen the results 
Byerly and Stephens slightly lower energy. For lead similar peak 
5-6 Mev. observed, but since there are only seven tracks represented this 
peak, may not statistically significant. The results for bismuth, reported 
Price, show similar peak 6-7 Mev., while his results for silver are 
inconclusive. 

tempting assume that this broad peak has the form 


that the energy the peak corresponds the energy the giant resonance 
relative the photoneutron threshold. When weighted with the 
peak given equation [6] shifted about one Mev. towards the low- 
energy end. Equation [6] based the assumption that the shape the 
cross-section curve for the photoemission these neutrons approximately 
the same that observed for the total (y, reaction; this consistent with 
the view that even when compound nucleus formed, the primary absorp- 
tion process involves one nucleon only. Further detailed speculation does not 
seem warranted until further experimental evidence forthcoming. 


IV. CONCLUSIONS 


From the angular distributions possible conclude that heavy 
nuclei, the evaporation photoneutrons from compound nucleus much 
more important than direct ejection. The results for carbon and beryllium 
probably indicate large direct effects these nuclei. The relative neutron 
yields agree with previous results. 

The energy distributions copper and lead can largely accounted for 
evaporation process governed constant nuclear temperature rather 
than temperature which increases with the excitation energy the resid- 
ual nucleus. The temperature found for copper 1.2 Mev. and for lead 
1.0 Mev., both slightly higher than the temperatures found from the inelastic 
scattering 14-Mev. neutrons. The high-energy tails the energy distribu- 
tions are attributed direct effect, possibly the Courant type, while 
broad peaks about Mev. may due toa special direct process. 
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THE ZERO-POINT KINETIC ENERGY LIQUID 


ABSTRACT 


The zero-point kinetic energy liquid helium has been calculated from the 
interatomic potential, the latent heat vaporization, and atomic distributions 
derived from neutron diffraction measurements. Calculations were carried out 
for two liquid temperatures and several published interatomic potential func- 
tions. resulting values the lie between 9.0°K. 


and 12.6° 
experimental value the zero-point kinetic energy liquid helium can 
obtained? from the radial distribution function, the interatomic potential, 
and the latent heat vaporization. the energy system helium atoms 
taken zero when they are rest and well separated, and the departure 
from ideal gas neglected for the moment, the internal energy per atom 
gaseous helium temperature and the enthalpy per atom 


where the pressure and the volume per atom. The non-ideal character 
the gas was taken into account using 
dB/dT)/U 


where the second virial coefficient per mole, the molar volume, and 


stands for the correction term. 

The kinetic energy per atom liquid helium temperature can 
temperature small compared the Debye characteristic temperature. 
Under these conditions nearly equal the (6). 


The potential energy per atom 


where the radial distance from atom, p(r) the atomic density, and 
V(r) the interatomic potential helium. 
The enthalpy per atom the liquid equal 


where the volume per atom the liquid. The latent heat vaporization 


per atom 
[3] 


received July 29, 
from Division Energy Canada Limited, Chalk River, 
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Equations [1], [2], and [3] yield 


convenient write equation [4] the form 
with 


independent the interatomic potential chosen for the calculation; 
which represents the potential energy the same way represents the 
kinetic energy, depends directly the potential function and the atomic 
distribution. Calculations have been carried out with the potential functions 
listed Table where the notation follows Reference (3) and using 


TABLE 
List OF POTENTIAL FUNCTIONS 
Desig- 


determined experimentally (2) 4.24°K. and nominally 2°K. (average 
results 1.65, 1.95, and 2.25°K.). The values are given Table II. 


TABLE 
VALUES 
Helium 
degrees degrees degrees degrees degrees 
3.33 0.00 —4.24 
4.24 7.07 0.26 —2.05 


Latent heat contributions were deduced from 23.2 joules/gm. 2°K. and 21.8 
joules/gm. 4.24°K. (5). The values required for were taken from 
Reference (4). 

The resulting values 2°K. and 4.24°K. are given Table III. Because 
the potentials increase very rapidly decreases below the finite root 
V(r) small errors p(r) this region can contribute appreciably the 
integral. Inspection the experimental rp(r) curves suggested that the non- 
zero values for 2.2A are spurious. Accordingly most the integrals were 
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TABLE III 
VALUES FOR VARIOUS POTENTIALS 


degrees degrees degrees degrees 
MR-1 15.9 11.7 14.4 12.4 
MR-5 13.1 8.9 12.0 10.0 
S-K 13.4 9.2 11.8 9.8 
LJ-1 14.3 10.1 12.4 10.4 
LJ-1* 13.3 9.1 12.1 10.1 


calculated with p(r) for 2.2A. check this assumption, Morse 
potentials were fitted MR-5 and LJ-1, and the wave mechanical density 
diatomic combination was calculated for each these potentials. Compari- 
son the calculated and observed p(r) justified the cutoff 2.2A for 
but showed that perhaps the observed p(r) should used 2.2A 
for LJ-1. The results retaining the observed p(r) atr 2.2A are shown 
starred Table III. These are preferred over the unstarred values for LJ-1. 

More figures are given than the accuracy warrants order show the 
variation caused the different potentials. Errors are difficult assess but 
may well degrees. There some doubt about the latent heat the 
higher temperature. instead Kistemaker’s result use lower value (1) 
20.3 joules/gm. the effect increase for 4.24°K. 0.5°. 

The quantity represents the kinetic energy the liquid the tempera- 
ture the measurement. Because includes some temperature excitation 
slightly larger than the zero-point temperature and expected increase 
with temperature. The 4.24°K. values Table III exceed the 2°K. values 
average 0.7° which more than would predicted for simple structure. 
For example, liquid helium behaved like simple Debye solid with charac- 
teristic temperature 26°K. value commonly assumed for helium) about 
one half the heat content constant volume would contribute the tempera- 
ture excitation. Expressed degrees manner comparable with the heat 
contents for the Debye model are 0.035° 2°K. and 0.626° 4.24°K., giving 
predicted change (0.626—0.035)/2 0.3° between these two tem- 
peratures. Some all the difference between 0.7° and 0.3° may arise from 
experimental error. the other hand the energy associated with the lambda 
transition may contribute the difference. 

upper limit 0.5° for the amount which 2°K. larger than the 
zero-point temperature obtained from the integral specific heat under 
saturated vapor pressure from 2°K. This possibly gross overestimate 
and view the experimental errors the 2°K. values Table III may 
considered represent also the zero-point temperature. 
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THE SPECTRUM THE 


ABSTRACT 


Two band systems attributed the molecule have been excited dis- 
charge tube. The first these system red-shaded bands with heads 
3489, 3556, 3626, and 3696 The second extensive system 
bands the visible and near ultraviolet region. High resolution spectra have 
been analyzed and the rotational and vibrational constants have been determined. 


INTRODUCTION 


The electronic structures silicon and carbon are similar and the Siz mole- 
cule might expected show spectrum similar that the well-known 
spectrum. The only published account spectrum attributed 
short note some flame bands which were photographed low resolution 
Downie and Barrow (1). Recently this laboratory two new band systems have 
been excited electric discharge and has been shown that these bands are 
due The method excitation and the analyses the bands are discussed 
this paper. 

EXPERIMENTAL 


Silane few millimeters pressure was admitted into discharge tube with 
aluminum electrodes and discharge passed through the gas. This process 
deposited grayish film (presumably Si) the walls the tube and the 
electrodes. The residue the SiH, was then pumped out and few millimeters 
pressure xenon admitted the tube. means small centrifugal pump, 
the xenon was circulated over heated mixture active uranium and uranium 
hydride (4) which served the dual purpose purifying the xenon and maintain- 
ing small constant partial pressure hydrogen the discharge tube. When 
discharge low current density was passed through the tube, the bands, 
described below, were observed. 

The spectrum was photographed first Hilger medium quartz spectro- 
graph and later 35-ft. concave grating spectrograph. Because the bands 
were weak, wide slit was used with the large spectrograph and the resolving 
power was limited about 150,000. This resolving power proved quite 
sufficient resolve the fine structure the bands. For the stronger bands 
exposure time one hour was sufficient but for the weaker bands even 
exposure time thirty hours proved insufficient. 

Besides the very strong bands the system SiH, which will 
discussed later paper, the low resolution plates showed four double- 
headed red-shaded bands with their short wavelength heads 3489.1, 3556.5, 
3625.9, and 3695.5 addition, number weak and almost headless bands 
could seen extending from 3980 into the red. Subsequent analysis the 
two band systems from high resolution plates showed that the stronger bands 
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Contribution from the Division Pure Physics, National Research Council, Ottawa Ontario. 


Issued N.R.C. No. 3791. 
801 


f 
i 


802 CANADIAN JOURNAL PHYSICS. VOL. 


are due transition Sig and the weaker bands are due 
transition the same molecule. 


THE BAND SYSTEM 


3489, 3556, and 3697 are similar appearance and was once apparent 
that these bands form lower state progression. The structure the bands is, 
however, most unusual since they show only two distinct heads, but large 
number branches, many which end abruptly. The fact that some the 
branches the bands show (see (2), 268) indicates that the 
bands arise from homonuclear molecule with zero nuclear spin. Close the 
strong band heads, two weaker heads can seen (see Fig. and the positions 
and intensities these heads leave doubt that they arise from the less 
abundant isotopes silicon. Thus was quite apparent that the bands arose 
from transition (or molecule, but, because the unusual 
appearance the bands, was not all apparent what type electronic 
transition was involved. was only detailed investigation the common 
combination differences the upper state the bands that was shown that 
each band possesses three and three branches and that they can readily 
interpreted bands resulting from transition. The peculiar appear- 
ance arises from the fact that only portions two the subbands are present 
and only the subband longest wavelength essentially complete. The wave 
numbers the lines the (0, 1), (0, 2), and (0, bands and their rotational 
numberings are given Table The excellent agreement the combination 
differences the upper state leaves doubt that this assignment correct 
and that the large number missing lines arises from some unusual sort 
predissociation. 

For state with Hund’s case coupling, the three spin components 
and can treated independently and each component can 
assigned effective value designated B*. The component with the largest 
effective value (that is, the component whose energy levels lie highest) 
designated the component. the band found that the subband 
lying longest wavelength has the largest effective values both the upper 
and lower electronic states and therefore involves the components. 

well known that, for homonuclear molecule with zero nuclear spin, the 
A-doubling the rotational levels state shows the bands 
transition staggering the lines the branches. For state 
(case a), the A-doubling largest for the component. the observed 
bands, the subband lying longest wavelength has much the largest staggering 
band and, therefore, most probably the subband. appears, there- 
fore, that both the upper and lower electronic states, the component 
the component or, other words, both states are inverted. 

The spin coupling constants cannot determined accurately from 
transition which only transitions with are observed. Rough 
values can found from the approximate equation 
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TABLE 
WAVE NUMBERS THE LINES THE SYSTEM Siz 

(0, Band 
02.34 
01.46 
00.51 
96.19 
94.98 
93.69 
92.33 
90.91 
44.01 
41.41 
38.91 
36.18 
33.59 
30.93 
28.24 
25.47 
22.67 
19.88 
16.99 


87.61 


98.3 


13.97 
94.76 13.59 
93.70 13.44 
92.28 13.00 
91.14 12.79 
89.67 12.26 
88.48 12.00 
86.89 11.41 
85.61 11.08 
83.97 10.36 
82.62 09.96 
80.95 09.20 
79.46 08.72 
77.69 
76.16 
74.30 
72.71 05.78 
70.76 04.77 
69.09 04.07 
28009 45.24 65.29 02.16 
44.01 01.02 
04.77* 42.85 00.13 
00.3 57.28 
95.72 37.61 53.03 95.53 
36.18 50.69 94.13 
90.98 34.75 
88.48 33.20 46.18 91.39 
31.68 43.99 89.91 
83.51 30.04 
80.95 28.44 
78.37 
75.74 
72.96 
70.14 
04.85 
02.06 
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TABLE (continued) 


WAVE NUMBERS THE LINES THE SYSTEM Siz (continued) 


(0, Band 
27896 .04 
92.84 
89.90 
86.59 
83.60 
80.20 
77.13 
73.61 
(0, Band 
65.80 
65.05 
64.25 
62.45 
61.50 
60.47 
59.41 27570.81 
58.28 70.60 
57.10 70.35 
55.86 
54.56 60.59 
53.17 59.40 
58.46 
57.18 
56.17 
54.85 
53.79 
49.75 
44.11 
42.76 
41.07 
36.40 
12.24 34.56 
11.33 33.01 
10.30 31.09 
29.45 
71.04 
07.20 25.78 
66.82 23.72 
64.89 04.94 21.94 
46.21 62.72* .67 19.81 
44.73 60.46 02.49 17.95 
01.18* 43.19 58.18 01.18* 15.75 
96.19 53.79 11.52 
51.51 97.13 
90.96 36.72 49.12 
35.00 46.75 94.21 
44.25 
82.95 31.48 41.7 
80.21 


76.26 
75.95 


804 
6 f 
75.95 
75.49 
75.49 
74.95 
74.79 
74.22 
74.01 
73.40 
73.09 
72.38 
72.03 
71.28 
70.81 
69.99 
68.59 
67.97 
65.25 
64.50 
63.40 
62.58 
61.38 
60. 
59.21 
58.18 
56.77 
27.76 
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TABLE (concluded) 


(0, Band 

30.96 

30.12 

29.26 
28.31 
27.35 
25.19 .51 
24.05 37.28 
22.87 
21.62 
20.24 
25.77 
43.21 
22.28 42.87 
21.31 42.87 
19.91 42.39 
18.89 42.39 
17.41 41.97 
16.36 41.81 
14.95 41.29 
13.71 41.08 
12.16 40.48 
10.91 40.17 
39.52 
39.19 
04.95 
37.26 
01.81 36.77 
00.04 35.95 
96.61 34.40 
19.59 33.77 
18.34 93.14 32.82 
17.02 91.46 32.09 
15.72 89.48 31.07 
12.95 70.99 85.70 29.14 
66.15 10.08 81.76 26.98 
63.79 21.81 67.45 79.89 
61.34 19.57 66.15 
58.85 17.41 64.84 
56.40 03.91 14.95 
53.78 02.29 12.72 61.84 
51.21 00.63 
48.59 
45.93 

Overlapped. 


value; seen but not measured. 


For the (0, band was found (by comparing values 46) that 
the difference the effective values the and components 
the lower electronic state and the upper. These differences combined 
with the observed values give 71.6 and 27.6. The value 
probably the best value that can obtained from the observed experimental 
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data. the other hand, sufficiently small that for longer 
have case coupling and equation not valid. more accurate value 
can determined finding value which, when substituted into 
equations (2, 235) for the energy levels, will give the observed differences 
between the and the component. When calculated 
accurately since not only this method finding quite inaccurate, but, 
the bands, predissociations the component make impossible 
compare very accurately the values the two components. The dif- 
ference the positions the origins the and subbands should 
given 2(A’— The observed difference the band origins 98.1 
while the value calculated from the above values 110.2 

The values both electronic states can best determined from the 
subband since here long and branches are observed. For the lower state 
where large, the effective value the component was found 
assuming 


From this value and from A”, the true value was determined Eq. 
For the upper electronic state where small, the true value was ob- 
tained directly fitting the observed values expression for 
calculated from the formula using the previously obtained value 
Bvalues were obtained from the (0, and (0, bands and was found 
from the differences the values. The measurements the band 
are less accurate. The rotational constants are given Table II. 


TABLE 
MOLECULAR CONSTANTS 


numbers are and values. 

vibrational numbering not known and these values are obtained assuming that the 
lowest observed level the level. 

¢The vibrational numbering this state not certain. 

The origins the subbands could determined quite accurately and 
fairly good values could also obtained for the origins the subband. 
The mean these two subband origins was taken the band origin. The 
vibrational constants the upper electronic state could not determined 
since only one vibrational level was observed. The lower state values 
could determined quite accurately for and from the (0, 1), 
(0, 2), and (0,3) bands. These data are sufficient determine and 
which are given Table II. The heads the isotopic bands are useful 
assigning the vibrational numbering and giving positive identification 
the emitter the bands, but are too weak measured accurately and 
therefore not give good value for 
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The observed abrupt terminations number the branches the bands 
are unusual and are not understood. certain that the upper electronic 
state responsible for the effect. probable that the effect arises from some 
type predissociation and below discussed such. The component 
least affected, but even this component the predissociation marked above 
47. overexposed plates, however, lines the branch arising from 
levels with between and can seen, but the branch breaks 
off sharply second time. the component, the predissociation affects 
all the low rotational levels and causes the subband headless. 
Also, all levels above are predissociated. The component 
these predissociations are probably the nature accidental predissociations, 
the fact that only one vibrational level the upper state observed seems 
indicate that the higher vibrational levels are completely predissociated. 

addition the striking predissociation, unusual intensity alternation 
occurs the lines the subband. This effect especially marked low 
values where every other line much weaker, but for greater than about 40, 
the effect almost disappears and successive lines have almost the same inten- 
sity. This effect may also arise from predissociation. state pre- 
dissociated state, then only one set the components affected (see 
(2), 417). For homonuclear molecules with zero nuclear spin such 
one set components accounts for alternate lines bands. Thus weak 
predissociation state may reduce the intensity alternate lines the 
component. 


IV. THE BAND SYSTEM 


The high resolution plates showed number simple bands extending over 
the region from 3980 into the red. large difference between the upper and 
lower state value causes the bands nearly headless. These bands are 
quite weak, and between 4090 and 4290 where the very strong SiH bands 
occur, bands this system could observed even though appears that 
number should occur this region. One the bands shown Fig. 

The bands have simple and branch structure, but the analysis shows 
that alternate lines (those even value) are missing. This structure implies 
that the bands are bands homonuclear molecule zero spin. low 
values the lines the bands break into two three components which are 
just barely resolved and again high values the lines become broad 
they are composed number unresolved components. This behavior 
exactly what might expected for bands transition. The values 
the upper and lower electronic states are about 0.17 and 0.24 respec- 
tively. The magnitude these values together with the type transition 
sufficient limit the emitter the bands Siz. The bands not involve 
either the well-known states and the strong bands not appear 
our discharge tube. Thus seems certain that the new bands are emitted 
the molecule. 

The rotational analysis the bands presented difficulty. number 
the bands had upper lower state common that the correct numbering 
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the lines was not doubt. The triplet structure the lines was not 
ciently well-resolved measured, and the three components are measured 
single line. The wave numbers the lines the bands are given Table 
III and the band origins and and values are given Table IV. 


TABLE III 


(m +4, n) (m +4, n+1) (m +1, n+2) 
P(N) R(N) P(N) R(N) P(N) R(N) P(N) 


-91 24618 .76 
-15 25120.99 618. 
89 A 5 


26 
23330. 66 -66 23098 .94 22825.87 22836. 

9 326.93 -60 095.24 820.86 83 
23309. 43 091.01 815.45 
1 317.91 809 .33 
312.60 080.88 802.79 
306.75 075.08 795.76 
300.38 4 068.71 788.20 
293.45 2.3 061.77 780.11 
286 .00 -46 054.32 771.45 
277 .98 24. 046.32 762.32 
: 037.84 752.64 
60.37 028.77 742.43 
019.18 731.69 


HE 08 .53 009.01 720.44 
8.51 24996. 8% 18.16 -96 29. 39.52 22998 .34 
5 987.12 
975.31 


(m +4, n +5) . : (m+1, n +4) (m, n +4) 
P(N) R(N) P(N) R(N) 


.29 

.26 22639 .08 22348 .39 22072.56 22081 .66 
5 635 .92 a 068. 078.73 
632.21 22558.81 ; 2. 063. 075.28 
627.90 548. 058. 071.33 
623.13 542.6: : 052.35 066.88 
617.77 536. 2. 28.36 046. 061.93 

611.92 529.25 22546.39 6 22.76 039. 
605.50 521. 298. 031.9: 050.48 
598.57 A : 024. 044.00 
591.05 26. 2.92 037 .04 
2 27% 5. 26 029.59 
021.63 


Though the relative vibrational numbering the bands readily obtained, 
the absolute numbering not all certain. Unlike the bands, the 
bands have strong band heads which show isotopic structure. Also, 
the lines the bands are quite weak and lines the isotopic molecules 
(which have only about one twentieth the intensity) have been observed. 
Towards longer wavelengths the bands become progressively weaker until 


N R(N) P(N) R(N) 
3 25118 
5 117 
7 114 
9 111.79 117.90 615.14 
ll 107.99 115.49 24605.23 612.76 Ld 
13 103.62 112.46 600.97 609.8 
15 098.68 108.85 596.15 606.3 24 
17 093.25 104.69 590.76 602.1 58 
19 087.17 099.92 584.78 597.5) 42 
21 080.51 094.58 578.23 592.3) 77 
23 073.28 088.64 571.12 586.4 58 
25 065.44 082.13 563.42 580.1 87 
27 057.04 075.00 555.16 573.1 60 
29 048.03 067.29 546.30 565.5) 82 
31 038.45 059.02 536.88 557.4 57 
33 028.25 050.13 526.86 548.7 73 
35 017.49 040.64 539.4 38 
37 006.15 030.58 529.5) 50 
39 24994 16 019.90 493 .37 519.0) 10 
41 748.13 
43 oF 737.71 
45 9 726.74 
47 715.21 
49 957 .73 458 .08 194.73 
51 181.87 
3 
5 
7 
9 
11 2263 
13 63 
15 62. 
17 62 
19 61 
21 
23 60 
25 59 
27 58 
29 
31 
33 
35 
37 476.82 502.09 253.38 278.46 
39 466.31 492.92 242.87 269.24 
41 455.24 483.20 231.83 259.53 
43 443.72 220.33 249.32 
45 431.63 462.20 208.21 238.56 
47 450.95 227 .32 
49 437 .62 
51 
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TABLE 
ORIGINS AND ROTATIONAL CONSTANTS THE BANDS 


23341 .23 

22646 

22088 


values were taken over from the better developed bands and used the determination 


from the observed lines. 
values were taken over from the better developed bands and used the determination 


they can longer observed. the violet end the band system, one 
the strongest bands (here designated the m+4, band) found 3980 
and bands are observed shorter wavelengths than this. The abrupt 
termination the band system the violet suggests that all levels above 
the level the upper electronic state are predissociated are not 
excited and that the level the level the lower electronic state. 
Thus probable that but the true value unknown. 

The vibrational and rotational constants the two states together with 


those the states are given Table 


DISCUSSION 


The electron configurations Sig may expected similar those 
the molecule which have been discussed Mulliken (3). Because the 
relative positions the energy levels change slightly the principal 
number increases from the lowest electronic configuration (see 
below) one which slightly excited and the resulting ground state 
expected state (see (2), 343). Since the nucleus obeys 
Bose statistics and since the even rotational lines the observed 
bands are missing, the bands must arise from either 
transition. seems reasonable therefore associate the observed lower 
state with the ground state. The fact that the electron configurations show 
‘that the ground state the only low-lying state and that 
transitions are expected from low-lying states gives strong support the 
assumption that the lower state the observed bands the ground 
state. accepted that the lower state the ground state, then the 
upper state can uniquely assigned definite electron configuration 
since there only one state arising from low-lying configurations. The 
observed bands can arise either from from transi- 
tion. transition can expected from low-lying states and 
therefore assumed that the bands are due transition. The electron 
configurations give only one low-lying state, but somewhat higher energy 
two states. Because the small change the values the observed 
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2.8 2.1 

3.7 2.1 

3.0° 2.0° 

3.0° 2.0° 

2.0° 

n+5 

m,n+3 2.0° 

3.1 2.0 

m,n+4 
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transition, probable that these bands are analogous the Swan band 
which also shows small value From these assumptions, the 
probable electron configurations the observed states are then follows: 


KKLL 
KKLL 
KKLL 
KKLL 


With our present limited knowledge the assignment the observed 
states the configuration given above not certain, but there appear 
other low-lying configurations which can account for the spectrum equally 
well. 

the lower state the bands the ground state, seems probable, 
then these bands may occur astrophysical sources. The large value 
for this transition causes the bands nearly headless and spread 
over large wavelength range. Therefore, even the bands occur such 
sources, they will difficult find. 


indebted Dr. Herzberg for numerous helpful suggestions and 
Miss Lila Howe for assistance the calculations. 
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HIGH RESOLUTION RAMAN SPECTROSCOPY GASES 


ABSTRACT 


The pure rotational Raman spectra allene, allene-d,, and were 
the second order ft. grating. Two plates each spectrum 
were analyzed yielding the following values for the rotational constants: 


the length the carbon-carbon double bond allene significantly shorter than 
that ethylene. 


INTRODUCTION 


Although the structure the allene molecule has been the subject several 
investigations electron diffraction and infrared spectroscopy, the geometri- 
cal parameters are not accurately known. Values for the distance and HCH 
angle have yet determined directly, and the existing values for the 
distance are not satisfactory agreement. The most accurate data available 
present were provided recent investigations the infrared spectra 
allene and allene-d,. These investigations, however, are not altogether without 
complications. Lord and Venkateswarlu (9) have carried out detailed study 
the perpendicular bands, and their results together with several plausible 
assumptions yield value for the HCH angle. Also, high resolution 
studies parallel bands Overend and Thompson (10) lead value 
1.30+0.005 for the length the carbon-carbon double bond allene. 
interest that this value considerably smaller than that found for the car- 
bon-carbon bond length ethylene (1.344 A), thereby suggesting that the 
length the carbon-carbon double bond dependent environment, 
conclusion already well established for the single bond (see, 
for example, Ref. 5). Since this result considerable importance for valence 
theory, investigation the rotational Raman spectrum allene was 
undertaken. 

The pure rotational Raman spectra allene, allene-d,, and 
were photographed with high resolution. The rotational lines were well resolved 
and the analysis was straightforward, resulting accurate evaluation the 
C=C distance. This bond length shorter than that found ethylene 
agreement with the observation Overend and Thompson. Also, relation- 
ship between the distance and HCH angle was obtained, although was 
not possible evaluate these parameters independently. 


Contribution from Pure Physics, National Research Council, Ottawa, Canada. 
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Allene was prepared the dehalogenation pure 2,3-dibromopropene 
with zinc dust alcoholic solution (7). and mol. 
per cent and mol. per cent purity respectively were prepared new 
synthesis Leitch and Morse the Division Pure Chemistry (8). the 
present investigation, the gas pressures were cm. for allene and cm. 
for allene-d, and 

The pure rotational Raman spectra the three gases excited the 
4358 line were photographed with the apparatus described this series 
(14). ft. concave grating spectrograph was used giving dispersion 
6.75 The slit width was 0.03 mm. Exposure times were the order 
15-20 hr. with Eastman Kodak 103 a-0 photographic plates. Second- and 
third-order iron lines provided the necessary wavelength standards. 


OBSERVED SPECTRA 


Photographs the pure rotational Raman spectra allene and allene-d, 
are shown Fig. Approximately lines are observed each side the 


TABLE 


OBSERVED WAVE NUMBER SHIFTS THE ROTATIONAL RAMAN LINES ALLENE, 
AND 


Allene Allene Allene-d, 
55.072 


values are averages the Stokes and anti-Stokes lines measured two plates, except the 
values marked 

lines measured one plate only. 

*Lines blended with grating ghosts. 


| 
‘ 


q 
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exciting line both spectra with only slight overlapping grating ghosts. 
The main series lines the branch corresponding transitions 
The line spacing about 1.2 and 0.9 for allene and allene-d, 
respectively and can seen that the rotational lines are clearly resolved. 
few lines midway between the main series lines can also seen the 
neighborhood the exciting line and these are part the branch (AJ 1). 

The rotational Raman spectrum similar appearance 
the spectra shown Fig. that is, the spectrum simple and the lines are 
sharp and clearly resolved, having spacing about The slight 
asymmetry the molecule did not produce any noticeable com- 
plications the observed spectrum. 

Two plates each spectrum were evaluated the usual way (15) and the 
observed wave number shifts for the rotational lines allene, and 
are listed Table Overlapping grating ghosts interfered with 
the measurement only few the lines and these are indicated the table. 
For all three molecules, the values obtained from two plates agreed 
within +0.02 


ROTATIONAL ANALYSIS 


(a) Allene and 

The allene molecule has the symmetry Va) and hence prolate 
symmetric top. The Raman selection rules for rotational transitions sym- 
and and branches both Stokes and anti-Stokes sides the undis- 
placed line. the spectra the allenes, however, only small fraction the 
total intensity the branches, and hence need only consider the 
branches for the present analysis. discussion the intensity distribution 
included the end this section.) 

The wave number displacements the rotational lines the branches 
are given 


where the centrifugal stretching constants and are small compared 
the rotational constant The line spacing very nearly 4Bo, 
and when for any line divided the average value the observed 
spacing, the rotational numbering unambiguously determined. Values the 
constants were obtained plotting |Av|/(J+3/2) against The 
graphs are shown Fig. From the intercept, one obtains 
and from the slope, The value 6D, entirely negligible 
otherwise the rotational lines would have substructure least would 
appreciably broadened. Hence the intercept essentially The values the 
constants obtained for allene are 


and for allene-d, 


814 


Fic. 


The values give the effective moments inertia for the ground states: 
gm. cm.? for allene, and 


using the atomic constants DuMond and Cohen (1). 
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x 


graph for the ground states allene and allene-d,. lines blended 
with grating ghosts.) 


cm: 
+0.04 
Allene 
+0.02 
fe) 
-0.02 
x x 
—0.04 
x 
x 
+0.02 x °° 
° 


307 


Deviations the observed wave number shifts the rotational lines from those 


calculated with the derived constants. lines blended with grating ghosts.) 


J+3/2 
1.1870 
x 
x (e) 
1.1860 
x 
Allene 
1.1850 x COeO., x 
x — 
0.930 
0.928 
0.927 
x 
° 10 20 30 40 
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The rotational constants were then used calculate the wave number dis- 
placements the rotational lines and comparison with the observed values 
shown the graph Fig. Almost all the observed —calculated values 
lie within +0.02 agreement with the estimated accuracy +0.03 
The scatter the values for allene-d, larger than for allene and 
probably due the presence small amounts partially deuterated mole- 
cules. 

Although calculation from force constants was not carried out, was 
tempting apply formula the type When the and 
values just derived were substituted this equation, the following values were 
obtained: 

974 for allene, and 818 for 


These are within 10% the observed frequencies 
1076 and 874 (9). 

The values the rotational constants derived here differ slightly from those 
obtained from the parallel infrared bands Overend and Thompson (10), 
d,, and for both. These differences may arise from complica- 
tions the infrared spectra which not permit clear resolution the lines, 
observed the Raman spectra, more probably may due the 
difficulty determining the correct rotational numbering such closely- 
spaced spectra. 

The observed intensity distribution the rotational Raman spectra the 
allenes quite regular with continuous decrease intensity the outer 
extremities; there intensity alternation observed the parallel 
infrared bands allene (10). This difference appearance readily explained 
considering the temperature the gas samples the two experiments. 
the latter study, the gas was cooled that large fraction the 
molecules would the rotational levels with where the statistical 
weight factors are the ratio 7:3 for even:J odd (3). the present ex- 
periment, the temperature was about 50°C. and rotational levels with 
were appreciably populated. Since intensity alternation occurs for the 
subbands with intensity alternation predicted for the Raman 
spectrum. Also, must pointed out that molecules which have been ther- 
mally excited into the lowest vibrational levels will contribute the observed 
Raman spectrum and affect slightly the intensity distribution and line posi- 

Placzek and Teller (12) have calculated the intensity distribution the 
rotational Raman spectrum symmetric top molecule and have shown that 
the intensity the branch decreases with decreasing value the ratio 
the least and largest moments inertia For allene, this ratio less 
than 0.1 and the intensity the branch would expected extremely 
small. can also shown from the equations Placzek and Teller that the 
branch lines with will have the highest intensity. Both expectations 


qualitative discussion given Section paper this series (15). 
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are borne out here; only few lines the branch can seen the spectra 
allene and allene-d, and these are fact confined the neighborhood 
the exciting line. 


(b) 

The isotopic molecule has the symmetry and has slightly 
asymmetric structure. However, since the general appearance the rotational 
spectrum simple, this molecule must very nearly symmetric top and 
hence was assumed that the wave number displacements the branch 
lines can expressed equation similar [1] above, but with replaced 

constants reproduce the observed displacements within +0.02 


MOLECULAR STRUCTURE ALLENE 


The structure and the structural parameters the allene molecule are shown 
Fig. The moment inertia about any axis perpendicular the 
axis and passing through the center atom given the simple expression 


where Mp. When the numerical values for the moments inertia 
allene and are substituted equation [2], one obtains 


and 


Fic. The structure and structural parameters the allene molecule. 


The latter value gives relationship for the distance, and the 
HCH angle, 26, namely 


Equation [3] represented graphically the shaded strip Fig. 
The width the shaded strip gives estimate the error the values the 
parameters based the above experimental errors. 

Initially had been hoped that the investigation would 
provide sufficient information determine and HCH independently. 


A =C 20 “A 
4 
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1-40 
08 
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noe use 


Fic. Graphical determination the distance and HCH angle allene. The 
shaded band B-B represents the values that are consistent with the values for allene and 
represents the values consistent with the infrared value (9). 


the graph (Fig. this would represented line crossing some 
point, giving the values these parameters consistent with the value Bo. 
However, the present data give line which falls entirely within the shaded 
area thus confirming the right-hand side Eq. [3], but without providing 
further useful information. This result can understood one recalls that 
Eq. [3] evaluated from the difference the moments inertia allene and 
find the second difference the moments inertia the three isotopic 
molecules. For the three allene molecules, the second difference very small 
number and has rather large uncertainty. With the present data, one can 
ethylene (2)) and obtain the value ZHCH 117+1.5° from the graph, 
Fig. 

summary, the results the present investigation the rotational spectra 
the allene molecules lead the value 1.309 and are consis- 
tent with the values 1.07 and ZHCH 117°. 


DISCUSSION 


The value for the carbon-carbon distance allene found here, 1.309 
good agreement with the value 1.30+0.005 found Overend and Thomp- 
son (10) from the infrared spectra. While the electron diffraction result 
Wierl (16) for this distance, also agreement, the value 
quoted Pauling and Brockway (11) for allene and ethylene 
not accord with the present results. shown Overend and Thompson 
and fully confirmed the present results, the double bond 
allene significantly shorter than that ethylene.* show the dependence 


value for the carbon-carbon distance ethylene was calculated from the most 
recent Raman value (13) und the infrared value (2). 


; 
= 
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the bond length environment, one may write: 


1.309 


confirmation this short bond length comes from the microwave spectrum 
ketene. Johnson and Strandberg (6) found the value 2.475+0.003 for the 
sum the carbon-carbon and distances. ketene the en- 
vironment the bond similar that and one assumes the 
value 1.162 found Herzberg and Herzberg (4) for 


value for the bond ketene obtained. Thus 


the spectroscopic investigations allene and ketene have established 
marked dependence the length the double bond adjacent 
bonds (5). this connection, would especially desirable investigate 
the structure butatriene which has one double bond with 
environment different from that allene, namely bond length 
considerably shorter than that allene would expected for this bond. 

The only evaluation the least moment inertia, for the allene mole- 
cule that obtained Lord and Venkateswarlu (9) from analysis the 
perpendicular infrared bands. From their value one 
this relation included Fig. can seen that their value and 
the new are consistent with the values and 
ZHCH 115+1°. Both values appear slightly low. definitive value 
the distance and HCH angle must await further high resolution studies 
the infrared‘ and Raman spectra allene. 


grateful Dr. Herzberg for helpful discussions and Dr. Leitch 
and Mr. Morse for preparing the deuterated allenes. 
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SOME EXACT SOLUTIONS THE LORENTZ INVARIANT 
PROBLEM THE MOTION TWO ELECTRIC FLUIDS' 


ABSTRACT 


Methods generating exact steady-state solutions for the flow super- 
imposed positive and negative fluids are developed for the cases linear and 
cylindrical motion. two fluids are assumed subject body forces 
arising from the total and fields. shown that the only axially sym- 
metric solutions are cylindrical symmetry, i.e. rotating spheres rings 
charge cannot satisfy the equations where the two fluids overlap. 


The problem the dynamics ionized media has been attacked various 
methods arises different fields investigation. This paper approaches 
the general problem not attempting find solutions with initial 
boundary conditions, but seeking any exact solutions the Maxwell 
equations coupled the dynamical equations motion. The resulting solu- 
tions may interpreted limiting cases ion cloud motion, may 
interpreted the sense the electrodynamics” Buneman 
(1). 

idealize the problem representing the positive and negative charges 
positive and negative perfect electric fluids, subject only body forces 
arising from the total fields and and not subject forces arising from 
pressure gradients viscosity. The general procedure assume some 
velocity field simple symmetry. From the Maxwell and dynamical equations 
then find the conditions the velocity components which give possible 
solutions. These conditions usually resolve into differential inequalities, be- 
cause the charge densities the two fluids computed from V-E and 
must essentially positive essentially negative quantities. 


GENERAL EQUATIONS, M.K.S. SYSTEM 


Symbols 

velocity positive fluid. 

=velocity negative fluid. 
electric field. 
=magnetic induction. 
scalar potential. 
charge density positive fluid. 

charge density negative fluid. 
rest ratio positive fluid. 
rest ratio negative fluid. 


received original form March 30, 1954, and, revised, August 1955. 
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electric permittivity free space. 
magnetic permeability free space. 


Lorentz invariant form the dynamical equations governing the steady- 
state motion the two fluids subject only body forces due and can 
put the form (2) 

div(pe) 


The Maxwell equations are 


The continuity equations and are not independent one another, 
since the conservation total charge follows from the Maxwell equations. 


II. STEADY-STATE LINEAR FLOW 
the simplest example consider that and have only components. 
For brevity write for and for this section. The dynamical 
equations give 
a 0a/dx = PE, a4, B 0B/dx —WNE, Ba, 


Then otherwise the and components and all vanish. With 
most constant, which can taken zero this problem. Then 
terms the scalar potential 


Substitution the div equation 


with the solution F(a). 
The conditions far imposed satisfy all equations except the two Maxwell 


equations now the form 


(curl 
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The first these solved for and substituted the second. The result- 
ing equation then 


The substitution g(a) gives 


can now generate any number solutions the steady-state linear flow 
problem. First, select any solution the two-dimensional Laplace 
equation; next, select any function F(a). This determines a(y, and 
From these the fields and charge density are derived. the linear 
steady-state flow problem only appears; reference made and 
individually. 

One other type non-accelerated linear flow was examined. was assumed 
that the flow was parallel straight lines any plane perpendicular the 
z-axis, but the direction the parallel flow lines might vary from one such 
plane another continuous way. was found that such solutions are 
possible with any actual variation flow direction with 


III. STEADY-STATE AXIAL SYMMETRY SPHERICAL COORDINATES 
possibility the flow the superimposed positive and negative fluids 
circles about the z-axis; that is, assume The remaining 
are zero. The dynamical equations reduce 


From these equations 


first alternative leads the trivial solution The second alternative 
states that cylindrical coordinates (r, and that therefore 
complete cylindrical symmetry the only symmetry possible. cannot, for 
example, have rotating spheres overlapping charge, polarized toroidal 
rings rotating about the axis. 


IV. CYLINDRICAL SYMMETRY 


circles around the z-axis, Bz, and are all zero, well 
all The remaining and are then designated simply and The 
dynamical equations become 
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and the only necessary Maxwell equations are 
d(rE,)/dr 


with all quantities functions only. Eliminating and means the 
dynamical equations and solving for and have and expressed 
non-linear functions and and their derivatives. The condition now 
imposed that and chosen such functions make and every- 
where positive. Solutions this problem have been obtained two special 
cases. 


Set and 1/P+1/N and consider the case which 
positive. The expressions for and become 


Addition the two inequalities gives (1/a) d(a/r)/dr and this condi- 
tion satisfied the inequality for also seen satisfied, since both terms 
are then positive. The equation for may written 


The denominator the left side always positive. Then taken any 
positive function and the resulting differential equation solved for the 
inequalities for and are simultaneously satisfied, because (1/a) d(a/r)/dr 
then positive. All solutions obtained this special case have con- 
stant magnetic fields along the addition the fields produced the 
moving fluid itself. 


Case II. 


The resulting equation and 


a 
q = 
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from which follows that 


log and integrate get allowed These solutions also show con- 
stant applied magnetic fields. 

Note that any exact solution for the fields and dynamical variables may now 
subjected Lorentz transformations, the cylindrical flow, e.g., going 
helical flow. 

Though the solutions far obtained not appear describe any usual 
physical situations, felt that the methods used may extended situ- 
ations less simple symmetry, and any case such solutions may used 
zero order solutions perturbation problems, exact solutions against 
which check approximation methods. 

The solutions given are all steady state, but the same analysis may ap- 
plied the general time-varying case, and this problem now being worked 
on. 


The author wishes express his appreciation for the assistance Mr. 
Bedient some parts this work. 
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EQUATIONS MOTION GENERAL RELATIVITY! 


ABSTRACT 


has been shown that both the equations motion charged particle 
gravitational field and the field equations can obtained from one variational 
principle suitably generalizing Dirac’s classical theory electrons. 


INTRODUCTION 


addition the field equations general relativity, which connect the 
distribution matter and energy with the resulting gravitational field de- 
scribed the potentials g,,, the motion particle this field must 
described additional equations, the equations motion. Although has 
been shown the case weak fields (2, that the latter can derived from 
the former, complete union has not been achieved. This least partially 
due the fact that two separate sets variables have employed. The 
field equations are functions the g’s and their derivatives, while the equa- 
tions motion employ co-ordinates (and their derivatives with respect 
parameter) themselves. This situation fundamentally different from other 
field theories, say, the scalar meson field, where the same set equations 
play the role field equations and equations motion and therefore only one 
set field variables necessary. 

Although both, the field equations and equations motion, can obtained 
from appropriate variational principle, has not been possible use only 
one describe both. Recently Dirac (1) has shown how Lorentz’s equations 
motion and field equations can both obtained from one varia- 
tional principle for stream electrons characterized velocity field 
and this note his method generalized for general relativity. 


FIELD EQUATIONS AND EQUATIONS MOTION 
field equations are given 
where R,, are the components the curvature tensor 


and 7,, the components the energy-momentum tensor. For the electro- 
magnetic field the latter are given 


where the F,, are the components the field tensor, which can defined 
terms the electromagnetic potentials follows: 
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addition, the satisfy Maxwell’s equations, which the absence any 
charge and current distributions are 


[4] 


Here denote tensor densities defined The other set 
equations automatically satisfied writing the field the form 
[3]. Both the field equations and Maxwell’s equations can obtained from 
variational principle (5): 


where for have take 


where are the Christoffel symbols, and for 


Variation with respect then results the field equations [2] where each 
quantity expressed tensor density, while variation with respect 
gives equations [4]. 

Similarly the equations motion free particle are given 


which can derived variational principle 


where use has been made the definition the line element ds? g,, dx’. 
the following will necessary use the covariant form the equations 
motion, which given 
this point might interesting point out that the equations motion 
can expressed Hamilton’s canonical form (4), might have use 
them that form when going the Hamiltonian form the theory. This 
achieved defining scalar inertial function which will become 
numerically equal the mass the particle consequence the equa- 
tions motion. Defining the equations 


consequence dx, dx,), obtain the canonical equations 
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LORENTZ’S EQUATIONS MOTION 


Lorentz’s equations motion for charged particle moving with velocity 
tional field are 


their covariant form 


and for stream continuous function and thus 


where 


From here follow Dirac’s treatment (1) and finally obtain for the 
potentials 


THE ACTION PRINCIPLE 


are now the position write down the variational principle, which 
give both the equations motion and the field equations. Let take 
for the action 


where for use the function previously defined (Eq. [6]) and the F,, are 
derived from the potentials [18]. independent field variables take 


[20] 


over the whole obtain 


*The difference sign due our alternative definition F,, (see Eq. 
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The condition that shall stationary after making use [18] gives 


Equation [22] becomes Maxwell’s field equations identify the charge- 
current density 


[25] 


Equations [23] become the field equations for particles electromagnetic and 
gravitational fields after using equation [21],* and, finally, equations [23] may 


and give, with the help [18], 


which are identical with [16] and thus lead back Lorentz’s equations 
motion. Thus the action principle gives both the equations motion and the 
field equations correctly. 

subsequent paper shall discuss the Hamiltonian form the field 
equations and give the corresponding action principle. 


conclusion, wish thank Professor Dirac for valuable dis- 
cussions. 
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plays the role the mass and thus identical the inertial function defined Section 
Hence becomes 
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GAMMA RAYS FROM THE DEUTERON BOMBARDMENT 
BORON 


ABSTRACT 


Bombardment with 1.4 Mev. deuterons has been found produce 
gamma rays energy 4.46+.04, 4.75+.03, 6.52+.04, 
6.78+.07, 7.29+.04, 8.27+.09, and 8.87+.02 Mev., measured with three- 
crystal pair spectrometer. Tentative assignments transitions and 
agree well most cases with the known energy level schemes. gamma rays 
very high energy were observed, even with single-crystal spectrometer; con- 
sequently upper limit has been placed the cross-section 
0.95 Mev. bombarding energy. 


INTRODUCTION 


Bombardment with deuterons produces simultaneously the mirror 
nuclei and through the reactions B!°(d, and 
the charge symmetry postulate that and forces are equal except for 
electrostatic forces correct, and many-body forces are either negligible 
charge independent, then the energy level diagrams mirror nuclei should 
identical the energies, spins, and parities the states. The ground states 
should differ energy the mass difference and the Coulomb energy 
the proton according 


where 
neutron mass, 
proton mass, 
proton charge, 
1.45 cm. the nuclear radius, and 


mass number. 


The energy level diagrams and (2) are shown Fig. with the 
ground states the same height facilitate comparison the level spacings. 
From equation [1] the ground state should Mev. above that 
the experimental energy difference 1.98 Mev. (2). This means 
indication that the charge symmetry postulate incorrect, since the electro- 
static term equation derived the assumption uniform charge 
distribution throughout the nuclear volume; furthermore, the nuclear radius 
not well-defined quantity. 

Unfortunately, there yet insufficient information, particularly 
spin and parity assignments, permit detailed comparison the two level 
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SAMPLE AL.: DEUTERON BOMBARDMENT BORON 829 


812 


Fic. Energy levels and Transitions shown were found the present work. 


schemes, but the pattern levels certainly similar. The first five excited 
states show marked resemblance each other spacing, with those 
lying the average 250 kev. below the corresponding levels The only 
certain spin assignment those shown Fig. that the ground state 
(8). There is, for example, considerable difference between the spin and parity 
assignments listed Ajzenberg and Lauritsen (2) and those listed Jones 
and Wilkinson (11). The energies the levels are assigned largely from the 
measurement energies proton and neutron groups from deuteron bom- 
bardment (7, 10, 14). Because the energy resolution attainable poorer 
for neutrons than for protons, possible that some levels are not yet resolved 
for example, the 6.46 Mev. level probably corresponds the 6.76, 
6.81 Mev. pair levels 


Yo + 
7-39 
5-03 
446 So 
214 
2 
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Proton groups corresponding all the known energy levels 
10.32 Mev. have been analyzed magnetically with only kev. uncertainty 
energy Van Patter al. (14) and Elkind (7). The relative intensities the 
proton groups measured 90° for 1.51 Mev. deuteron energy are listed 
Table The energies the neutron groups were measured with uncertainty 
kev. Johnson (10), who used the photographic emulsion technique. 
bombarding energy 3.6 Mev. found groups corresponding all the 
levels shown Fig. His values for the relative intensities the 
neutron groups observed the direction the deuteron beam are given 
Table The energy dependence the cross-sections both reactions 
basically non-resonant character, and the angular distributions some 
the groups have shapes indicating Butler stripping (2). 

Previous investigations gamma rays originating from the reactions 
showed only five gamma transitions (1), although moderate bombarding 
energies many more are possible. three-crystal pair spectrometer has been 
used the present experiment attempt determine the gamma-ray 
spectrum more completely, and obtain information the relative 


intensity the transitions. The magnetic pair spectrometer measurements 


Bent al. (4), reported since the experiment described here was completed, 
are listed Table for comparison. 


THE THREE-CRYSTAL SPECTROMETER 


With crystals having dimensions the order two inches and gamma rays 
the energy range 2.5 Mev., the pulse height spectrum from sodium 
iodide single-crystal spectrometer corresponding single gamma energy 
contains three peaks because the escape from the crystal one both 
quanta from annihilation positrons produced pair events. While the 0.51 
Mev. spacing the peaks may useful energy calibration, spectrum 
containing several gamma rays may difficult interpret because the 
multiplicity overlapping peaks. neutrons intensity similar that 
the gamma rays are produced the same reaction, large background tends 
obscure the gamma-ray peaks. This due largely neutron capture 
the crystal, which results the prompt emission gamma rays the 
energy range Mev., followed 2.0 Mev. beta decay, 1.6 Mev. beta 
decay coincident with 0.4 Mev. gamma ray. The beta decay half-life 
min. This background can cause, through pulses magnitude 
corresponding Mev. gamma rays. typical single-crystal spectrum 
obtained with in. diam. crystal from deuteron bombardment 

The three-crystal pair spectrometer, the expense factor 100 1000 
efficiency, completely eliminates the first difficulty mentioned above, and 
reduces the second level tolerable most experiments. The principle 
the three-crystal spectrometer has been discussed several authors (3, 15). 
The spectrometer used the present experiment differs detail, but not 
principle, from those described the above papers. Figure block diagram 
the complete spectrometer, excluding the high voltage and supplies. 


CENTER CRYSTAL 


SAMPLE DEUTERON BOMBARDMENT BORON 831 


MEV. 


COUNTS PER CHANNEL THOUSANDS 


PULSE HEIGHT, VOLTS 


14-06 20-26 


14 
CHANNEL NUMBER 


Fic. Single-crystal spectrum from 0.95 Mev. 230 microcoulombs. 
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FOLLOWER FOLLOWER MIXER 
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Fic. diagram the three-crystal spectrometer. Power supplies are not shown. 
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Differential discrimination each side channel was used reduce the num- 
ber unwanted true coincidences due the effect (12) 
and scattering single Compton quanta from one side crystal the other. 
Accidental coincidences were also, course, decreased the differential 
discriminators, and were further reduced the use triple, rather than 
double, coincidence mixing. addition, the side crystals were shielded from 
the target in. lead, and all three crystals were shielded from the back- 
ground X-rays from the Van Graaff generator in. lead. 

The side channel discriminators were adjusted pass the total energy pulses 
0.51 Mev. quanta viewing the spectrum the annihilation radiation 
from source with oscilloscope receiving pulses from the collector 
the side channel photomultiplier. With the oscilloscope trace triggered the 
output from the side channel discriminator receiving the corresponding last 
dynode pulse, the discrimination levels were adjusted until only the bright 
full energy group pulses was visible the screen. The one microsecond 
resolving time the triple coincidence mixer was sufficiently short make 
accidental coincidences negligible, since the permissible counting rate was 
restricted the effect high counting rates the gain stability the center 
channel photomultiplier (6). The coincidence output pulse triggered uni- 
vibrator which produced v., gating pulse. Center channel pulses 
were amplified Northern Electric No. 1444 amplifier and passed through 
the gated biased amplifier channel Marconi kicksorter. 1.5 
delay, the form 1000 ohm cable, was inserted the center channel 
locate the pulses within the gated interval. 
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The performance the spectrometer has been very satisfactory; its resolu- 
tion was found about Mev., shown Fig. energy cali- 
bration spectrum the 6.13, 6.94, and 7.1 Mev. gamma rays from the 0.873 

resolution 6.5% was obtained for the 2.62 Mev. line from radiothorium, 
shown dotted Fig. Comparison Fig. with Figs. and shows that the 
neutron rejection the spectrometer was very good indeed. The pulse height 
scales the various spectra not agree because different injection points 
the center channel circuit for the calibration pulses from standard pulse 
generator were tried, avoid troublesome slow drift the gain the center 
channel. the final configuration, although only the center channel crystal 
and photomultiplier were outside the calibration loop, gain drifts persisted. 
correlation between instability gain and counting rate was found, agree- 
ment with Reference (6). 


EXPERIMENTAL PROCEDURE 


The targets, kindly supplied the Isotopes Division, A.E.R.E., Harwell, 
were enriched better than 95% Deposited platinum backings, they 
were 100 thick, approximately kev. thick 1.5 Mev. deu- 
terons. They were clamped copper support insulated permit measure- 
ment the beam current. coolant was used, since some heating was desir- 
able reduce carbon build-up the targets. 3/16 in. diam. beam-defining 
stop cut from 1/16 in. molybdenum sheet was placed in. from the target. The 
spectrometer was located the deuteron beam, the distance from the 
target the center the center crystal being cm. 

inserting standard pulses some point the center channel, the later 
runs the grid the cathode follower, and adjusting the gain and bias the 
Northern Electric amplifier-biased amplifier combination, the desired range 
pulse heights was made cover the kicksorter channels. This setting was 
checked the beginning and end each two hour bombardment. Gamma 
rays known energy were used for energy calibration, usually those from 
(6.13, 6.94, and 7.1 Mev.) and radiothorium (2.62 Mev.). This 
was repeated least twice each day, two more runs intervening between 
calibration runs. check the linearity the center channel gain, several 
other energy points were used occasionally. The center channel could 
operated single-crystal spectrometer disabling the gate, permitting 
determination the photoelectric peak the 2.62 Mev. line from radio- 
thorium, and those the 0.51 and 1.28 Mev. lines from This was not 
repeated often, the linearity was always very good, least far Mev. 
Checks with the 17.6 and 14.8 Mev. gamma rays from the 440 kev. resonance 
indicated that the gain was linear 17.6 Mev., but the peaks 
were broad and asymmetrical because bremsstrahlung losses 
effects the rather small center crystal. 

Even though power supply and heater voltages were carefully stabilized and 
monitored, and the center channel counting rate was kept less than 8000 counts 
per second limiting the beam about five microamperes, small gain drifts 


834 CANADIAN JOURNAL PHYSICS. VOL. 


persisted. analyzing the results, the peaks particular spectrum were 
abnormally broad, the spectrum was not used determining thé gamma-ray 
energies. the peaks spectrum were normal width, but the energies 
all gamma rays determined from the peaks deviated the same direction 
from the respective mean values other determinations, the energy scale was 
the mean value for the energy the most prominent gamma 
ray present the spectrum; this procedure was only necessary for one spec- 
trum retained the final analysis. 
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SAMPLE DEUTERON BOMBARDMENT BORON 
RESULTS AND DISCUSSION 


Typical spectra, obtained the 1.4 Mev. deuteron beam and covering 
the range gamma-ray energies from 9.5 Mev., are shown Figs. 
and The peaks are labelled according actual gamma-ray energy, 
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determined from each spectrum, rather than pair peak energy. The integrated 
beam current for each run indicated the captions the figures. Direct 
comparison intensities between the various spectra only rough, because 
the loss target material during bombardment. 

The energies, found averaging the results from several spectra, and 
relative intensities the gamma rays assigned transitions and 
are listed Table The assignments are indicated arrows Fig. The 
energy errors listed the table are standard deviations several determina- 
tions each gamma energy, ranging from three determinations the 4.75 
Mev. gamma ray twelve the 7.29 Mev. gamma ray. The measurements 
Bent al. (4), listed Table for comparison, were corrected for Doppler 


TABLE 


|EUTRON 


shift, whereas the present measurements were not corrected. Doppler correc- 
tions for center mass motion (for 1.4 Mev. bombarding energy and observa- 
tion 0°) range from kev. for the 4.46 Mev. gamma ray kev. for the 
8.87 Mev. gamma ray, the order the standard deviations. However, 
several the neutron and proton groups are emitted predominantly forward 
(5, 10, 13), tending reduce the mean Doppler correction the corresponding 
gamma rays. 

The relative intensities listed Table were obtained measuring the 
areas under the peaks the spectra, subtracting somewhat arbitrary back- 
ground, normalizing the integrated beam current, and applying correction 
factor for the variation the pair cross-section with energy. correction, 
the order 30% for the higher energy gamma rays, due bremsstrahlung 


6.87 16.98 
4°75 4:74 
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loss and wall effects the center channel crystal, has not been applied, since 
the uncertainties intensity were this order. Agreement between successive 
runs was not good, owing loss target under bombardment and possible 
variations the distribution the deuteron beam the target, the inten- 
sities from overlapping spectra were normalized the use gamma rays 
common both. All intensities Table including those the proton and 
neutron groups, are normalized directly comparable. The gamma-ray 
intensities Bent (4) were measured for thick targets and Mev. bom- 
barding energy; the difference conditions perhaps explains why there only 
qualitative resemblance the intensities from the present experiment. 


TRANSITIONS 


The intensities the proton groups listed Table were measured 90° 
and 1.51 Mev. bombarding energy (14); because differences the angular 
distributions the groups, only qualitative comparison with the gamma-ray 
intensities possible. The difference the intensity the proton group and 
gamma ray corresponding the 4.46 Mev. state may explained 
the proposed cascade through this state from the 9.19 Mev. state. Serious 
discrepancies arise for the three highest states 8.92, 9.19, and 9.28 Mev. The 
8.92 Mev. ground state transition much weaker than the proton group 
exciting it, the 9.19 Mev. state decays apparently only 4.73 Mev. cascade 
gamma ray much less intense than the corresponding proton group, and 
gamma rays were attributed the 9.28 Mev. state. few early runs, with 
wider side channel windows, and thus greater efficiency, there was some indi- 
cation radiation energy greater than Mev., but these suffered from poor 
resolution and gain drifts. The proposed decay scheme does not agree 
well with that Jones and Wilkinson (11) from measurements the reaction 

The 8.27 Mev. gamma ray has been assigned the 8.57 Mev. state 
although the energy discrepancy three times the standard deviation the 
energy measurements. This discrepancy may arise part from the low inten- 
sity the radiation, and perhaps from conflicting counts due 8.87 Mev. 
photons which fail produce the full pulse height because bremsstrahlung 
loss and wall effects, and produce peak below the pair peak. similar effect, 
though much smaller proportion, has been observed with the clean 6.13 Mev. 
gamma ray from 


TRANSITIONS 


The neutron group intensities Table were measured and 3.4 Mev. 
bombarding energy (10); the difference bombarding energy and the strong 
forward peaking some the groups make comparison with the gamma-ray 
intensities rather hazardous. has been found this laboratory (13) that, 
1.5 Mev. deuteron energy, the neutron groups corresponding states 
6.46 and 4.77 Mev. are much more intense than other groups, with the 
possible exception the ground state group which was not under observation. 
the 4.75 Mev. gamma ray observed has been attributed part the ground 
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state transition from the 4.77 Mev. state well the 4.73 Mev. 
transition between the 9.19 and 4.46 Mev. states The 7.01 Mev. gamma 
ray attributed Bent al. (4) the 6.87 Mev. state was not observed 
this experiment. Mev. gamma ray was assigned the 5.49 Mev. 
transition between the 7.39 and 1.90 Mev. states analogous the 5.87 Mev. 
transition from the 7.99 Mev. state reported Bent al. (4), although 
comparison neutron and gamma-ray intensities not encouraging. 

Four peaks corresponding gamma-ray energies 2.57, 3.09, 3.52, and 6.13 
Mev. have been observed (Figs. and 7). The 2.57 Mev. gamma ray may 
from the 2.38 Mev. transition between the 9.19 and 6.81 Mev. states 
Carbon accumulated the target and beam stop was responsible, through the 
reaction for the 3.09 and 3.52 Mev. gamma rays, since their 
intensities were much smaller from new target and freshly cleaned stop. 
The 3.52 Mev. gamma ray, probably from the 3.68 Mev. state C!*, was seen 
only with very dirty targets, and had about one-fortieth the intensity 
the 3.09 Mev. gamma. The excitation curve these gamma rays also showed 
resonance near Mev. deuteron energy, agreement with the assignment 
The 6.13 Mev. peak may attributed the intense 6.12 Mev. gamma 
ray from the reaction Using the yield values given Reference 
(2) and assuming isotopic abundance for the 3.09 Mev. gamma-ray 
peak should have intensity larger factor 170 than that the 6.12 
Mev. gamma ray, allowing for the energy variation the pair cross-section. 
The intensity ratio the two peaks was found experiment approxi- 
mately 60—fair agreement, considering that the intensity comparison was 
between different runs, that the amount carbon had probably changed. 
There may some contribution this peak from bremsstrahlung losses and 
wall effects the 6.52 Mev. gamma ray, discussed above connection with 
the 8.27 Mev. peak. 

The relative yields the 6.12, 6.52, 6.78, and 7.29 Mev. gamma rays were 
determined four bombarding energies between 0.8 and 2.2 Mev. The 
excitation curves, shown Fig. within the rather poor accuracy the 
intensity determinations, rise smoothly throughout this energy range, indi- 
cating that the reactions are non-resonant type. 

evidence was found, either during the initial single-crystal survey, 
the later measurements with the three-crystal spectrometer, high energy 
gamma rays from deuteron capture leading excited states 
The excitation energy deuteron energy 0.95 Mev. 26.1 Mev., well 
above all single-crystal background except cosmic rays. Assuming that total 
count three times background corresponds the smallest detectable gamma- 
ray intensity, and assuming efficiency 20% order allow for the severe 
bremsstrahlung losses and wall effects, this indicates that the cross-section 
the reaction less than cm.? 0.95 Mev. bombarding 
energy. 

These results provide additional data the and nuclei, and 
general confirm the mirror character, although, addition those discussed 
above, there are several conflicting details the data. The 4.23 Mev. ground 
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Fic. Yield curves the 6.12, 6.52, 6.78, and 7.29 Mev. gamma rays. 


state transition should have been detectable, since the intensity should 
have been comparable that the 4.46 Mev. transition The fact that 
gamma rays 5.87 and 7.01 Mev. reported Bent (4) were not detected 
the present work may the differences target thickness and 
bombarding energy, but this does not explain why they did not detect the 
5.35 Mev. gamma ray which was quite evident the present work. Measure- 
ments the angular distributions some the gamma rays, for example 
those 6.52 and 6.76 Mev., could aid comparing the mirror levels. Since 
has been found (13) that the neutron groups are strongly peaked forward even 
deuteron energies 1.5 Mev., preferable measure the gamma-ray 
energies 90° rather than 0°, reduce neutron background. The experi- 
ment reported here has demonstrated the usefulness three-crystal spectrom- 
eter accurate measurements gamma-ray energies intense neutron 
flux. 
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THE NEW ISOTOPES AND 


ABSTRACT 


Two new neutron deficient isotopes (59 min.) and (27 min.) have been 
produced the reaction Mev. the McGill University 
synchrocyclotron and the subsequent growth silver from the cadmium. 
These isotopes have been investigated means 180° spectrograph high 
resolution, lens spectrometer, and scintillation spectrometer. Chemical 
methods identification have also been employed. has been found emit 
positron spectrum end point 2.70+.01 Mev. and 556.2 kev. which 
thought series with the positron spectrum. The origin 
energy 118.4 kev. which appears converted palladium present 
unknown. tentative decay scheme proposed. appears decay prima- 
rily K-capture since positron spectrum was observed which could 
assigned this isotope. Four energies 66.7, 83.6, 123.6, and 134.2 kev., 
have been definitely assigned this isotope the basis differences 
shown conversion lines; the basis relative line intensities and K/L 
ratios the two low energy have been shown represent transitions. 
conversion lines were found whose origin and element conversion are 
unknown. 


INTRODUCTION 


Various neutron deficient isotopes cadmium may produced proton 
bombardment silver. The lowest mass cadmium isotope which formed 
depends the energy the bombarding protons. Thus when silver 
bombarded with protons energy less than about Mev. the main re- 
action while above Mev. the reaction 
initiated; the nucleus stable, and the excited product the 
reaction probably too short lived observed the 
techniques here employed. The product the 3”) reaction has been dis- 
cussed previous paper (11). the bombardment energy further in- 
creased above Mev. the reaction begins take place, 
and this paper concerned with investigation the radiations emitted 
the isotope and its daughter product 

The was produced proton bombardment metallic silver McGill 
University’s 100 Mev. synchrocyclotron. All bombardments were carried out 
nominal proton energy Mev. where the activity was produced 
high yield. The two isotopes and were investigated means 
180° permanent magnet spectrograph, lens spectrometer, and scintilla- 
tion spectrometer manner similar that used for the investigation 
(11). Asin the previous investigation all sources were prepared evapora- 
tion the active cadmium from the irradiated silver target. Consequently the 
radiations consisted mixture those emitted and its daughter 
The previous investigation, however, had served identify the con- 


Contribution from the Radiation Laboratory, McGill University, Montreal, Quebec. 

Preliminary results have been reported Paper No. 76, Royal Society Canada, June 
Meeting, Quebec, 1952. 

2Present address: Defence Research Board, Suffield Experimental Station, Ralston, Alberta. 


841 


842 CANADIAN JOURNAL PHYSICS. 


version lines, positron spectra, and y-radiations from and 
and thus was hoped that any new observations could assigned the 
decay chain Analysis some the data, however, indi- 
cated that the isotope might also have been present some extent. The 
fact that the two cadmium isotopes, (55 min.) and (59 min.), had 
similar half-lives did not complicate the investigation greatly. Assignment 
and positron spectra particular mass 104 isotope was made 
where possible the basis K—L differences shown conversion lines, half- 
life considerations, and chemical methods. 


PREVIOUS WORK 


has not been reported previously. 

The only references are very brief and non-specific. Enns (3) found 
two activities half-lives and min. after proton bombardment palla- 
dium. These assigned and respectively. Lindner and Perlman 
(14) reported min. activity obtained amongst the spallation products 
antimony and concluded that there was ‘‘considerable orbital electron capture 
Shore, Brown, and Becker (1) extracted min. silver activity from palla- 
dium which had been bombarded Mev. deuterons and found internal 
conversion lines energies 530, 550, and 740 kev. They assigned these lines 
the and lines kev. y-ray to7) and the line 
kev. y-ray respectively, both converted palladium. definite 
assignment was made the mass number the silver isotope involved, but 
masses 104 and 102 were cited possibilities. Recently this activity has been 
assigned Haldar and Wiig (7) who prepared activity identical 
half-life bombardment silver with high energy protons. These authors 
also showed that this isotope decayed emission positrons 1.3 Mev. end 
point. 

MEASUREMENTS WITH THE 180° SPECTROGRAPH 


The values and energies the conversion lines detected the 180° 
spectrograph are listed Table together with their most probable assign- 
ments. Approximate relative intensities the most prominent lines are also 
included, these being derived from K/Z ratios measured the lens spectrom- 
eter and from visual estimates the intensities weak lines observed the 
photographic plates. The added binding energies were obtained from the most 
recent values X-ray level energies (9). was pointed out above, the radi- 
ations the higher mass cadmium and silver isotopes were also present 
sources prepared Mev. bombardment energy; only the lines which were 
found after bombardment the target Mev. and which were not found 
Mev., see Reference (11), are listed Table Thus the table should 
contain only those lines emitted the decay number 
the unassigned weak lines might due will noted that num- 
ber the lines are assigned being converted definite element (as 
Pd); this has been done only where both and conversion lines were ob- 
served. Two doubtful lines are indicated such. 
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TABLE 
INTERNAL CONVERSION LINES anp DETECTED 180° SPECTROGRAPH 


Approximate 


Hp, relative 
kev. Assignment intensity kev. 

834.9 58.01 500 
854.7 60.65 

991.9 80.23 83.61 
1009.2 82.85 83.59 
1080.5 94.02+0.15 118.38 
1169.4 108.70 134.24 
1204.8 114.77 118.38 
1254.3 123.45 

1260.8 124.61 

1294.3 130.63 134.18 
1421.4 154.37 

1441.9 

1783.0 
533.8 
3111.1 552.5 556.1 
3120.7 

3138.6 559.7 

3271.6 594.9+0.6 

600.5 

3294.1 (?) 

8.8 604.8+0.7 

629.1 
631.2 
3514.2 
3532.1 (?) 664.5 
669.1+0.8 

684.6 
733.7 

3821.9 


Figure reproduction spectrograph plate taken the region below 
1100, and shows the most prominent conversion lines which were observed, 
those due the y-rays energies 66.7 and 83.6 kev. emitted 
(converted silver). For the 66.7 kev. y-ray the only line which was ob- 
served was the case the 83.6 kev. y-ray all three lines were visible 
the original plate, with being very much the strongest; visual estimate 
the relative intensities indicated 

the liries above 1100 the most prominent were those the 556.2 kev. 
emitted the decay For this y-ray the ratio was visually 
estimated 7+3. The energy value obtained for this y-ray accord with 
values 550+10 (12) and 560+10 (15) found for y-ray investigations 
the decay these latter investigations also serve confirm the present 
assignment mass 104. 

the weak lines the ones most readily observed were those energies 94.02 
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and 98.09 kev., followed those energies 108.70 and 123.45 kev. The first 
three these were assigned lines, and the assumed lines related them 
were weak that reliable visual estimate the K/L ratio could not 
made. would appear that these cases K/L All other lines were very 
weak and many could observed only after repeated exposures many 
sources the same plate. 


MEASUREMENTS WITH THE LENS SPECTROMETER 
4.1 Conversion Lines 

The decay the conversion lines the 66.7 and 83.6 kev. y-rays was 
measured the lens spectrometer and values and 58+2 min. respec- 
tively were obtained for the half-lives. The decay rate the line the 66.7 
kev. was subject greater uncertainty than that the line the 
83.6 kev. y-ray since the counting rate the peak was just over twice the 
rather high background (positron) count this region. both cases the 
and conversion lines were weak that reliable decay rate determination 
over sufficiently long period time could not made. 

The K/L ratios for these two y-rays were obtained planimetric measure- 
ment the areas under the lines; the case the line the 83.6 kev. 
was necessary subtract the small contribution due the partially 
resolved and lines. For the 66.7 kev. y-ray this ratio was 10+3, and for the 
83.6 kev. y-ray The ratio was ~13. These ratios will 
referred later aid the assignment the multiple orders the y-rays. 

the other conversion lines found means the 180° spectrograph the 
only one which could detected above the intense positron continuum the 
lens spectrometer was the conversion line the 556.2 kev. y-ray. Since 
previous measurements the 180° spectrograph had indicated that this y-ray 
was converted palladium, determination its decay characteristics could 
give valuable information which would aid the assignment half-lives and 
positron spectra. After subtraction the positron background under the 
conversion line the data were sufficient indicate 27+4 min. growth and 
then 59+4 min. decay. This evidence must considered along with that 
provided the positron spectrum (see below). 


4.2 Positron Spectrum 

The positron continuum was investigated and found complex. This was 
expected owing the presence well the mass 104 isotopes. 
The first run the spectrometer showed the already known 1.69 Mev. spec- 
trum and also new spectrum extending 2.7 Mev. This latter 
spectrum showed growth and another run number points this con- 
tinuum were followed for eight hours establish the decay characteristics. 
Analysis showed that this spectrum eventually went into simple 59+2 min. 
period, and subtraction the actual counting rate from the extrapolated 
portion the final decay curve led value 27+1 min. for the rate 
growth. Such growth and decay curve shown Fig. will seen that 
the two straight lines, extrapolated backwards, meet point within min. 
the time which the cadmium was evaporated (chosen zero time), 
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— POINTS 


2721 MIN 1, 2592 2 MIN 


; Ag 104 Cd 104 


s 


COUNTS PER SECONDS 


120 180 240 300 360 420 480 
TIME MINUTES 


Fic. and decay the positron spectrum emitted observed the lens 
spectrometer (at 1.78 Mev.). The time which the cadmium was evaporated was chosen 
zero time. The rate decay eventually controlled the decay the parent, min. 
The standard deviation shown. 


indicating that one the periods was not present initially. Since conversion 
lines had already been found half-life ~59 min., and these showed growth, 
was evident that the min. activity was daughter product the 
min. activity and should therefore assigned This assign- 
ment partially confirmed the observed characteristics the 556.2 kev. 
y-ray; this y-ray was converted palladium and appeared show the same 
growth and decay periods the positron spectrum which being considered. 

With knowledge the decay characteristics the positron spectrum, 
further runs were made provide data for Fermi analysis and accurate 
determination the end point. Owing the relatively short half-lives in- 
volved and the complex nature the decay, and considering the time 
necessary sweep through the whole spectrum, the complete spectrum was 
investigated two sections with fresh source for each run. the first run 
many points were counted the region 1.5 3.0 Mev. and three below 1.5 
Mev. the second run points the region below 1.8 Mev. and three points 
above 1.8 Mev. were counted. All points were followed for period eight 
hours and growth and decay curves drawn for each. Counting rates corre- 
sponding times were read from these curves the usual manner. complete 
Fermi plot the whole spectrum was made combining the two sections and 
normalizing the common points. Values the Fermi function were taken 
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from the National Bureau Standards publication AMS 13, and the screening 
correction was incorporated. Such plot shown Fig. The main feature 
this plot the high energy spectrum with end point 2.70+.01 Mev. which 
emitted The spectrum end point 1.69 Mev. due the 


150 


FERM PLOT 


5S 
ENERGY MEV 


Fic. Fermi plot positron spectrum emitted evaporated cadmium source prepared 
Mev. The standard deviation for the high energy spectrum larger than the plotted 
points. 


upward curvature beginning about 0.8 Mev. similar that found the 
investigation and cannot considered new spectrum. This 
upward curvature was originally (10) interpreted positron spectrum 
0.93 Mev. end point and was assigned because did not show the 
growth and decay characteristics the spectrum. Further work, how- 
ever, established that the general form the curvature was similar that 
found the spectrum (11), and indicated that this spectrum low 
that there evidence for spectrum end point 1.3 Mev. Thus the decay 
chain only one positron spectrum evident, that 


MEASUREMENTS WITH THE SCINTILLATION SPECTROMETER 


Further investigation the gamma radiations was carried out means 
scintillation spectrometer which had been built Martin and 
Breckon (16) and Henrikson (8). single channel analyzer was used obtain 
pulse height distribution, the source being covered with Mev. aluminum 
absorber absorb all positrons. 

the high energy region new peaks were found which had not already 
been established belong The region below 1.4 Mev. shown 
Fig. The most prominent peaks are those 1.30 Mev., 930 kev., and 270 
kev. already assigned that 740 kev. which probably due 
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930 KEV 


1.30 MEV 


PULSE HEIGHT VOLTS 


Fic. Scintillation spectrometer pulse height distribution the y-rays from 
and (source produced Mev.). Successive points have not been corrected for half-life. 


mixture the y-rays associated with the internal conversion lines energies 
684.6, 733.7, and 743.6 kev. listed Table the annihilation radiation 511 
kev.; and the large peak due the 93.0 kev. y-ray addition the 
effect the 556.2 kev. y-ray and possibly the y-ray (about 559 kev.) 
associated with the internal conversion line 533.8 kev. may seen hump 
the high energy side the annihilation radiation peak. Other slight humps 
the curve indicate the presence other y-rays, but with the profusion 
this region further analysis was not warranted. The peak assigned 
the 93.0 kev. y-ray appears show structure, but with the poor resolution 
this region nothing definite can said about the presence peak due the 
66.7 83.6 kev. y-rays. The hump the high energy side this peak 
unexplained but was also found later run after bombardment 
Mev. and might belong might scattering effect. 


CHEMICAL INVESTIGATIONS 


order verify that the min. activity assigned did fact belong 
silver, attempt was made isolate chemically. This was rather im- 
portant since the gross decay evaporated cadmium sample, produced 
Mev., showed signs such activity, the decay curve being resolvable 


one 
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only into 6.7 hr. and 57-60 min. components. Consideration the laws 
radioactive decay shows that daughter activity maximum time 
tmax after isolation from its parent given 


=> a 


where half-life parent, half-life daughter. For the case 
same obtained for the case min. and min. Evidence has 
already been presented the effect that the former the case. 

The chemical procedure was essentially that dissolving the evaporated 
cadmium dilute nitric acid with added cadmium and silver carrier, precipi- 
tation silver means hydrochloric acid remove any active silver 
which might have been evaporated, and filtration. The filtrate was then allowed 
stand for about min. order allow grow from this being 
about the time necessary for the reach its maximum activity. little 
silver nitrate was then added and the precipitate was collected filter 
stick, redissolved ammonium hydroxide, and then precipitated again 
acidification. portion the precipitate was dissolved ammonium hy- 
droxide, and drop two pipetted onto glass slide. This was then dried 
and counted. 


10° 


@ EXPERIMENTAL POINTS (LESS BACKGROUND) 
© AFTER SUBTRACTION OF 6.7 HR COMPONENT 
© AFTER SUBTRACTION OF S9 MIN COMPONENT 


SECONDS 
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6,7 HR 


TIME MINUTES 


Fic. Gross decay chemically separated measured with end-window Geiger 
Counter. The standard deviations are shown for the original points and for the min. com- 
ponent 
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The gross decay the silver chloride precipitate, measured with end- 
window Geiger counter, shown Fig. The initial min. half-life 
very good agreement with that obtained the lens spectrometer. Although 
attempts were made purify the silver chloride precipitate from the very 
large excess cadmium activity, there was still noticeable amount (usually 
about 20% computed the time final separation) min. activity 
present, well still longer activity. Since this 55-65 min. activity was 
present such large amount, must assumed that was due mainly 
silver isotope (which must grow from the evaporated cadmium). yield 
arguments might expected that Mev. the production 
will appreciable, perhaps greater than that and possible that 
this 55-65 min. activity due the assumption that the primary 
yield could twice that this energy, the presence this 
amount indicates that might have lifetime the order 
min. Part this min. activity may also accounted for the 
known adsorption property silver chloride precipitate for metal ions, 
pointed out Bradt (2), but this will minor effect. The longest 
period, will seen graph, appears approach half-life the 
order 6-8 hr. (suggestive the presence Accordingly, line 
6.7 hr. half-life was drawn the assumption that the experimental decay curve 
approached such line asymptotically (the fact that the 6.7 hr. line goes 
through the last experimental point accidental, this point appears low 
compared previous points). While the contribution the 6.7 hr. activity 
(assumed very small, serves indicate that some the cad- 
mium activity may have been adsorbed the silver chloride, else that the 
separation was not complete. Subtraction this 6.7 hr. component from the 
experimental points yielded line half-life min. for the second com- 
ponent the decay. the basis half-life this component could due 
pointed out above, with small admixture adsorbed (55 
min.) and (59 min.). 

connection with the suspected presence should pointed out 
that present, would not affect the determination the end-point 
the positron spectrum, nor the half-life measured the lens 
spectrometer, since both these determinations depend data taken beyond 
the end point the positron spectrum 1.3 Mev. (7). The absence 
evidence for the positron continuum the Fermi plot Fig. does 
not give conclusive proof that this isotope was not present, since expected 
from ratio considerations that the number positrons will 
small compared those the other isotopes present. 

Since and both have isomeric states with half-lives the order 
sec. and emit y-rays 93.0 and 87.9 kev. respectively, the discovery 
y-ray comparable energy, 83.6 kev. (and also converted silver), raised 
the question whether there might similar short-lived state measur- 
able period here also. The fact that was even mass made this hypothesis 
unlikely. Nevertheless, search was made using scintillation spectrometer 
and one channel fast counting apparatus (16) available the laboratory. 
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The procedure followed was precipitate out all the silver from nitric 
acid solution the evaporated cadmium manner similar that described 
previously, and then quick chemical separation the pre- 
cipitate after had been allowed grow for varying lengths time (of 
the order minute). The elapsed time from the separation the 
until counting was started was less than minute. 

Since the y-ray had energy only 83.6 kev. the lower bias level had 
set just above noise level, and was found absorber could not 
used over the source. Thus the counting rates obtained were very high the 
rather wide channel that was used. The first run gave initial decay sec., 
but subsequent runs gave values ranging from sec. which seemed 
proportional the time the was allowed grow before being 
precipitated. all these cases the counting rate decreased not more than 
factor two the counting period. Finally, about seven hours after the cyclo- 
tron bombardment, when most the min. and min. had 
died out, the separation procedure was repeated order see whether the 
44.3 sec. could detected, and consequently the method was valid. 
This run was very successful and gave period sec., with the counting 
rate decreasing more than factor the counting period. 

These results not appear indicate definitely the presence any 
activity half-life sec. which could assigned metastable state 
fact felt that the periods 21, sec., etc., which 
were obtained, could have been due being counted top very 
high background, the latter causing pile-up pulses, could have been due 
deterioration the photomultiplier dynode voltages caused the neces- 
sarily high counting rates. Also the uncertainty the measured period might 
have been caused form temporary fatigue the photomultiplier itself, 
such fatigue being known appear high counting rates and disappear 
again lower rates. This problem has not been investigated further, but 
there any short-lived activity present, either would appear 


DISCUSSION RESULTS 


7.1 Mass Assignment 


The mass assignment the min. rests considerations the 
systematics nuclei this region the periodic table. 
nucleus only two mass units removed from the stable For this reason 
expected that this isotope would have readily measurable period. 
the other hand, continuation the trend (1.3 yr.), 
(6.7 hr.), and (55 min.), would expected have shorter half-life 
the order 2-20 min. might expected short-lived. 

Rather more substantial arguments may obtained consideration the 
threshold energies for the (p, and (p, reactions 
These may calculated (17) Mev., Mev., and Mev. respectively. 
Experimental results indicated that the activity which was assigned 
began appear nominal proton bombardment energy just above 
Mev., which would favor assignment 


a _ 
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Since the min. was observed grow from its mass assign- 
ment depends that its parent, and very little additional discussion 
necessary. One the isomers has very similar half-life, 24.5 min., 
but the fact that its positron end point 2.04 Mev. very different from the 
2.70 Mev. observed for the min. activity assigned removes all 
doubt that the two could have been confused. Also the observed growth, the 
method production, and the energetics preclude assignment mass 106. 
The fact that the 556.2 kev. y-ray assigned probably identical with 
y-ray found the decay (12, 15) also lends support the proposed 
mass assignment. 


7.2 Multiple Order y-Rays 

The relative intensities and internal conversion lines, 
well intensity ratios, provide information concerning the multipole 
order the y-ray with which the lines are associated. 

obtain information about the multipole orders the 66.7 kev. and 83.6 
kev. y-rays use made the experimental fact that the conversion lines 
was very much the strongest both cases. Comparison with theoretical 
shell ratios (5) indicates that both these transitions are this assignment 
further confirmed for the 83.6 kev. y-ray the fact that the line was 
found more intense than the line. This assignment also accord 
with graphs presented Tralli and Lowen (24) concerning the 
Further information provided the measured ratios, for the 66.7 
kev. y-ray and 7.5 for the 83.6 kev. y-ray. Comparison with empirical K/L 
ratios (6) would appear confirm that both the y-rays represent 
transitions. detailed examination theoretical K/L ratios for and 
transitions has been made Swann and Hill (23) who plotted the ratio 
conversion coefficients (20) and the total conversion coefficients 
Gellman, Griffith, and Stanley (5). These authors point out that the theo- 
retical K/L ratios show significant dependence, and examination their 
plotted curves shows that for the values (33.2 and 26.5) applicable 
the two y-rays question K/L ratio might expected for 
transition the region (for Ag, 47). Such K/L ratio within 
the limits the experimentally determined ratio for the two y-rays. Higher 
multipolarities may excluded the basis both ratio and lifetime. 

For the case the 556.2 kev. y-ray the value 3.8, and this 
range difficult assign the type transition without ambiguity solely 
the basis the K/Z ratio (estimated ~7). From the curves given 
Swann and Hill will seen that the region small for the 
K/L ratio curves for and transitions appear approach each other and 
indicate K/L ratio ~9. Indeed the observed K/L ratio cannot used 
with certainty eliminate any the possible multipolarities (6). The fact 
that the conversion lines only was observed does not provide choice 
between any one M1, transitions this energy (5). Consideration 
the observed intensity the conversion line this y-ray relative the 
2.70 Mev. positron spectrum and the conversion coefficient expected 
this energy (21) indicated that the number y-rays amounted least 
50% (if E2) the number positrons; hence likely that this y-ray 
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series with the 2.70 Mev. positron spectrum and leads the ground state 
Thus accord with the evidence (6, 22) that the first excited states 
nuclei can expected have spin and even parity reason- 
able assign the 556.2 kev. y-ray transition. This assignment also 
accord with the observation (18) that, for low and high energy, conversion 
the shell becomes increasingly important for transitions. 

For the case the three remaining y-rays for which both and conversion 
lines were observed the very low intensity the conversion lines makes 
multipolarity assignment difficult. For the 118.4 kev. y-ray (converted 
palladium) assignments are most probable (5). For the 123.6 and 
134.2 kev. y-rays (converted silver) transition perhaps likely 
(5, 18). 


7.3 Decay Schemes 


Conclusions decay schemes may drawn from considerations 
positron spectra, y-ray intensities, and theoretical predictions. the decay 
chain convenient start from the stable end product 
Pd!°, 

Since nucleus its ground state might expected 
have spin and even parity. The theoretical (4) K-capture positron ratio 
for the decay ~0.18, calculated the basis the 2.70 Mev. posi- 
tron spectrum. Thus for positron emission the predominant process 
(except perhaps very much higher levels) and y-rays accompanying 
capture other levels might expected difficult detect the tech- 
niques employed. was mentioned above, since the 556.2 kev. y-ray the 
decay far the most intense y-ray which could definitely 
assigned being converted palladium, logical conclude that 
follows positron emission and leads the ground state The position 
the first excited state, 2+, 556.2 kev. compares closely with the 
corresponding state 513 kev. and consistent with the systematics 
these states this region the periodic table decay scheme based 
these ideas shown Fig. 

The positron end point 2.70 Mev. taken conjunction with the experi- 
mentally determined half-life min. gives log product 5.4. The 
value used here was determined from the graphs given Feenberg and Trigg 
(4) and includes the theoretical correction for K-capture. the absence 
further information regarding other modes decay this figure may taken 
minimum value for log ft. The apparent absence other positron spectra 
intense y-rays which could assigned this decay would indicate that the 
log value 5.4 might expected near the correct value. This implies 
that the transition allowed; assuming the correctness the rule for the 
first excited states nuclei, the ground state the nucleus 
may expected have spin and even parity. The absence 
any detectable positron emission the ground state however, 
indicates that the spin the basis the shell model the nucleus 
(13), the configurations 7/2, and are indicated for the 


| 
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Fic. Proposed decay scheme for 


odd proton and odd neutron orbitals respectively The parity assign- 
ment even excludes the configuration. The remaining configurations 
can yield spin values although simple application the empirical 
rules would favor spin (19). 

The decay scheme for presented Fig. accord with that 
proposed for the decay sec. (12, 16). The y-ray energy 118.4 
kev., which was assigned palladium the basis two conversion lines whose 
separation was equal the K—L, separation palladium, does not fit into the 
proposed scheme. may that this y-ray represents transition between 
two states higher energy the difference may fortuitous 
and the two conversion lines may actually represent lines y-rays converted 
silver whose origin and position are present unknown. The second 
excited state may reasonably expected lie about Mev. 
(12, 15, 22), and possible that one the high energy conversion lines 
whose assignment was not obvious might represent transition from the 
second the first excited state y-ray which could have represented 
transition from this second excited state the ground state was not 
observed. Since the position such second excited state was uncertain 
has not been included the proposed decay scheme. 

The apparent absence positron transition the decay makes 
the assignment spins and parities difficult any decay scheme which 
proposed, and thus provides corroborating evidence for the state 
assigned above the ground state approximate argument may 
used show that positron transition the decay may indeed 
expected weak and that the y-rays which were observed probably 
arise from K-capture. Using the energies the decay and the 


2mc 
2.70 
(24) 
ow 
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decay, together with the experimental correction term the 
semiempirical mass formula found for nearby masses, possible construct 
approximate mass parabolas and estimate that the energy available for the 
decay about 0.7-0.8 Mev. This value could error almost 
factor two, but nevertheless shows that positron emission, not actually 
forbidden energetically, will negligible compared K-capture. 

While there are numerous suggested relationships between the energies 
the observed conversion lines, the complexity the y-ray spectrum and also 
the possibility that might have been present contribute the observed 
radiations make unwise attempt construct decay scheme for 
will noted that the energy differences the pairs conversion lines 600.5 
and 533.8 kev., and 626.7 and 559.7 kev., are very close the energy the 
66.7 kev. y-ray, and the differences the lines 743.6 and 659.9 kev., and 684.6 
and 600.5 kev., are very close the energy the 83.6 kev. y-ray. This ap- 
parent relationship might indicate that these lines are converted silver and 
belong the decay. also tempting postulate that the 66.7 and 
83.6 kev. y-rays are series, appears possible virtue the presence 
the weak conversion line energy 124.61 kev. which, assumed the 
line y-ray converted silver, would give the exact y-ray energy 
expected for crossover transition. Many other relationships may obtained 
between the observed conversion lines between y-ray energies obtained from 
the conversion lines the assumption that conversion takes place silver, 
but further data required corroborate any decay scheme which 
might proposed. was stated above, some the conversion lines which 
were observed might due There also the possibility that some 
them may represent weak transitions the decay which were not 
observed the previous investigation. 

order verify and extend the above decay scheme for and 
obtain more evidence for scheme for would desirable perform 
coincidence experiments, but the very large number y-rays the energy range 
500-770 kev. would necessitate the use rather high resolution. Nevertheless 
valuable information could obtained even with moderate resolution, par- 
ticularly coincidences were sought between y-rays and the adequately 
intense conversion lines the 66.7, 83.6, and possibly the 556.2 kev. y-rays. 
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STANDARD MUTUAL INDUCTANCE! 


ABSTRACT 


primary standard mutual inductance the Campbell type has been 
constructed the National Research Council serve ultimate standard 
reference for electrical quantities Canada. The electrical value computed 
from its measured physical dimensions. detailed description warranted 
view its national importance. The value the standard henry 
with error not larger than p.p.m. 


DESCRIPTION 


Two the first considerations making any standard are stability and 
accuracy. This particularly true the objective produce primary 
standard. Consequently the earlier experience the National Physical Labora- 
tory and the United States National Bureau Standards has been followed 
and both laboratories have been most helpful and co-operative throughout 
the program. Modifications were made the original techniques when 
seemed probable that the final result would improved. 

The Campbell mutual inductance (3, represented Figs. and The 
primary has two equal single layer coils 100 turns bare copper wire wound 
helical groove quartz cylinder. The cylinder roughly cm. diam- 
eter, and the windings, having mm. nominal pitch, occupy axial length 
cm. each and are separated distance cm. The secondary has 
459 turns insulated wire wound groove marble ring the mean 
diameter which approximately cm. The secondary located midway 
between the two primary windings and coaxial with them. 

well known the main advantage design that the mean 
diameter the secondary coincides with the neutral circle (the circle zero 
magnetic field between the primary coils when the latter are excited series). 
This means that small errors physical dimensions the secondary have 
relatively small effect the calculated inductance compared the effect 
similar dimensional errors the primary. 

marble base, together with marble supports which hold the secondary 
ring and provide horizontal and vertical adjustment the latter, completes 
the instrument. The ends the electrical leads from the primary coils must 
one straight line within tolerance +0.1 mm. The inductor housed 
laboratory having controlled temperature and humidity, and built 
without use magnetic materials and avoiding closed electrical loops large 
metallic areas. 

Our constructional methods were very similar those described the U.S. 
National Bureau Standards (6, 7). The primary quartz cylinder, obtained 
before the war from ‘‘Quartz (Paris), was lapped true cylinder. 
The thread hold the primary winding which extends the entire length the 
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Fic. Primary inductor. 
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Fic. Principal dimensions inductor. 


cylinder was also lapped its final shape starting from initial helical saw 
cut 1/100 in. square. The process was stopped when photographic profiles 
showed that the primary wire would supported along two generating lines 
without tendency roll touch the bottom the groove. 

Measurements made preliminary coil wound the full length the 
cylinder showed there was optimum position for the final coils, the cylinder 
being remarkably uniform. The final winding the two 100-turn coils was 
then made, using bare copper wire which was drawn through series dies. 
The last diamond die was checked for quality with Taylor Hobson projector 
and shown circular. Samples the wire taken from the beginning and 
end each primary helix were measured. addition the relative diameter 
the wire was continuously checked during the winding with Zeiss microm- 
eter. The uniformity the readings showed there was significant differ- 
ence between the diameters the samples and the helices. (The wire the 
helices carried one half turn beyond the electrical end the coils simplify 
securing the end the wire.) 

The base, the supports for the secondary, and the secondary ring itself were 
all made the National Research Council Central Workshops. The measure- 
ments the preliminary coil mentioned above were accurate enough 
compute the position the neutral circle the precision required cut the 
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groove the secondary ring. The mean diameter the groove should not differ 
from that the neutral circle more than 0.01 mm., but high grade machine 
shop practice was quite adequate achieve this. 

constructed six separate secondary rings used with the same pri- 
mary. Three these secondaries will yield mutual inductors mh. nominal 
value, two 15.8 mh., and one mh. All them can serve computed 
standards but present the best secondary winding giving mh. taken 
the standard. Five secondary forms are constructed Alabama marble and 
one form Italian marble. The extra secondaries are intended for other work 
and will not discussed detail here. All were wound using the same tech- 
nique. Each layer the winding separated from the next insulating fish 
paper 0.005 in. thick and the resulting winding very uniform, important 
point because only the outer wires each layer can seen during measure- 
ment the diameter the coil. Electrically the secondary coils terminate 
coaxial connector. 

Although were using established techniques, the laboratory had make 
all the necessary laps and special micrometers, and modify many existing 
gauges carry out the work precisely. 


MEASUREMENT 


All measurements the primary coils are required the highest precision 
attainable. The limiting factor the accuracy probably the fact that are 
attempting define the position the primary current measurements made 
hard drawn copper wire, relatively soft material, and any significant 
improvement can come only from complete change approach. 


2.1. Diameter the Primary Coil 

substitution method developed Moon (12, 13) was adopted. fixed 
and moving anvil micrometer are secured diametrically opposite one 
another large invar ring. The moving anvil motor driven, and auto- 
matically stopped when predetermined measuring pressure (150 gm.) 
reached. large dial attached the rotating anvil reads directly microns 
and measurements can repeated with ease and certainty. (See Figs. and 4.) 

The reference standard quartz end gauge, having optically polished 
spherical ends, the center the sphere being the center the gauge. The 
gauge placed shaped support suspended above the primary cylin- 
der. When the invar ring elevated pick the the quartz end gauge 
automatically placed the same position each time with respect the 
micrometer anvils. 

The primary quartz cylinder mounted with its axis nearly vertical immedi- 
ately below the quartz gauge. The axis tilted compensate for the helical 
angle and this must allowed for diameter measurements, but the method 
has the advantage keeping the micrometer axis horizontal. Diameters were 
measured each wire along four pairs generating lines 45° intervals 
around the cylinder. 

The measurements were made groups turns each, the time interval 
being divided into halves, symmetrical with respect the mid-point. Every 
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UPPER HELIX 


100 


TURNS 


Fic. individual turns, compared average radii (upper helix) and (lower 
helix). the mean that used for inductance computation. 


wire was touched three times the first half and three times the second 
half the group. The cylinder was rotated slightly after each measurement 
that the anvils did not touch the wire twice the same spot. The gauge was 
measured three times each half group. The quartz cylinder was turned end 
for end after the first helix was completed, facilitate measurement the 
other coil. 

minimize thermal effects the whole equipment was enclosed box and 
two thermometers read before and after each set measurements. Normal 
agreement between the two was few tenths degree. The change tem- 
perature for whole day operation did not exceed one one and half 
degrees. The thermal coefficient the quartz gauge was found experimentally 


TURNS 
(b) LOWER HELIX 
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Fic. Ring comparator, showing quartz gauge its Y-shaped support, which has point 
line-plane contact with invar ring. 

Fic. Longitudinal comparator. The only portion rider assembly visible mirror 
(near right-hand microscope objective). 
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good agreement with accepted values found the tables (0.5 0.7 
p.p.m. per degree C.). Assuming that the cylinder has the same coefficient, all 
errors due difference temperature between the gauge and the cylinder 
should have averaged out. The remaining error estimated not exceeding 

measure the compression the copper wire pressure 150 gm. 
exerted the micrometer anvil short pieces the same wire were fixed 
auxiliary gauge and from measurements made the Interferometry Section 
the compression was found the order 0.1 0.2 The value 0.15 
was adopted the most probable and employed correction. 

The errors arising from the ring comparator itself, insofar they are 
mechanical origin, are the order 0.2 0.3 The sum all these errors, 
i.e. all those due the comparator, amounts about 0.3 0.5 

The mean diameter the primary winding (which used the inductance 
calculations) the average several thousand measurements that the effect 
accidental errors the final result may considered very small. 


2.2. Length Quartz Gauge 

The length the quartz gauge used above was measured interferometrically 
for Baird and given 0.302,6087 The error the length the 
gauge may estimated between 0.2 and 0.4 Measurements this gauge 
about year earlier tend confirm this, the difference between the two 
determinations was only 0.1 

the same time estimate was made the error that might appear the 
diameter measurement owing the quartz gauge being misaligned the ring 
comparator, but was shown this would not exceed 0.1 any misalignment 
that could escape detection. 


Axial Lengths Primary Helices 

The Abbe longitudinal comparator method was found most convenient for 
axial measurement the primary helices (Fig. 5). According this method 
scale and generating line the cylinder were placed the same straight 
line and observed with two microscopes rigidly connected one another. 
One pointed the scale while the other pointed the turns the winding. was 
preferred move the microscopes together rigid mounting rather than the 
heavy cylinder and scale because obvious difficulty preserving their 
relative positions. The distances measured were the axial length each 
winding and the distance between them. 

Some useful preliminary data were obtained pointing the microscope 
the lines light reflected the wires from the axial illumination the 
microscopes. This method can relied upon only few microns however. 

the precise measurements the definition the position the wire was 
obtained from mechanical brought into contact with the copper 
wire analogous the system described Volet (11) for similar work 
self-inductance standard. fine line (Fig. defines the center the 
wire. the order 1.5 mm. and the should vertical, very 
sensitive mechanism necessary level the rider (or least ensure that 
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Quartz 


Vert. adj. 


Hor. adj. 


(b) 


Fic. (a) Cross section rider. 
Rider adjusting device. light tube fastened rider rests bearings 


and Rider and mirror can raised thread over pulley 


always the same relative position the wire). The levelling mechanism 
shown Fig. sensitive level bubble adequate indicator correct 
positioning the rider but for the interferometric observations small 
mirror was rigidly fastened the rider which automatically assumed the 
proper position when the mirror was set correctly. 

The microscopes (Fig. are joined two invar rods. These rods are secured 
one end only rigid aluminum casting sliding the ways lathe bed. 
Checked autocollimation the ways are straight, horizontally and vertically, 
one second arc. 

The distances between turns were determined terms the intervals 
52% length standard, which was calibrated against invar meter, which 
checked periodically against the N.R.C. standard meter, the latter being known 
directly terms the international prototype. 


2.4. Interferometric Measurement Axial Dimensions 

recent improvement Baird (1) interferometric technique, materially 
shortening the observation time required for optical determination 
length, made practical apply these axial measurements. Accordingly 
series measurements the displacement the rider were made 
Baird and one wherein the displacement was obtained optically. These 
measurements were simultaneous with those described above. The rider, with 
its associated errors, common both methods. 


K 
Rider 
Horizontal projection 
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Baird’s system applied this case shown schematically Fig. The 
collimated light from 198 source is, after has passed through beam 
splitter reflected from fixed reference mirror and from mirror 
rigidly secured the rider, producing interference fringes the telescope 


Fic. Schematic diagram interferometer. 


The air pressure cell the path SmS controlled the piston 
essential feature his method the use auxiliary oscillating piston 
that applying the principle flicker photometry the optical paths from 
the two mirrors can adjusted make the order interference exactly 
integral. 

The difference setting piston when the mirror displaced 
indication the fraction the displacement exceeds whole number half 
waves. Using three colors (the lines 5461, 5777, 5791) and applying the 
method exact fractions (5) the displacement can accurately and quickly 
determined. (The required preliminary knowledge the measured distance 
within was easily supplied the comparator 

this measurement the beam light first made accurately parallel 
the direction travel the microscopes, then the rider set wire 
midway between the wires measured. The mirror set perpendicular 
the light beam and then adjusted give sharp fringes nearly zero 
order. Now when the rider moved the wires whose spacing required and 
adjusted that mirror produces good fringes, this automatically ensures 
that the rider always the same relative position with respect the wires 
and that the engraved line will focus the microscope. Separate settings 
were used measure the two coils and the spacing between the coils 
minimize the broadening fringes due long path differences. The greatest 
path difference this work was cm. 

ensure stable conditions for the optical measurements the comparator 
was placed separate enclosure inside temperature-controlled room. 
air pump evacuating the hot air produced the mercury lamp causes slow 
exchange air between the room and the enclosure. The optical measurements 


ch 

q 
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were made first observer outside the enclosure. Then, immediately 
following, the wires were pointed with microscopes. the later observations 
lasted only few minutes the presence the observer inside the enclosure 
did not affect the temperature the comparator the enclosure. mercury 
thermometer placed close the scale was read after the pointings and all 
results reduced 20°C. Temperature changes did not exceed rate 0.01°C. 
min. and reasonable expect that the reduction 20°C. gave rise 
only very small errors which were partly eliminated repetition readings. 

The sources error the optical measurement include: obliquity the 
measured direction with respect the axis the helix, tilt mirror with 
respect the required direction multiplied the distance the rider from the 
center obliquity the interfering beams from finite aperture, error 
piston setting, and corrections due the atmosphere. All these are very 
small easily corrected for and unlikely that their total contribution 
exceeds 0.1 0.2 

more serious source error due change conditions between settings 
the rider. due both changes the dimensions the comparator and 
interferometer and the change pressure the gas the air cell. The 
latter resulted from differential changes different parts the air system 
caused insufficient temperature control the pistons. These effects should 
produce random errors and the scatter readings indicates their magnitude. 
The scatter was the order 


2.5. Measurement the Secondary 

The diameters the secondaries were measured laboratory comparator 
against invar meter bar temperature-controlled room. The secondary 
windings were viewed through radial slots cut 90° intervals the marble 
form. Pointings were made directly the wire. The radial depth deduced 
directly from pointings made the inside and outside layer when determining 
the coil diameter. The axial length (or width) the coil was measured 
with good hand micrometer. Adequate accuracy for all these dimensions 
can obtained this way. 

RESULTS 

The limit precision the Campbell inductor set the errors assigned 

the measurements the primary helices, the mean radius and the axial 


distances 

The mean radius deduced from ring gauge measurements both helices, 
corrected for the slight tilt during measurement and for the compression 
the wire, is: 


0.150,984,28 


When describing the diameter measurement the following errors were 
included: 


Error comparison diameter against reference 
gauge 0.2-0.3 


Error due temperature 


— 


866 CANADIAN JOURNAL PHYSICS. VOL. 


Error length gauge 0.2-0.4 
Error wire diameter 0.1-0.2 
Error due misalignment gauge 0.0-0.1 
Sum the errors 0.6-1.2 


reasonable conclusion seems that even some small cause error 
was overlooked, the upper limit the total error the diameter about 
micron and therefore the radius, half micron.. 

The axial measurements, (length lower helix), (separation be- 
tween coils), (length upper helix), are used determine the distances 
beginning and end the coils from the neutral plane, viz. }BC and 

Because the extra half turn the wire beyond the electrical end the 
winding, these distances could measured along three lines displaced 90° 
from one another around the cylinder. The lengths the two coils agree within 
0.6 0.7 The individual scale measurements and the interferometric mea- 
surements agree about 0.5 Some the discrepancy may due obser- 
vational technique and some slight imperfections the coil. The systematic 
part the error measuring the displacement the rider about 0.2 and the 
scatter individual settings about the same, one left with the over-all 
impression that the rider technique used here probably defines the position 
the wire about half micron. 

The measurements are: 


Scale Interferometer Average 


that 
m., 


Considering the maximum deviation measurement, one arrives 
upper limit for the probable errors of: 


The pitch was taken uniform. Actually every 10th turn both helices 
was measured one line and significant departure from the ideal case was 


observed. 
CALCULATION MUTUAL INDUCTANCE 


The value the inductance computed from the measured physical 
dimensions and assuming value permeability vacuo. Jones’ formula 
(10) used and three corrections were applied allow for (a) the finite cross- 
section the secondary, the finite diameter primary wire, and (c) the 
magnetic permeability the quartz form. Direct measurements Bailey 
small sample taken from the quartz cylinder whose density was 2.1 give 
value for the volume susceptibility. 
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Jones’ formula gives the mutual inductance between helix and circle 
the plane the end the helix, and must applied twice for circle 
displaced from the end the helix. Hence 


Here measured diameter secondary 0.257,600 m., 

axial distances from secondary circle m., 

pitch 0.000,999,654, 


are complete elliptic integrals the first, second, and third kind, 
the modulus and parameter where and are evaluated from the dimen- 
sions by: 

Bartky (2) gave method for reducing such integrals arithmetico- 
geometrical series (already partly established King) and the Bureau 
Standards (16) applied the method the computation similar mutual 
inductance, giving very complete examples. 

The computation yields 


Doubling this because there are two primary coils and multiplying 460, the 
original number turns the secondary, gives 
0.010,029,264,8 


One turn less will obviously bring the result closer mh., but complete 
calculation required for each primary coil account for the removal one 
turn the secondary because this turn displaced from the neutral axis and 
the electrical system longer symmetrical. The amount allow for the 
removal the 260th turn found h., value which 
smaller than simple proportionality would indicate. The true 
value was checked means second order Taylor expansion with respect 
the radius and position the 260th turn. 

Hence our inductance becomes 0.010,007,463,2 with the following 
corrections: 

(A) cross-section secondary h., 


(B) cross-section primary wire 
(C) magnetic susceptibility quartz 


The first these (A) briefly summarized below, the other two are taken 
from the accompanying article (15). (B) represents the effect due non-uni- 
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form distribution current the primary wire; (C) represents the allowance 
made for the difference magnetic permeability vacuum and the 
material the primary former. 

therefore conclude that the value the N.R.C. mutual inductor 


henry. 


From the estimates the errors (diameter primary and axial 
positions the coils) can arrive order magnitude the resulting 
error the following simple analogy. Consider two coaxial circles 
another. Their mutual inductance and 


therefore seems reasonable believe that the probable error smaller 
than parts per million from errors the primary helices. 

The radius the secondary (A) will probably not affect the final result but 
errors the dimensions the cross-section the secondary coil might pro- 
duce error large three parts per million. 

Summing these linearly conclude that the error definitely smaller 
than p.p.m. emphasized that obtaining this number not only 
assumed that circumstances are most unfavorable and all errors act the 
same direction but also assumed that each individual error has been 
assigned relatively large value. Hence p.p.m. felt very definitely 
upper limit rather than probable error. 

The effect temperature the value this mutual inductance very 
small. Following similar reasoning that Dye and Hartshorn (9) the 
temperature coefficient would the order few parts per degree 
The temperature the room now controlled +1°C. and separate 
enclosure provided for the inductor that closer control its temperature 
can arranged. 

Direct comparison this standard with those other laboratories would 
difficult and for the present will depend the comparison between labora- 
tories resistances, after the unit resistance has been determined (14) 
terms the calculated inductance. 


FINITE CROSS-SECTION THE SECONDARY 


(4) and (8) earlier work the effect finite cross-section 
the secondary was evaluated using formula due Searle. Wilmotte 
(17) further reviewed the problem and obtained equivalent solution. 

Fraser this laboratory, starting with the assumption that the secondary 
composed group concentric circles different radii and different axial 
positions and number equal the number turns the secondary, 
arrived formula expressed terms the number turns each 
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layer the secondary and the number layers the winding. The total 
correction for the finite size the secondary given 


result that good the fourth order small quantities for symmetrical 
winding, and being the radial depth between centers the wires, the 
inner and outer layers, and the axial width between centers the exterior 
wires. 

This formula gives equivalent results those derived from 
Applying Fraser’s formula the present case, where 


and 
find that 


which represents relative value 15.5 p.p.m. 

This correction can error depending errors measurement the 
radial (6H) and axial dimensions the secondary coil. these are 
large +0.05 mm. the resulting error can shown large 
2.7 p.p.m. The dimensions should somewhat better known than this assump- 
tion that correction the value the mutual inductance should smaller 
than p.p.m. used above. 
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CONSIDERATIONS THE PRIMARY WINDING 
THE CAMPBELL STANDARD 


ABSTRACT 


The effect the presence the quartz cylinder the mutual inductance 
the standard examined using the expression for the energy the 
circuit made the primary and the secondary connected series. The analysis 
and the computation show that the volume susceptibility the quartz, equal 
(in m.k.s. rationalized units), produces decrease the 
mutual inductance the order parts the second section the effect 
the finite cross-section the primary wire examined. The assumption that 
the current concentrated the central filament the wire was several past 
cases considered sufficient. The computation given here shows that this effect 
not negligible and the actual N.R.C. standard amounts parts 


EFFECT THE MAGNETIC SUSCEPTIBILITY QUARTZ 


The energy circuit placed medium magnetic permeability 
and carrying current given, the m.k.s. rationalized system units, 
the following expression 


where the self-inductance the circuit, the magnetic field produced 
the current, and the magnetic induction given 
[2] 


Let now consider circuit composed two practically closed loops and 
interconnected two parallel conductors (actually they are twisted 
concentric). Considering negligible the effects these wires, and following 
the method previously used for self-inductance standards (1) can write 


Fic. but practically closed, loops. 
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and designating the self-inductance loops and and their mutual 
inductance. Resolving the total magnetic field into partial components 
created the loops and have 

[4] ” 
may written the following form: 


W uf Jor, B’) dv+ Jar. B”) dv+ Jor, B”) dv 


Comparing [3] and [5] and noticing that the third and the fourth integrals are 
equal, have 


portion space has permeability different from then 


designating the space outside v’. 
Assuming that the permeability free space (or practically air), the 
relation between and 


being the volume susceptibility the medium v’. Then 


The second integral considered due the susceptibility 
the volume v’, 


This expression easily extended the case several media v’, etc. 

will now establish, for Campbell type mutual inductor, the expression 
for magnetic field which would suitable for numerical calculations. 

Let designate one the primary helices turns wound 
cylinder radius (Fig. 2). The magnetic induction (or the field 
produced will the same all points circle radius (the helix 
being considered equivalent current sheet). The mutual inductance 
between and given the following integral: 


| 
| 
| 
| 
| 
| 
| 
| 
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Fic. position one the helices (C’) and the secondary the Campbell 
inductor. auxiliary circle used compute the magnetic field any point the space. 


being the current the helix and the vector potential produced 
any point Because cylindrical symmetry the integration above can 
performed and yields 


designating the scalar magnitude the vector potential. From this finally 
the field calculated the relation 


coordinates may choose and being the distance between and 

then 

being expressed terms and its derivatives. The mutual inductance 
and 
[15] 


the Jones formula [16a] distances and are considered algebraical 
quantities, positive direction being (arbitrarily) chosen the inductor’s 
axis Fig. The formula Jones and its derivatives apply all positions 
either outside inside the helix. 


[16a] 


with 


c' 
z 
| 
| 
{ 
| 
| 
| 
i 
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2 2 


Combining [14] with [15] and [16] find that 


For the calculation [17] the Nagaoka and Sakurai tables are 
sufficient. note them given the end this section. 

actual Campbell inductor the primary made two coaxial helices, 
nominally identical. The field produced one them computed the 
formulae [17] and [18]. The total field, i.e. the field produced both helices, 

The secondary Campbell inductor coil small square cross-section 
(11 mm. X11 mm.) containing 459 turns wire. The actual calculation 
was made for one single turn lying the neutral plane the inductor. Calling 
the radius the secondary and the distance from its plane the plane 
have apply the following formulae where designates the mutual 
inductance between and 


From these formulae and using the relation being the 
vector potential created any point the circle p), the values field 


af m ) 1 [ ] 
With these expressions form, cylindrical coordinates, the expression 
for AM: 


The volume made the volume the quartz cylinder which the 
primaries are wound and the volume the marble ring containing the 
secondary coil. The latter lying the neutral circle primary helices, 


} 
| 
4 
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reasonable expect that the magnetic field the vicinity the secondary 
will very this marble ring relatively far from the primaries, the 
portion space filled with marble was entirely neglected our numerical 
calculation, and the integral was extended only the volume occupied 
the quartz cylinder. The dimensions the quartz cylinder are given Fig. 


Fic. Cross-section the quartz cylinder showing points such (x, where the magnetic 
field was actually computed (dimensions cm.). 


(a) Field produced the primary windings.—The radial component H,’ and 
the axial component H,’ produced one helix were computed formulae 
[17] and [18] using the following set values and 

0.1125 (inside the quartz, close the inner surface), 


0.1252, 0.1379, 0.1475 m., 
0.1507 (inside the quartz, close the outer surface), 
0.15, 0.16, 0.17, 0.18, 0.19, 0.20, 0.21, 0.22, 0.25, 0.30, 0.35, 0.40, 


0.50, 0.60 


The results these computations for the radial and axial components the 
total field are represented Figs. and 

(b) Field produced the secondary coil.—This field being much smaller 
than the primary field, was found unnecessary take many values 
and extend the computation beyond 0.30 Results are given Figs. 


and 

The double integration was made various preliminary operations 
graphical evaluation areas, and the final computation also means 
planimeter. The result was amp.): 


i , ” 1.31X10° 


The volume susceptibility® had been determined experimentally means 


the present application the standard the current the secondary equal zero. has 


effect the mutual inductance. 
5In tables actually existing the susceptibility given the non-rationalized system. such 


system would equal i.e. times smaller than the rationalized system. 


| 
j 
| 
| 
] 
| 
| 
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HELIX 


Fic. Radial component the magnetic field produced primary helices. 
Fic. Axial component the magnetic field produced both primary helices. 


magnetic balance sample quartz taken from the cylinder and had 
been found equal 
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Ampere 
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HELIX 


Fic. Radial component the magnetic field produced the secondary winding. 
Fic. Axial component magnetic field produced the secondary winding. 
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The mutual inductance between the primary (both helices) and the second- 
ary Sis The relative value correction therefore: 


The magnetic susceptibility quartz decreases the mutual inductance 2.34 
parts million. 
Note Nagaoka and Sakurai Tables 


these tables (2) the quantity given, being the mutual induct- 
ance circles radii and (Table No. 140). This quantity func- 
tion defined [16]. and tables are the classical c.g.s. system, 


being simple number, independent units. this paper, the similar 
relation takes the form 


therefore 


Now the formula [17] can put into the form 


a, GeP (ap) pv (2/e)(Q: Q:). 


The values and are therefore the ones given and tables for 


EFFECT THE CROSS SECTION THE PRIMARY WIRE 


The computation the mutual inductance Campbell type inductor 
means the Jones formula based the assumption that the primary 
current concentrated the central filament the wire. The error arising 
from this assumption relatively small, and this due the fact that the 
distance between the primary helix and the secondary coil very large 
compared with the cross-section the wire. Since Moon’s method was applied 
the construction mutual inductors the quality primary windings notably 
improved that the calculation these inductors the distribution current 
over the cross-section can longer neglected. Two types current dis- 
tribution may considered: 

(1) Ina perfectly homogeneous ring (Fig. the local current density 
point its cross section may taken inversely proportional the 
distance between and the axis the ring: 


This expression rigorously valid only for true ring (toroid), but con- 
sidered good approximation for helical conductor small pitch. Such 
distribution generally called ‘‘natural’’. 


: | 
| 
| 
i 
| 
| 
| 
| 
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Fic. Diagram representing the cross-section circular ring, y-axis being its axis 
symmetry. 


(2) ring made originally straight wire bent into circular loop, 
elastic stresses should also taken into account. Inner filaments are com- 
pressed resistivity, while the outer filaments are stretched. extensive 
investigation was made the National Bureau Standards the influence 
elastic stresses the resistivity wires (3). development works 
Bridgeman and Kotter. The result all these studies that 
can accept the assumption that the resistivity filament proportional 
its length, and thus readily establish the relation 


being constant dependent only the total current the ring. can 
computed the integral 


The computation mutual inductance between circle and ring finite 
cross-section accomplished two steps. 

and short current sheet radius and length (Fig. 9). 


Fic. Relative position the secondary and the current sheet C2, the latter being con- 
sidered composed infinitesimal sheets dy. 


i 
} 
| 
-- 

| 
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designates the total current the sheet, the current differential sheet 
will equal 


Applying the Neumann formula for the energy between and find the 
following expression: 


The integration with respect can performed and the remaining expres- 
sion can reduced the classical formula Jones. 
Let use the notation: 


gives 


and 


and, therefore, 


Terms higher order are neglected. 

The expression for energy becomes 

2 
0 

and can written 
[29] w= = 
The symbol designates the mutual inductance between (circle radius 
and the current sheet designates the mutual inductance between 
and the central filament (in which the total current would concen- 


trated). 
evaluate the last term differentiate twice f(y) and find 
with 


Substituting [28] obtain 
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Both integrals can found Bierens Haan Tables (p. 91, Table 58, No. 
and 93, Table 59, No. 4): 


They yield the following expression for suitable for numerical computa- 
tions: 


[32] 


where, usual, designates the quantity 

Complete numerical computations will found the end this note. 
may readily seen, however, that the main term and only small 
second-order correction. For inductors like those used primary standards 
the correction term negligible. This will simplify considerably the 
following computation the effect the radial distribution current. 

Let consider (Fig. 10) circular ring mean radius the cross-section 
being circle (small) radius and carrying the total current The ring 


Fic. showing how the effect the current distributed over the cross-section 
the ring can reduced the effect current concentrated the diameter TU. 


can subdivided into infinitely thin coaxial sheets radii The current 
one these sheets thickness and length 


equal 


= 

j 

j 
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considered function the radius the sheet. Using [29] have 
the mutual energy between and the sheet being their mutual induct- 
ance): 

[35] AW = (a4 Aiz)m. 


Let now assume that flows sheet the same length but radius 
The mutual energy AW, this case would equal 


The relation between and AW, can expressed the series 
[37] AW = indi] (—a)( ) + (—a)*( )+ 


which only first and second order terms (x—a) will retained. 
Integrating between points and gives 


Using [34] and its integral 


obtain 


a+p 


a—p 


Designating (as usual tables integrals) the trinomial 


and also introducing the notation 
[41] W = Wottite + A.M, 
have 


at+p 
a—p 


at+p at+p 


a+p 


| 
[36] = Atiz2)mo. 
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both cases, (1) and (2), all integrals involved may performed using 
Peirce’s Dwight’s Tables. the case when 


obtain after few elementary transformations 


Both brackets are obviously small and, therefore, only approximate value 
for required for correct computation A,M and will assume 
for this purpose that [43] may taken equal (mutual inductance 
between coaxial circles radii and being the distance between their 


[43] 


The first derivative mo’ with respect well known: 


with 


easily computed with sufficient accuracy means Nagaoka and 
Sakurai tables for all necessary values 

The second term being order magnitude smaller than A,M, 
there need establish elaborate formula for This second 
derivative was, therefore, computed numerically for three values namely 
10, 15, and cm., and all intermediate values were found graphical 
interpolation. 

the last operation this analysis, will examine the integration [38] 
the expression which may written follows: 


[46] = 1, my = 1, Ate mo +1, Als my’. 

easy notice that my’, given [44], independent and therefore 
a+p at+p 

a—p a—p 


The first term obviously equal and represents the mutual energy 
between two coaxial circles radii and located distance from one 
another, the circle carrying the total current the ring. 

The numerical computation the second term using any well-known in- 
tegration method (e.g. Simpson’s) would tedious work, the result which 
would add only very little the actual value Wo. The value 
obviously maximum for infinitesimal sheets having the largest length 
i.e. (sheet through the center cross section). If, for example, 


| 
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0.145 

0.000,000,04 uh. 


The error [46], when neglected with respect my’, the 
order 0.3 parts million. the other hand, decreases zero when 
approaches a—p that 0.3 p.p.m. appears the upper limit 
the error The second term [47] was therefore simply assumed 
negligible. 
The computation was performed using the following dimensions the pri- 

mary inductor described above: 

0.257,69 m., 

0.150,98 m., 

0.0003 m., 


[48] 


each helix being treated set 100 rings regularly spaced between 0.10 
and 0.20 The computation was simplified forming groups turns, 
associating turns each group with the mean value and 
Actually was computed for 0.10, 0.20 while 
was computed only for 0.10, 0.15, 0.20 Graphical interpolation was 
used for intermediate values. 


(in m.) Mean for group Mean for group 
0.11 17.781 17.1170 0.87 0.805 
0.12 16.453 15.8795 0.74 0.690 
0.13 15.306 14.7615 0.64 0.600 
0.14 14.217 13.7185 0.56 0.535 
0.15 13.220 12.7635 0.51 0.480 
0.16 12.307 0.45 0.430 
0.17 11.472 11.6890 0.41 0.385 
0.18 10.706 10.3555 0.36 0.340 
0.19 10.005 9.6825 0.32 0.310 
0.20 9.360 0.30 


The substitution these values into [43] followed summation (with 
respect turns) gives the following numerical results: 


h., 


With 459 turns the secondary winding the total correction 


mutual inductance 0.010 this represents 3.6 parts million. 


| 
| 


| 
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NOTE 


NOTE EVIDENCE FOR 2.158 MEV. LEVEL 


study the radiations from 24.6-min. Cu®, Nussbaum al. (2) have 
recently observed 0.85+0.03 Mev. y-ray coincidence with the well-known 
1.33 Mev. y-ray transition Ni®, from which they deduce the existence 
2.18 Mev. level They have also observed y-ray 2.13+0.04 Mev. 
which they interpret crossover transition from the proposed 
Mev. level the ground state Ni®. 

report some time ago Fluharty and Deutsch (1) y-ray from the 
decay with energy given 1.63 Mev. 2.23 Mev., and 
intensity about per disintegration, lends support the conclusions 
Nussbaum and co-workers. The result Fluharty and Deutsch was 

became interest see whether the y-ray could detected manner 
which would enable its intensity and energy determined accurately. 
source cobalt consisting cylinder mm. long and mm. diameter, 
irradiated the NRX reactor strength about 150 mc., was placed 
double lens magnetic spectrometer adjusted for line width 0.83% 
momentum half maximum intensity. The source was oriented with its 
axis coincident with the magnetic axis the spectrometer. uranium con- 
verter mgm./cm.? thick and mm. diameter was affixed the front 
the source without interposition any absorbing material. The detector 
consisted crystal, sufficiently thick absorb completely electrons 
energies less than Mev., coupled Dumont 6292 photomultiplier tube 
means lucite light pipe. employing pulse height selection the 
counting circuit, the background counting rate was reduced very low 
level despite the very intense source used. 

The result survey the region above 1.88 Mev. shown Fig. 
The photoelectron line 2.043 Mev. interpreted the line due 
y-ray energy 2.158+0.005 Mev. The energy the line was determined 
comparison with the energy the photoelectron line from the 1.33 Mev. 
y-ray, using the method extrapolated high energy edges. 

While attempt was made observe corresponding line, neverthe- 
less certain that the line observed does not correspond the line 2.064 
Mev. y-ray, since much more intense photoelectron line would then have 
been observed 1.949 Mev. Further, the end point the Compton electron 
spectrum occurs about 1.930 Mev. expected for y-ray energy 2.158 
Mev. 

The intensity the 2.158 Mev. y-ray was found 


*Issued A.E.C.L. No. 
address: Jewish General Hospital, Montreal, Quebec. 
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PHOTOELECTRON LINE FROM Mev. ¥-RAY 
CONVERTER THICK 


2158 Mev. 
END POINT 
COMPTON ELECTRONS 
L930 Mev. 
7800 8000 8200 
CM) 


Fic. secondary electron spectrum above 1.88 Mev. The line interpreted 
the photoelectron line due 2.158+0.005 Mev. y-ray. The vertical bars through the 
experimental points represent standard deviations. 


per disintegration, intensity calibration the spectrometer using the 
1.38 and 2.76 Mev. y-rays from source the form NaCl crystal 
approximately the same dimensions the source, and similarly situated 
relative the converter. 

seems likely that the 2.158 Mev. y-ray observed corresponds the 
Mev. y-ray found Nussbaum al. The latter authors propose 
assignment for the 2.158 Mev. level. this case, the level can fed 
only 0.346 Mev. y-ray proceeding from the 2.504 Mev. level Ni®, 
since feeding 0.656 Mev. cannot account for more than very small 
part the observed intensity per disintegration. |Accordingly 
search was made for 0.346 Mev. y-ray. This y-ray was not detected, and 
upper limit per disintegration was found for its occurrence. Because 
the high intensity the Compton electrons, proved impracticable 
search for the 0.826 Mev. y-ray (corresponding the 0.85+0.03 Mev. y-ray 
found Nussbaum and co-workers). 

should pointed out that the author has not demonstrated beyond 
question that the observed 2.158 Mev. y-ray proceeds from decay 
and the possibility remains that due some unknown radioactive im- 
purity the sample. felt however that this possibility rather remote. 
The source was several months old prior the experiment, and search 
through the literature has failed indicate contaminant, with lifetime 
the order days longer, which could account for the presence the 
2.158 Mev. y-ray. 


pleasure acknowledge the helpful discussions and suggestions made 
Dr. Elliott during the course the work. 
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LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month proof will sent the authors. 


Neutron Scattering and the Frequency Distribution the 
Normal Modes Vanadium Metal* 


Slow neutrons incident crystalline solid may scattered either elastically the crystal 
whole, inelastically with the production annihilation one phonon the crystal, 
inelastically with production and/or annihilation more than one phonon. For monatomic 
crystal composed atoms which scatter incoherently the incoherent part the differential 
cross section (per steradians per atom) for the elastic process (2) 


and the incoherent part the differential cross section for production one phonon angular 
frequency belonging one normal mode (2, 

where gine the incoherent bound scattering cross section, the number atoms the 
crystal, the mass one atom, exp(—2W) the factor, the incoming 
propagation vector the neutron, and the outgoing propagation vector the neutron. 

Equation [2] has been shown apply Debye crystal (2) and cubic crystal (3). 
also applies randomly oriented powder composed non-cubic crystals, provided the 
factor sufficiently close unity. The cross section for annihilation one 
phonon simply [2] multiplied the Boltzmann factor 


800 
VANADIUM 
= 
2 
2 
200 


20° 30° 40° 50° 60° 
ANGULAR POSITION COUNTER 


Fic. 1. Counting rate of 0.0493 ev. neutrons scattered at 123° by vanadium at 84°K., 296°K., 493°K., asa func- 
tion of the setting of the analyzing spectrometer and therefore of the energy. The curves are displaced vertically 
relative to each other. 


Part this work was presented the New York Meeting the American Physical 
Jan. 27-29, 1955 (Phys. Rev. 
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The incoming and outgoing energies the neutron and respectively) are related the 
frequency the phonon The cross section per energy interval dE’ the 
outgoing neutron 


where g(w) the frequency distribution the normal modes without regard direction 
polarization. Measurement the energy distribution incoherently scattered neutrons thus 
provides method for direct measurement g(w) (3). 

Measurement energy distributions neutrons scattered vanadium metal* several 
temperatures have been made. Monoenergetic neutrons mean energy 0.049; ev. were scat- 
tered symmetrically backwards angle 123° slab vanadium metal 0.35 scattering 
mean free paths thickness. The energy the scattered neutrons was measured means 
crystal spectrometer, using aluminum crystal (1). Fig. the counting rates obtained 
are shown function the angular position the counter the analyzing spectrometer. 
The energy corresponding this position also shown. 

the energy the incoming neutrons strong peak found having width corresponding 
the resolution the apparatus. This peak decreases with increasing temperature, anda Debye 
temperature for vanadium 366°+25° can inferred from its temperature dependence 
using [1] and the Debye model expression (4) for the factor. 


(a) 


Ol 02 .03 04 .05 .06 .07 .08 09 O10 
ENERGY (ev.) 


(b) 


Ol 02 .03 .04 .05 06 .07 .08 .09 0.10 
ENERGY (ev.) 


Fic. 2. (a) The points represent the measured partial differential cross section of vanadium metal at 84°K., 
296°K., and 493°K. for 0.0493 ev. neutrons scattered at 123°. Standard deviations are shown for the 84°K. cross 
sections. The solid curves are the one-phonon cross sections computed for a Debye temperature of 348° and folded 
with the resolution function. The dotted curves are the unmodified cross sections. Estimated resolution functions 


are shown at the bottom of the diagram. ¢ 
(6) Measured and calculated partial differential cross sections of 0.0622 ev. neutrons scattered at 95° by vanadium 


metal at 83°K. 


both sides the elastic peak inelastically scattered neutrons are observed whose inten- 
sity increases with increasing temperature. Converting the counting rates Fig. energy 
distribution and correcting for the spectrometer efficiency (1) the points Fig. are ob- 
tained. The cross sections have been put absolute basis using the elastic peaks and 
equation [1]. the calculations account was taken multiple scattering which involved 
interchange zero one phonons. 
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The upper cutoffs the spectra imply Debye temperature 320°+30°. The lower cutoff 
could not obtained because contamination the second order spectrum the analyzing 
spectrometer. Another spectrum was therefore measured using higher initial energy (0.0622 
ev.) and shown Fig. The low energy end the spectrum has cutoff which sharp 
within the resolution the experiment. The location the cutoff implies Debye tem- 
perature 350°+20°. 

Assuming Debye frequency distribution 

wo, 


where the one-phonon spectra were computed for 348°, the weighted mean 
the above values. They are shown Fig. calculated and folded with the estimated 
resolution function the apparatus. Agreement the calculated spectra with the measured 
cross sections reasonable the lowest temperature, becoming worse the higher tempera- 
tures. Rough estimates the cross sections for multiple phonon transfer suggest that the dif- 
ferences probably arise from this process. 

The results these preliminary experiments can thus satisfactorily accounted for 
Debye model for The experiments are being continued under better conditions 
intensity and resolution, and hoped soon obtain the essential features the frequency 
distribution this method. 


1. Brocxnouse, B. N. Phys. Rev. 99: 601. 1955. 

2. Casse.s, J. M. Progress in nuclear physics. Academic Press, New York. 1950. p. 185. 
3. Piaczex, G. and VAN Hove, L. Phys. Rev. 93: 1207. 1954. 

4. Weinstock, R. Phys. Rev. 65: 8. 1944. 
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CORRECTIONS 

Vol. 29, 1951. Inthe paper The magnetic dipole over the horizontally stratified earth James 
Wait, pp. 577-592, 
(1) the denominator Equation (38) should read 
(2) the sign preceding the integral Equation (52) should sign. 

Vol. 32, 1954. Calculated diffraction patterns dielectric rods centimetric wavelengths 
Froese and James Wait, pp. 775-781, the relation following equation 
should read 
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